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FOREWORD 


During  the  six  years  that  this  contract  with  the  Air  Force  Office  of 
Scientific  Research  has  heen  in  effect,  some  spectacular  developments  in 
quantum  electronics  have  taken  place.  In  195^,  while  the  first  draft  of  the 
technical  proposal  was  in  preparation,  the  term  "quant\im  electronics"  had 
not  yet  heen  coined.  The  stated  purpose  of  the  proposed  program  was  merely 
to  carry  out  some  basic  investigations  of  chemical  and  lattice  defects  in 
crystalline  solids,  using  electron  spin  resonance  as  the  tool. 

By  the  end  of  1956,  however,  Bloemhergen' s  suggestion  of  the  solid-state 
maser  became  very  widely  known,  and  a  short  time  later  the  feasibility  of  the 
suggestion  was  reported  by  Scovil,  Feher,  and  Seidel  of  the  Bell  Telephone 
Laboratories.  At  the  Lincoln  Laboratory,  McWhorter  and  Meyer  discovered  maser 
action  in  chrome  cyanide  around  the  middle  of  1957^  an<i  on  December  20,  195T> 
the  workers  at  The  University  of  Michigan  Willow  Run  Laboratories  observed 
maser  action  in  pink  ruby. 

Development  of  the  ruby  maser  as  a  device  took  place  very  swiftly.  Dur¬ 
ing  these  months  after  the  discovery  of  maser  action  in  ruby,  Townes  and 
Giordemaine,  using  a  tested  pink  ruby  crystal  that  had  been  supplied  by  The 
University  of  Michigan  workers,  announced  the  successful  operation  of  the 
radiotelescope  maser  at  the  Naval  Research  Laboratory.  A  similar  device  was 
built  later  for  The  University  of  Michigan  radiotelescope.  More  recently, 
the  ruby  maser  was  used  by  the  workers  at  the  Jet  Propulsion  Laboratory  for 
the  precision  determination  of  the  Astronomical  Unit  (149,589,500  Km).  At 

xi 


the  same  laboratory,  a  program  evaluating  the  use  of  ruby  masers  for  satellite 


communication  is  in  progress,  and  it  is  well  known,  of  course,  that  a  ruby 
maser  is  part  of  the  Telstar  communication  system. 

The  pink  ruby,  which  previously  had  been  used  primarily  for  phonograph 
needles,  was  very  soon  found  to  have  another  important  application  in  a  de¬ 
vice  now  popularly  called  the  laser.  Soon  after  theoretical  suggestions  were 
offered  by  Schalow  and  Townes,  T.  H.  Maiman  of  the  Hughes  Research  Laboratory 
and  Collins  and  his  co-workers  at  the  Bell  Telephone  Laboratories  reported 
the  predicted  phenomenon. 

The  program  carried  out  under  this  AFOSR  contract  is  a  classical  example 
of  the  importance  of  basic  research  in  technological  developments.  As  men¬ 
tioned  earlier,  our  initial  objective  was  merely  to  explore  the  uses  of  elec¬ 
tron  spin  resonance  as  a  tool  to  study  defects  in  hard  crystalline  materials, 
to  which  only  limited  attention  had  been  paid  at  that  time.  It  was  this  pro¬ 
gram  that  led  this  principal  Investigator  to  examine  the  merits  of  such  mate¬ 
rials  as  ZnS,  MsO,  CaCoa,  and  AI2C3  as  maser  materials,  and  thus  to  the  ob¬ 
servation  of  maser  action  pink  ruby. 

The  technological  developments  in  masers  and  lasers  have  in  turn  already 
raised  a  number  of  questions,  the  answers  to  which  can  come  only  from  more 
intensive  basic  investigations.  One  such  question  is  the  very  old  one  of  en¬ 
ergy  transfer  and  transformation  mechanisms  in  solids.  We  have  noticed  that 
the  mechanism  invoked  to  explain  laser  action  in  solids  is  different  from 
that  needed  for  luminescence  and  scintillation  phenomena,  A  program  to  study 
these  phenomena  has  been  started,  with  the  hope  of  gaining  deeper  insight 


xli 


to  energy  conversion  processes  in  solids. 


In  "bringing  this  contract  to  a  close,  I  wish  to  thank  Mr.  Charles  F. 
Yost,  now  of  the  Advanced  Research  Projects  Agency,  and  the  various  members 
of  the  Air  Force  Office  of  Scientific  Research  for  the  generous  support  that 
made  our  work  on  the  ruby  maser  at  Willow  Run  possible,  and  for  the  continua¬ 
tion  of  the  program  to  introduce  these  solid-state  concepts  and  techniques 
into  the  research  programs  of  the  Department  of  Nuclear  Engineering. 

Chihiro  Kikuchi 
December,  I962 
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PROPERTIES  OF  RUBY  AND  OTHER  SAPPHIRES 


C.  Kikuchi  and  S.  Karavelas 


I.  INTRODUCTION 


For  centuries  ruby  has  been  highly  prized  as  a  gem,  but  It  is  only  very 
recently,  since  the  discovery  of  maser  action  in  ruby,  that  the  industrial 
importance  of  this  material  has  come  to  be  appreciated.  Before  1957>  "the 
Linde  Company  was  the  principal  supplier  of  ruby  and  sapphires.  The  so- 
called  pink  ruby  was  manufactured  by  this  company  to  make  the  so-called  sap¬ 
phire  phonograph  needles. 

The  present  surge  of  Interest  in  ruby  and  related  materials  stems  from 
the  fact  that  ruby  has  been  shown  to  be  useful  as  a  maser  and  laser  material. 
Maser,  an  acronym  for  Microwave  Amplification  by  Stimulated  Emission  of  Radia¬ 
tion,  is  a  low-noise,  high-gain  amplifier,  for  which  ruby  is  the  critical 
material.  The  ruby  maser  was  used  in  the  Telstar  communication  system,^  and 
at  the  California  Technology  Jet  Propulsion  Laboratory  a  packaged  ruby  maser 
for  satellite  tracking  stations  is  under  evaluation.  In  scientific  applica¬ 
tion,  the  ruby  maser  made  possible  the  precision  measurement  of  the  Astro¬ 
nomical  Unit,  the  average  distance  of  the  earth  from  the  sun.  The  precision 
value  of  li+9,589,500  ±  500  km  was  obtained  by  the  workers  at  the  Jet  Propul- 
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sion  Laboratory,  after  making  careful  analysis  of  the  radar  echo  from  Venus. 

Ruby  is  also  used  in  lasers,  a  device  for  the  generation  of  intense 
coherent  optical  radiations.  The  word  Laser  is  an  acronym  for  Light  Ampli¬ 
fication  by  Stimulated  Elnission  of  Radiation.  The  industrial  and  techno¬ 
logical  uses  of  this  device  are  still  under  research  and  development.  It  is 
possible  that  the  device  will  find  applications  in  communication  systems 
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operating  at  optical  frequencies,  and  in  certain  Industrial  processes  such 


as  micro-cutting,  micro-etching,  and  micro -welding.  It  also  seems  to  be  use- 
full  as  a  medical  tool— in  coagulating  a  detached  retina  onto  the  eyeball, 
for  example. 

A  great  deal  has  been  said  about  the  usefulness  of  ruby  as  a  maser  and 
laser  material,  but  few  popular  expositions  attempt  to  explain  why  ruby  is 
useful.  Consequently,  it  will  be  the  purpose  of  this  discussion  to  point 
out  the  factors  that  make  ruby  behave  the  way  it  does. 
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II.  CRYSTAL  STRUCTURE 


Chemically,  ruby  is  aliaminum  oxide  which  contains  a  small  concentration 
of  chromium.  Aluminum  oxide  occurs  in  two  forms;  the  one  that  concerns  us  is 
known  as  a-AlgOa.  In  mineralogy,  this  form  is  called  corundum,  and  the  com¬ 
mercial  name  of  synthetic  corundum  is  sapphire.  The  commerfelal  sapphire, 
which  is  clear,  transparent,  and  colorless,  should  not  he  confused  with  the 
gem  sapphire,  which  is  blue.  (In  passing,  perhaps  it  should  be  noted  that 
the  technological  importance  of  sapphire  is  beginning  to  be  realized.  For 
exaii5)le,  according  to  a  recent  NASA  report,  sapphire  windows  will  henceforth 
be  used  on  satellite  solar  batteries.  Also,  sapphire  has  the  very  unusual 
property  of  high  thermal  conductivity  but  low  electrical  conductivity,  in 
contradiction  to  the  Weldemann-Franz  law,  which  asserts  that  good  electrical 
conductors  are  also  good  thermal  conductors.  This  particular  property  is 
used  in  the  laboratory  to  provide  good  thermal  contact  euid  good  electrical 
Insulation. ) 

There  are  several  varieties  of  rubles,  distinguished  by  their  colors. 

Gem  rubies  are  deep  red,  due  to  the  high  chromium  concentration,  whereas  the 
rubies  Important  for  masers  and  lasers  are  pink,  due  to  the  chromium  con¬ 
centration  of  0.1^  or  less. 

5  h 

G?he  crystal  structure  of  sapphire  is  rather  complicated.  '  It  can  be 
conveniently  generated  by  placing  AI2O3  molecules  at  the  comers  of  a  cube 
and  stretching  the  cube  along  one  of  the  body  diagonals.  Another  AlaOa, 
rotated  I80®  about  the  molecular  axis  with  respect  to  the  first  molecules,  is 
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placed  at  the  center  of  this  distorted  cube.  Upon  careful  examination  of 
this  structure,  it  will  be  seen  that  each  A1  atom  is  sandwiched  between  two 
groups  of  three  oxygens,  as  shown  in  Fig.  1.  The  three  oxygens  are  in  a 


Fig.  1.  Coordination  of  A1  atom  in  a-AlgOa. 

plane  about  1.37^  from  the  A1  site,  and  the  Al-0  distance  is  about  I.98A,  mak¬ 
ing  an  angle  of  1+6  °27'  with  the  crystal  c-axis.  The  other  three  oxygens  lie 
in  a  plane  O.80A  away;  the  Al-0  distance  is  about  1.8UA  and  the  angle  61+ “10'. 
The  relative  orientations  of  the  two  oxygen  triangles  are  not  quite  l80“.  If 
the  oxygens  are  projected  on  a  horizontal  plane,  the  result  shown  in  Fig.  2 
is  obtained. 
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III.  MERGY-LEVEL  DIAGRAM 


Our  next  task  is  to  analyze  the  effect  of  the  six  nearby  oxygens  when 

li+ 

an  impurity  ion  such  as  ,  Cr-^  ,  or  Mn  is  substituted  for  Al.  These 
three  ions  are  cited  because  each  has  the  Jd^  configuration,  so  that  energy- 
level  diagrams  of  all  three  have  many  qualitative  features  in  common.  In 
most  theoretical  analyses,  the  assumption  is  made  that  the  ion  consists  of 
an  inert  inner  ion  core  (the  argon  core)  and  that  the  electronic  properties 
of  the  ions  are  determined  by  the  three  electrons  in  the  outer  unfilled  3d 
shell.  The  validity  of  this  assumption  has  been  questioned  recently,  but  for 
the  present  discussion  we  shall  assume  the  existence  of  such  inert  inner 
filled  electron  shells. 

Earlier,  it  was  emphasized  that  the  crystal  structure  of  sapphire  is 
complex.  To  make  theoretical  analysis  somewhat  tractable,  it  is  generally 
assumed  that  the  oxygens  give  rise  to  a  crystalline  electric  field  which  has 
predominantly  cubic  symmetry  but  which  also  has  a  small  component  of  tri¬ 
gonal  symmetry.  Such  a  crystalline  electric  field  can  arise  if  we  first 
imagine  that  the  oxygens  are  placed  at  the  center  of  the  faces  of  a  cube 
surrounding  the  impurity  ion,  and  then  imagine  a  small  distortion  produced 
by  stretching  the  cube  along  the  body  diagonal.  If  these  assumptions  are 
made,  we  can  first  take  into  account  the  effect  of  the  cubic  component  of 
the  crystalline  electric  field,  and  later  consider  the  perturbation  of  the 
cubic-field  energy  levels  by  the  smell  trigonal  field. 


8 


The  details  of  this  quantum  mechanical  calculation  is  somewhat  involved, 
but  the  general  qualitative  features  of  the  energy-level  diagram  can  be  pre¬ 
dicted  readily  from  group  theory.  We  are  concerned  with  what  is  called  the 
crystal-field  spectra,  i.e.,  the  spectrum  which  arises  from  electrons  making 
transitions  within  the  Jd  shell.  For  a  cubic  field,  the  energy-level  diagram 
would  be  as  shown  in  Fig.  3:^ 


2T2 


Fig.  3-  Energy-level  diagram  for  cubic  field. 

The  niimerlcal  values  given  for  the  wavelengths  of  the  emitted  and/or  absorbed 
radiation  are  for  ruby.  For  the  other  two  isoelectric  ions,  V^"*"  and  Mn^'*’, 
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the  wavelengths  would  he  slightly  different. 

If  we  next  take  the  trigonal  C2rystalline  electric  field  into  account, 
each  of  the  indicated  levels  will  split  into  two  or  more  levels.  The 
level  will  split  into  two  levels;  the  separation  in  ruhy,  for  example,  is 
about  29  cm“^,  as  shown  in  Fig.  V. 


Fig.  k.  Split  in  level  by  trigonal  field. 

This  results  in  the  splitting  of  the  R  line  into  two  components,  referred  to 
in  the  literature  as  the  Ri  and  R2  lines.  The  former  is  the  ruby  laser  line. 
The  ground  state  also  splits,  as  Indicated  in  Fig.  5- 


Fig.  5*  Split  in  ground  level. 
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The  separation  of  the  levels  in  this  case  is  very  small,  only  about  l/lOO 
that  of  the  state.  But  this  very  small  splitting  is  Important  in  the 
ruhy  maser. 
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IV.  PROPERTIES  OF  OTHER  IONS 

In  the  preceding  section  a  few  of  the  properties  of  Cr^  in  sapphire 
were  described.  The  next  question  that  arises  is  whether  or  not  it  is  pos¬ 
sible  to  fabricate  other  materials  which  have  similar  properties.  To  answer 

o+  li+ 

this  question,  a  brief  discussion  of  v  and  Mn  will  be  presented.  The 
results  for  Mn^"*"  were  published  very  recently  by  S.  Geschwind  and  others  at 
the  Bell  Telephone  Laboratories."^  These  results  along  with  properties  for 
V^'*’  and  Cr^^  are  summarized  in  Fig.  6. 


t »  p  t“  3.4  ms  t*  .8ms 

Fig.  6.  Properties  of  V^^,  Cr^^,  and  Mn^^. 
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The  measured  values  of  Mn^'''  are; 

Ra  =  14,866  cm-i  (79°K) 

Ri  =  14,786 
2D  =  -.391^ 

g||  ~gj_  =  1-9937 

|A|i  I  =  IAj^I  =  70.0  X  10-4  cm"^ 

There  are  several  interesting  facts  to  note.  In  the  first  place,  the  R^  lines 
are  very  nearly  equal.  Another  very  remarkable  fact  is  that  the  results  for 
Mn^"*"  in  sapphire  are  very  close  to  those  of  Mn^'*'  in  lithium  titanate.  In  a 
paper  by  Lorenz  and  Prener®  the  values  k  =  6790^  (14,700  cm"^)  and  t  =  1.1  ms 
are  reported.  The  values  for  the  ground-state  splitting  are  not  kno\m,  be¬ 
cause  the  measurements  were  made  on  a  powder.  The  very  close  agreement  of 
the  numerical  values  for  sapphire  and  lithium  titanate,  despite  the  apparent 
differences  in  the  host  crystal  material,  stems  from  the  fact  that  coordina¬ 
tion  for  Mn^^  is  octahedral;  i.e.,  in  both  cases  the  Mn^"*"  ion  is  bonded  to 
six  nearby  oxygens. 

The  work  on  in  sapphire  was  started  in  our  laboratory  when  it  was 

noticed  that  is  produced  by  X-rays. 5'  Our  experiments  showed  that  normally 

vanadlijm  in  sapphire  is  present  as  but  that  after  irradiation,  part  of 

this  converted  to  The  chemical  impurity  responsible  for  the  stablllza- 

■2+  4+ 

tion  of  V  is  not  yet  known.  For  Mn  in  both  sapphire  and  lithium  tetanate, 
the  presence  of  Mg  impurity  is  needed. 

Vanadium  is  a  favorable  system  to  investigate  because  it  is  essentially 
lOO'jt  and  because  its  nuclear  magnetic  moment  is  rather  large,  due  to 
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the  fact  that  vanadium  has  an  odd-proton-even-neutron  nucleus.  The  nuclear 
magnetic  moment  can  then  be  used  as  a  signature  to  identify  vanadium,  and  for 
this  reason  vanadium  can  be  quite  easily  identified  in  the  different  oxida¬ 
tion  states  in  sapphire. 

With  these  preliminary  remarks,  let  us  return  to  the  ruby  energy-level 

diagrams  given  in  Section  III,  As  mentioned  before,  the  paramagnetic  ion  in 

ruby  is  Cr^  ,  whose  electron  shell  structure  is  15‘^25^2p°35"^3p  The  three 

3d  electrons  in  the  last  unfilled  shell  are  the  electrons  responsible  for  the 

Interesting  microwave  and  optical  properties  of  Cr^'*’.  (it  should  be  noted 
9+  I4.+ 

that  V  and  Mn  have  the  same  electron  shell  structure,  so  that  the  com- 

0+  3+  i|.+ 

ments  about  to  be  made  are  applicable  to  the  three  ions  V  ,  Cr^  ,  and  Mn  . ) 

It  was  mentioned  earlier  that  the  chromium  ions,  occupying  the  A1  sub¬ 
stitutional  sites,  are  surrounded  by  six  oxygens.  The  question  we  now  wish 
to  answer  is:  What  effect  will  these  oxygens  have  upon  the  Cr^'*’  energy  levels? 
Since  the  actual  arrangement  of  the  oxygens  is  quite  complex,  to  make  the 
analysis  more  tractable  we  shall  first  assume  that  the  effects  of  the  oxygens 
come  predominantly  from  the  octahedral  arrangement  of  the  oxygens.  In  other 
words,  for  the  first  step  in  the  analysis  we  shall  assinne  that,  with  the  Cr^'*’ 
ion  at  the  origin  of  the  coordinate  system,  the  oxygens  are  located  along  the 
coordinate  axes  at  the  same  distance  from  the  origin.  Later  we  shall  con¬ 
sider  the  effects  of  a  small  distortion  along  the  body  diagonal. 

Perhaps  the  graphical  representation  in  Fig.  7  will  make  the  physical 


lU 


problem  clearer. 


Fig.  7.  Octaiiedral  coordination. 

Electrons  1,  2,  and  5  are  moving  in  the  spherically  symmetric  coulomb 
field  of  +6  e.  Each  electron  has  orbital  and  spin  angular  momentum  li  and  si 
respectively.  There  are  mutual  coulomb  repulsive  forces  among  the  electrons, 
and  the  negative  charges  located  along  the  coordinate  axes  will  also  have  an 
important  influence  on  the  motion  of  the  electrons. 

The  problem  indicated  here  is  a  complicated  one;  consequently  it  will  be 
necessary  for  us  to  make  a  series  of  approximations.  In  order  to  avoid  the 
complexities  of  theoretical  arguments,  we  shall  first  consider  the  case  of  a 
single  5d  electron,  and  shall  then  show  how  additional  factors  have  to  be 
brought  into  the  analyses  as  we  proceed  to  the  5d  and  the  Jd^  cases. 

The  case  of  the  single  5^  electron  will  be  dealt  with  in  great  detail  in 
order  to  bring  out  many  features  of  the  theoretical  techniques.  The  single 
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3d  electron  is  in  the  central  spherical  syiranetric  coulomb  field,  but  its  orb¬ 


ital  motion  is  profoundly  affected  by  the  oxygens  surrounding  it.  Further¬ 
more,  the  electron  has  orbital  and  spin  angular  momentum,  so  that  the  coup¬ 
ling  of  these  vectors  will  affect  the  electron  energy  level.  The  Schroedinger 
equation  we  need  to  solve  is 

H\|f  =  (4-1) 

with 


s  -  —  +  V(r)  +  Fc(x^+y‘*^+z'^-  -  r"^)  +  F^(xy+yz+zx)  +  \^  •  s  (4-2a) 

10 

The  full  trigonal  part  of  the  Hamiltonian  up  to  fourth  power  is 


Vtr  =  Ft(xy+yz+zx)  +  F^ [xyz(x+y+z)  -  y  r^(xy+yz+zx)  ] 

+  F^(x+y+z)  +  F^'  xyz  +  F.|.  [x^+y^+z^-  ^  r^(x+y+z)  ]  (4-2b) 

J"*"  3"*‘ 

if  the  approximation  is  made  that  the  symmetry  of  the  A1  or  Cr  site  is 

Cav  This  is  equivalent  to  neglecting  the  angle  of  4'’22'  shown  in  Fig.  2. 

In  the  following  discussion  we  restrict  ourselves  to  the  form  (4-2a)  for  the 

sake  of  simplicity.  However,  one  important  characteristic  of  (4-2b)  must  be 

kept  in  mind,  namely,  that  it  has  no  center  of  inversion.  This  is  important 

when  selection  rules  for  electric  dipole  transitions  are  examined. 

There  are  three  steps  to  the  Schroedinger  equation,  as  follows: 

I.  -  +  V(r)  +  FcCx'^+y^+z^-  |  r^) 

2m  ^ 

2  3 

II.  ■  ^  ^  +  V(r)  +  Fc(x^+y'*+z'*-  “  r"*)  +  F^(xy+yz+zx) 

2  3 

III.  ■  ^  ^  +  V(r)  +  Fp(x'*+y'*+z'*-  r'*)  +  F^(xy+yz+zx)  +  •  s 
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The  purpose  of  Step  I  will  he  to  show  how  the  symmetry  properties  of  the 


Hamiltonian  can  be  exploited  to  construct  wave  functions,  and  also  to  point 
out  the  physical  reasons  for  the  importance  of  group  theory  in  attacking  these 
problems.  Consider  then  the  first  three  terms  in  (4-2a),  i.e., 

1^0  "  ■  SS  +  V(r)  +  Fc(x4+y4+z^-  ^  r4)  (4-5) 

We  notice  that  the  first  two  terms  are  invariant  under  all  rotations,  finite 
or  infinitesimal.  Consequently,  as  is  well  known,  if  the  last  term  were 
absent,  the  solution  of  the  Schroedinger  equation  could  be  represented  by  the 
spherical  hamonlcs.  The  corresponding  energy  levels  are  called  the  s,  p,  d, 
etc.,  states.  The  mathematical  significance  of  this  remark  is  that  the  wave 
functions  belonging  to  a  particular  energy  level  can  be  represented  by  poly¬ 
nomials  of  X,  y,  and  z  which  transform  among  themselves  under  infinitesimal 
rotations.  For  example,  the  wave  functions  for  the  d-electron  can  be  repre¬ 
sented  either  by  the  set  of  spherical  harmonics  of  degree  2 — Y2,  m  =  0,  ±1, 

±2 — or  by  the  polynomials  xy,  yz,  zx,  x^-y^,  5z^-r^. 

But  what  happens  if  the  symmetry  of  the  Hamiltonian  is  lower?  For  ex¬ 
ample,  in  (4-5),  the  Hamiltonian  is  Invariant  under  the  simple  substitution 
of  ±x,  ±y,  ±z  among  themselves,  but  not  under  infinitesimal  rotations.  In 
the  latter  case, 

X  — ^  X  -  ey 

y — ^  y  +  ey  (^-4) 

and  polynomials  in  mixed  powers  of  x  and  y  are  generated. 
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However,  the  Hamiltonian  is  invariant  under  a  group  of  transformations 
called  the  cubic  group.  This  is  the  group  of  operations  that  carry  a  cube 
into  itself.  Let  us  see  how  we  can  make  use  of  the  Hamiltonians  invariance 
under  the  cubic  group  to  generate  wave  functions  belonging  to  the  same  energy 
level.  For  this,  we  need  to  recall  that  if  ijr  is  a  wave  function  belonging 
to  the  energy  level  E  of 

Ht  =  E\|;  (4-1) 

and  if  Rg  is  one  of  a  group  of  transformations  that  leaves  the  Hamiltonian  un¬ 
changed,  i.e.. 


Rg(a|f)  =  H(Rgt)  (4-5) 

then  the  set  of  functions  Rg\|r  generated  by  the  group  of  operations  Rg  are 
also  solutions  belonging  to  the  same  energy  level  E.  To  illustrate  the  phys¬ 
ical  significance  of  this  result,  let  us  examine  the  effect  of  the  cubic 
crystalline  electric  field  upon  the  5<i  electrons.  As  mentioned  earlier,  the 
5d  wave  functions  can  be  represented  by  the  five  linearly  independent  mon¬ 
omials  of  the  second  degree:  xy,  yz,  zx,  x^-y^,  and  Consequently, 

let  us  suppose  that 

f(r)xy  (4-6a) 

is  a  solution  of  (4-1)  belonging  to  E.  Since  f(r)  is  a  function  independent 
of  the  angle,  in  order  to  exhibit  the  angular  dependence  hereafter  we  shall 
write  simply 

xy  (4 -6b) 

instead  of  f(r)xy.  According  to  the  comment  made  earlier,  if  Rg  is  an  opera- 
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tion  of  the  group  then  other  solutions  can  be  generated  by  the  operation.  For 
example,  (4-3)  does  not  change  if  the  substitution 

X  - )  X 

y - >  z 

z - > -y  (4-7a) 

is  carried  out.  Geometrically,  this  corresponds  to  a  90°  rotation  about  the 
x-axis.  Carrying  out  this  substitution  on  (4-6b)  we  obtain 

xz  (4-6c) 

Another  possible  substitution  is 


X  - z 

y - )  y 

z - >  X  (4-Tb) 

Geometrically,  this  represents  a  reflection  on  a  diagonal  plane  passing  through 
the  y-axis.  Applied  to  (4-6b)  this  operation  yields 

yz  (4-6d) 

which  along  with  xy  and  zx  belong  to  the  same  energy  level  E.  Also,  it  is 
clear  that  for  the  rest  of  the  substitutions  that  carry  a  cube  into  itself, 
the  functions  xy,  yz,  and  zx,  will  transform  among  themselves,  but  that  at 
no  time  will  they  transform  into  x2-y2  and  3z^-r^.  But  since  these  last  two 
polynomials  will  transform  among  themselves,  they  belong  to  another  energy 
level. 
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V.  TRIGONAL  FIELD  AND  SPIN-ORBIT  COUPLING 


In  the  last  section  it  was  shown  how  symmetry  properties  can  be  used  to 
show  that  the  3<i  electron  will  split  into  two  energy  levels,  with  2-  and  3- 
fold  degeneracy  respectively.  In  order  to  bring  out  other  properties,  we 
shall  consider  the  trigonal  case  in  detail.  A  detailed  discussion  of  this 
case  is  possible  because  there  are  only  six  operations  Involved. 

As  mentioned  earlier,  the  six  oxygens  in  sapphire  form  a  distorted  cube. 
This  can  be  realized  by  distorting  the  cube  along  one  of  the  body  diagonals. 
If  this  is  done  along  the  [ill]  body  diagonal  (see  Fig.  8),  the  additional 
term  F-t(xy+yz+zx)  shown  in  (4-2a)  appears  in  the  Hamiltonian.  What  effect 
will  this  term  have  upon  the  energy  levels?  In  particular,  we  shall  consider 
the  effect  of  such  a  distortion  upon  the  energy  level  associated  with  the 
functions  xy,  yz,  and  zx,  or  some  suitable  linear  combination  of  these  func¬ 
tions  . 


Fig.  8. 


According  to  the  above  comments,  we  need  to  look  for  the  group  of  sub¬ 
stitutions  or  operations  that  leaves  the  Hamiltonian  unchanged.  Clearly  the 
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Hamiltonian 


2m 


V(r)  +  F(,(x4+y4+z4-  ^  r4)  +  F-j-Cxy+yz+zx)  (5-1) 


will  "be  invariant  for  substitutions  that  leave  F^(xy+yz+zx)  unchanged;  and 
these  will  be: 


X 


X 


E: 


(5-2a) 


z 


z 


Z - ^  X 


(5-2b) 


z 


y 


(5-2c) 


Substitution  (5-2a)  is  the  identity  operation;  substitutions  (5-2b)  represent 
120°  and  2U0°  rotations  about  the  [ill]  body  diagonal,  and  substitutions 
(5-2c)  represent  reflections  in  the  three  face  diagonal  planes. 

From  the  set  of  functions  xy,  yz,  and  zx  we  note  that  the  function 

xy  +  yz  +  zx  (5-5) 

can  be  constructed.  Clearly,  this  fiinction  transforms  always  into  Itself 
under  the  operations  indicated  in  (5-2).  Physically  this  means  that  one  of 
the  three  functions  will  split  off  due  to  this  trigonal  field.  To  construct 
the  remaining  two  linearly  independent  functions  consider  the  geometrical 
procedure  shown  in  Fig. 
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Fig.  9*  Geometric  constructions  of  linearly  independent 
functions  (simplified). 

The  plane  normal  to  the  vector  _i+j+k  intersects  the  xy  plane  along  AB,  so  that 
a  vector  in  this  direction  can  be  written 

-i+i 

Referring  to  Fig.  9'b,  this  suggests  the  corresponding  function 

-yz  +  zx  (5-5) 

The  vector  normal  to  -i+j,  and  ^+j+k  is  given  by 


Fig.  10.  Geometric  constructions  of  linearly  independent  functions 


i 

1 

1 


J 

1 

1 


k 

0 

1 


i  +  J  -  2k 


(5-6) 


so  that  the  combination 

yz  +  zx  -  2xy  (5-7) 

is  suggested.  We  shall  now  verify  that  the  functions  of  (5-5)  (5-7)  trans¬ 

form  among  themselves  for  the  substitutions  (5-2)  and  we  shall  furthermore 
construct  the  matrices  representing  the  operations  of  (5-2) .  If  we  let 


fi  s  -yz  +  zx 
fg  s  yz  +  zx  -  2xy 


(5-8) 


then  E  carries  fi  into  fi  and  fg  into  fg.  Consequently  the  representative 
matrix  is 


M(E)  = 


1  0 


0  1 


(5-9) 


Consider  next  Ch  ^ .  In  this  case 


f  1 - >  -zx  +  xy  =  -  I  f  1  -  I  fs 

fa - >zx  +  xy  -  2yz  =  |  f^  -  i  fa 


(5-10) 


so  that  the  matrix  representation  is 


M  C 


(5-11) 
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Similarly, 


fi - ^  -xy+yz  =  -  I  fi  +  I  fa 

fa - ^  xy+yz-2zx  =  -  |  f i  -  |  fs 

1  1 

fi - >  -zy+yx  =  2  ^"1  -  2  ^2 

fa - >  zy+yx-2xz  =  -^fi-^f2 

1  1 

fi - )  -xy+xz  =  -  fi  +  2  fs 

fa - >  xy+xz -2yz  =  |  fi  -  |  fs 


fi - >  -xz+yz  =  -fi 

fa - f  xz+yz-2xy  =  fa 


(5 -12a) 


(5 -12b) 


0^1 

j  (5-12<l) 


Earlier  we  showed  that  when  a  suitable  linear  combination  of  xy,  yz,  and 
zx,  namely 


xy  +  yz  +  zx 


(5-5) 


is  taken,  this  function  is  such  that  it  always  transforms  into  itself  under 
the  operations  of  (5-2).  The  following  question  now  arises:  Is  it  possible 
to  construct  some  linear  combination  of  fi  and  fa  of  (5-8) — say  and  ga — 
such  that  gi  will  always  transform  into  itself  and  ga  into  itself?  The 
answer  to  this  question  is  no;  it  is  provided  by  group  theory  and  stems  from 
the  fact  that  the  set  of  numbers  obtained  by  taking  the  diagonal  sum  of  the 
matrices  is  identical  to  the  row  of  numbers  in  the  group  character  table. 

The  character  table  for  the  Cay  group,  with  which  we  are  concerned  here, 
is  as  follows: 
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The  symbols  appearing  at  the  top  of  the  column  represent  the  symmetry  opera¬ 
tions,  and  the  numbers  in  the  colijmn  below  give  the  trace  or  the  diagonal  sum 
of  the  representative  matrix.  The  numbers  under  E  give  the  dimensionality 
for  the  so-called  irreducible  representations  Ai,  Aa,  and  E.  (The  maxtrix 
representing  the  symmetry  operations  Ai  and  Ag  are  1-x-l  matrices,  whereas 
those  for  E  are  2-x-2  matrices.) 

Clearly,  then, 

xy  +  yz  +  zx  (5-3) 

transforms  like  A^. 

If,  next,  the  last  term  in  (k-2a.)  is  taken  into  account,  the  transfoma- 
tion  properties  of  the  spin  functions  have  to  be  considered.  The  appropriate 
character  table  is  the  one  for  the  double  Cav  group,  which  is  as  follows: 
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It  can  be  shown  that  the  spin  functions  for  S  =  l/2  transforms  like  E.  Figure 
11  shows  how  the  energy  levels  split  upon  introducing  the  perturbations  in 
(i4-2a)  in  succession  in  the  indicated  order. 


Spherical  Cubic  Trigonal  Spin-Orbit 


/ 

/ 

/ 

/ 

/ 

/ 


Fig.  11.  Split  in  energy  levels  upon  introduction  of  perturbations. 
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VI.  CHARACTER  TABLE 


In  the  previous  section,  we  mentioned  such  terms  as  character  table,  ir¬ 
reducible  representations,  etc.,  without  adequate  explanation.  Here  we  shall 
discuss  these  concepts  in  greater  detail,  with  particular  reference  to  the 
cubic  group.  A  clear-cut  understanding  of  the  character  table  is  necessary 
before  we  can  proceed  to  the  two-electron  Jd  and  three-electron  Jd-^  config¬ 
urations.  We  shall  present  a  discussion  of  what  is  meant  by  the  character 
table,  point  out  some  of  the  important  properties,  and  show  how  this  table 
can  be  used.  Perhaps  it  should  be  emphasized  that  the  character  table  is 
somewhat  like  the  multiplication  table;  if  we  know  how  to  use  the  table,  we 
can  use  it  to  work  out  multiplication  problems.' 

Consider,  then,  the  group  character  table  for  the  cubic  grovp  : 


The  symbols  across  the  top  of  the  table  represent  the  different  classer-  of 
symmetry  operations  and  the  number  in  front  of  the  symbols  gives  the  number 
of  operations  belonging  to  that  class.  A  "class"  of  operations  is  a  set  of 
similar  operations  and  can  often  be  obtained  more  or  less  intuitively.  For 


exanple,  for  the  group  of  operations  that  carry  a  sphere  into  a  sphere,  all 
rotations  of  a  given  angle  independent  of  the  axis  of  rotation  belong  to  the 
same  class.  For  the  cubic  group,  there  are  the  set  of  rotations  about  the 
body  diagonals.  The  two  rotations  of  120°  and  2k0°  about  a  body  diagonal 
will  cariy  a  cube  into  a  cube.  Since  there  is  a  total  of  four  body  diagonals, 
we  might  intuitively  guess  that  all  eight  such  operations  belong  to  the  same 
class.  In  the  above  table,  the  rotations  about  the  body  diagonals  are  in  the 
class  C3.  As  another  example,  consider  Ca,  the  class  of  l80°  rotations  about 
the  coordinate  axes.  Since  there  are  three  axes,  we  should  expect  three  op¬ 
erations  in  class  Ca,  as  Indicated  in  the  table.  Class  Ca  comprises  the  set 
of  180°  rotations  about  the  face  diagonals,  and  C4  the  class  of  90°  rotations 
about  the  coordinate  axes. 

Mathematically,  the  set  of  similar  operations  belonging  to  a  given  class 
is  generated  by  the  similarity  transformation 

=  RgRcRg^  (6-1) 

If  all  symmetry  operations  Rg  belonging  to  the  group  are  used,  a  set  of  dis¬ 
tinct  operations  Rq  will  be  generated.  This  is  the  set  belonging  to  the 
class  C. 

For  the  cubic  group,  there  are  five  classes,  and  the  total  nvunber  of 
symmetry  operations  is 

1  +  8  +  5+  6+  6  =  2k  (6-2) 

The  symbols  Ai,  Ag,  Ei,  Ti,  Ta  are  the  different  Irreducible  representa¬ 
tions  of  the  cubic  group.  In  order  to  see  what  is  meant  by  this,  consider  a 


28 


set  of  polynomials  of  degree  1.  There  are,  as  Is  well  known,  21+1  linearly 
independent  such  polynomials.  Suppose  these  polynomials  are  represented  hy 
•••  h  •••  in  which  n  =  2t+l.  If,  now,  some  operation  Rg  of  the 
cubic  group — and  these  operations  are  simple  linear  substitutions — is  applied 
to  any  one  of  the  functions  )i^i(x,y,z),  the  new  function  will  be  a  polynomial 
of  x,y,z  of  the  same  degree  t;  therefore  the  new  function  can  be  written  as 
a  linear  combination  of  the  original  set  of  functions.  Thus 

Rg(^i)  =  ^  Mij(Rg),:5j  (6-5) 

and  the  operation  Rg,  operating  on  the  set  (j^),  can  be  represented  by  the 
matrix 

M(Rg)  (6-4) 

The  dimension  of  this  matrix  is  clearly  (21+1)  x  (2i+l).  Also,  since  there 
are  24  operations  in  the  cubic  group  there  will  be  24  such  matrices,  one 
corresponding  to  each  operation  Rg.  This  set  of  matrices  is  said  to  be  a 
representation  of  the  cubic  group  because  the  symmetry  operations  of  the  group 
and  the  matrix  representing  the  operation  can  be  set  into  one-to-one  corre¬ 
spondence;  i.e.,  if  Ri,  Rj,  and  Rj^  are  any  three  operations  such  that 

=  %  (6-5a) 

and  if 

Ri  < - >  M(Ri) 

Rj  < - >  M(Rj) 

Rk  i - >  M(Rk)  (6-5b) 
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then 


M(Ri)M(Rj)  =  M(Rk) 


(6-5c) 


This  group  of  matrices  may  have  one  of  two  possible  forms.  It  is  possible 
that  all  the  matrices  in  the  set  will  simultaneously  have  the  form 


(6-6) 


or  can  be  put  into  it.  In  this  form  the  non-zero  blocks  occur  along  the 
diagonal  only.  If  the  matrices  have  or  can  be  put  into  this  form,  the  group 
is  said  to  be  reducible;  if  not,  it  is  said  to  be  irreducible.  For  example, 
for  the  3d  electron  case  discussed  earlier,  1=2,  the  linearly  independent 
functions  are  xy,  yz,  zx,  x^-yS^  and  ^7.^-r^ ,  so  that  24  operations  of  the 
cubic  group  will  generate  24  5-x-5  matrices.  Furthermore,  each  of  the  24 
matrices  will  have  the  form 
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Thus  the  2k  5-X-5  matrices  constitute  a  reduc ibie  representation  of  the  cubic 
group . 

Suppose,  now,  we  consider  the  group  of  2k  matrices  generated  by  Rg  operat¬ 
ing  on  the  set  xy,  yz,  and  zx.  The  matrices  will  be  3-X-5,  as  shown  earlier. 

It  will  be  seen  that  all  the  matrices  cannot  be  partitioned  into  the  form 
(6-6).  The  next  question  is:  By  taking  a  suitable  linearly  independent, 
linear  combination  of  xy,  yz,  and  zx,  is  it  possible  to  generate  a  set  of 
matrices,  all  of  which  will  have  one  of  the  forms  in  (6-8)? 


/  1  .  \ 


(6-8) 


Group  theory  guarantees  us  that  this  is  impossible.  The  set  of  nvunbers  ob¬ 
tained  by  taking  the  diagonal  sum  of  the  representative  ipatrices  will  always 
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be  identical  to  the  set  Tg. 


The  symbols  Ai,  Aa,  Ei,  Ti,  and  Ta  stand  for  a  set  of  matrices  that  are 
in  one-to-one  correspondence  with  the  symmetry  operations  of  the  cubic  group. 
These  five  sets  of  2k  matrices  are  irreducible  in  that  the  matrices  of  each 
set  cannot  all  be  put  into  the  form  (6-6)  simultaneously.  Furthermore,  group 
theory  assures  us  that  whatever  its  dimension,  any  set  of  2h  matrices  which 
represents  the  cubic  group  operations  can  be  decoii5)osed  or  reduced  to  a  linear 
combination  of  the  five  Irreducible  representations.  For  example,  for  poly¬ 
nomials  of  degree  100,  there  will  be  201  linearly  independent  polynomials j 
the  cubic  operations  will  generate  a  set  of  2k  matrices  of  rank  201.  But 
the  rank  of  the  highest  sub-matrix  cannot  be  more  than  three. 

However,  if  we  were  to  select  a  set  of  matrices  from  the  original  2k 
such  that  the  selected  set  satisfies  the  group  property,  then  this  set  could 
possibly  be  reducible.  For  example,  consider  the  rotations  about  a  body 
diagonal.  If  the  axis  of  rotation  is  the  (ill)  direction,  then 
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The  above  syimnetry  operations  constitute  the  group  C3,  whose  character  table 
is  as  follows: 


E 

2C3 

A 

1 

1 

E 

2 

-1 

xy,yz,zx 

3 

0 

=  A+E 

Mathematically,  this  means  that  if  an  appropriate  linear  combination  of  xy, 
yz,  and  zx  is  taken,  the  three  matrices  can  be  reduced.  As  indicated  earlier, 
the  appropriate  linear  combinations  are 

fo  =  xy  +  yz  +  zx 
f].  =  -yz  +  zx 
fa  =  yz  +  zx  -  zxy 

fo 

1  ^  1  ^ 

-zx  +  xy  =  -  2 

zx  +  xy  -  2yz  =  |  fi  -  |  fa 


33 


so  that 


Since 


(6-10) 


(6-11) 


(6-9a) 


all  three  matrices  can  be  partitioned  as  indicated  by  the  dotted  lines. 

The  characters,  or  the  traces,  satisfy  the  vertical  and  horizontal  or¬ 
thogonality  relations.  Let 


x(Ci,  r) 


(6-12) 


represent  the  charsicter,  or  the  trace,  of  the  matrix  which  represents  opera¬ 
tion  Rg  in  class  Ci  for  the  irreducible  representation  f.  Then  the  vertical 
orthogonality  relation  states  that 


X  ’^(Ci,r)x(Cj,r) 
r 


G 

n(Ci) 


6lj 


(6-13) 


3^ 


in  which  G  is  the  order  of  the  group,  or  the  number  of  symmetry  operations  of 
the  group,  and  n(Ci)  represents  the  number  of  symmetry  operations  in  class 
Ci.  For  example,  if  the  numbers  in  the  column  3C2  are  multiplied  by  the  cor¬ 
responding  numbers  in  6C4  and  added,  we  find 

1-1  +  1(-1)  +  2(0)  +  (-1)(1)  +  (-1)(-1)  =  0  (6-14) 

On  the  other  hand,  taking  the  sum  of  the  squares  of  the  characters  in  J)C2,  we 
find 

1^  +  1^  +  2^  +  (-1)^  +  (-1)^  =  8  (6-15) 

which  is  24/3,  This  orthogonality  relation  applied  to  the  class  E  gives 

1^+1^+22+32+52  =  24  (6-16) 

which  states  that  the  sum  of  the  squares  of  the  degrees  of  the  irreducible 
representations  is  equal  to  the  order  of  the  group. 

The  horizontal  orthogonality  relation  states  that 

^  n(Ci)x(Ci,ra)x(Ci,rp)  =  G6ap  (6-17) 

i 

For  exaji5)le,  if  we  take  the  representations  Aq  and  we  find 

1(1)(3)  +8(1)(0)  +3(1)(-1)  +6  (-1)(1)  +6(-l)(-l)  =  0  (6-18) 

On  the  other  hand,  taking  the  square  of  Ta,  we  find 

1.3^  +  8.02  +  3(-1)2  +  6(1)2  +  6( -1)2  =  24  (6-19) 
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This  theorem  is  important  in  determining  the  irreducible  components  of  a  re¬ 


ducible  representation.  Suppose  that  X(CiD)  is  the  character  of  the  matrix 
of  an  operation  in  class  for  some  representation  D  of  the  cubic  group.  If 
the  latter  is  reducible,  its  characters  are  given  by  the  sum  of  the  irreduc¬ 
ible  components.  Thus,  if  a( r)  represents  the  number  of  times  the  Irreducible 
representation  T  is  contained  in  D,  then 


X(Ci 


,D)  -  X 


a(ra)x(Ci,ra) 


(6-20) 


a 


If  we  multiply  this  equation  by  n( C^) X( Cj^ , Tp)  and  sum  over  classes,  we  obtain 

n(Ci)X(Ci,D)X(Ci,rp)  =  Y  n(Ci)X(Ci,rp)a(ra)X(Ci,ra) 
i  a,i 


=  I 

a 


therefore 


=  qY  n(q)X(Ci,D)X(Ci,rp)  (6-21) 

i 

In  many  instances,  the  Irreducible  representations  can  be  determined  by  in¬ 
spection. 

Another  important  theorem  is  concerned  with  the  reduction  of  a  product 
representation.  For  example,  this  problem  will  arise  when  we  consider  the 
two-  and  three-electron  cases.  Since  the  wave  functions  will  be  products  of 
wave  functions  of  the  Individual  electrons,  we  shall  be  concerned  with  the 
transformation  properties  of  such  product  functions.  The  theorem  simply 
states  that 
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x(Ci,rarp)  =  x(Ci,ra)x(Ci,rp) 


(6-22) 


and  the  reduction  of  this  is  carried  out  by  the  recipe  stated  earlier.  This 
theorem  is  useful  in  determining  the  condition  under  which  certain  integrals 
can  be  expected  to  vanish.  In  the  calculation  of  the  energy  levels,  we  need 
the  integrals 

f  ;FiV^|fJdT  (6-23) 

The  perturbation  potential  is  invariant  under  the  symmetry  operation,  and  so 
in  the  case  of  the  cubic  group  it  belongs  to  the  identity  representation  Ai. 
Now  it  can  be  verified  that 

Aifi  =  Ti  (6-24) 

where  rj_  is  Ai,  Ag,  Ei,  Ti,  or  Tg.  Furthermore,  it  can  be  shown  that 
for  i  j  does  not  contain  Ai.  Now  V\}rj  contains  but  the  integral  will 
contain  Ai  only  if  i  =  j;  furthermore,  the  integral  will  be  automatically 
zero  if  Ai  is  not  contained  in  the  product  representation.  This,  then,  tells 
us  that  in  the  perturbation  calculation  we  need  to  calculate  only  those  in¬ 
tegrals  which  connect  states  belonging  to  the  same  representation. 

The  theorem  given  in  (6-22)  is  also  useful  in  determining  the  selection 
rules  for  electric  dipole  transitions.  For  this  the  relevant  integral  is 


/ 


(6-25) 


in  which  i  and  J  are  the  two  states  involved  in  the  transition.  Now  r  trans¬ 
forms  like  Ti  so  that  the  integral  will  be  zero,  unless  the  product  of  \|fj^  and 
\|rj  (or  more  precisely,  Fj^Fj)  contains  Tj..  Using  the  cubic  group  character 
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tatle,  ve  find 

A2T1  =  Ts 

AqE  =  E 

A2T2  =  Ti 

ETi  =  Ti  +  T2 

Er2  =  Ti  +  T2 

aI  =  Ai 

Eg  =  Ai  +  Ag  +  E 

Ti  =  Ai  +  E  +  Ti  +  T2 

Ti  =  Ai  +  E  +  Ti  +  T2 

T1T2  =  A2  +  E  +  Ti  +  T2  (6-26) 

The  allowed  transitions  then  are 

A2  ^ - >T2 

E  < - >Ti 

E< - >  T2 

Ti  < - ^Ti 

T2< - >  T2 

Ti< - >T2  (6-27) 

We  see,  then,  that  the  levels  E  and  Tg  brought  about  by  the  splitting  of  the 
5<3  levels  can  give  rise  to  an  allowed  transition,  despite  the  fact  that  Ai  =  01 
Of  course  this  is  true  for  the  cubic  group  0,  which  has  no  center  of 
symmetry.  For  the  cubic  group  0^,  the  levels  derived  from  the  5d  orbitals  be- 
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long  to  the  even  (gerade)  representations  Eg  and  Tag,  whereas  r  transforms 
according  to  the  odd  (ungerade)  representation  Tiu,  so  that  Integrals  of  the 
form 


dr 


vanish  identically  and  electric  dipole  transitions  are  not  allowed  between 
states  of  these  levels. 
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VII.  THE  CONFIGURATION 


Earlier  we  indicated  that  the  3<i-electron  energy  level  is  split  into 
levels  belonging  to  the  cubic  irreducible  representations  E  and  Ts*  In  our 
notation  henceforth  we  shall  speak  of  the  tg  or  the  e  electron  depending  upon 
whether  the  electron  occupies  one  of  the  T2  or  E  orbitals.  Since  there  are 
three  Ta  and  two  E  orbitals,  it  is  clear  that  there  can  be  all  together 
2x3=6  electrons  occupying  the  T2  orbitals,  and  2x2=4  electrons  oc¬ 
cupying  the  E  orbitals.  Consequently,  we  may  speak  of  the  ta  and  the  e  sub¬ 
shells  of  the  d-electron  shell.  As  in  atomic  spectroscopy,  we  shall  con¬ 
sider  next  the  possible  electron  configuration  arising  from  two  equivalent 
ta  electrons. 

The  basic  Schroedinger  equation 

Hf  =  Eijr  (4-1) 

we  wish  to  solve  is  such  that 


H  =  Hi  +  Ha  +  - 

1-12 


(7-1) 


in  which 


Hi  =  -  ^  +  V(ri)  +  Fc(x|+y|+z|-  |  rf) 


(7-2) 


The  method  of  attack  is  to  expand  the  solution  in  terms  of  product  functions 
of  the  two  electrons.  For  the  tf  configuration,  there  are  nine  product  func¬ 
tions  obtained  by  multiplying  the  set  (xy) i,(yz) i,(zx) 1  into  (xy)2,(yz)2, 
(zx)2.  The  indices  1  and  2  refer  to  the  two  ta  electrons.  Earlier,  the 


symmetry  properties  of  the  Hamiltonian  were  exploited  to  obtain  information 
about  the  nature  of  energy  levels  and  wave  functions.  The  Hamiltonian  (7-1) 
is  invariant  under  the  operations  of  the  cubic  group,  so  that  some  of  the 
previous  arguments  can  be  used. 

However,  we  note  that  the  Hamiltonian  exhibits  an  additional  symmetry 
property  of  being  invariant  when  the  indices  1  and  2  are  permuted.  This 
means  that  the  solutions  of  (4-1)  must  constitute  the  basis  functions  for 
the  symmetric  group  Ss.  For  this  group,  there  are  only  two  operations:  the 
identity  operation,  represented  by  the  permutation  (l)  (2),  and  the  operation 
exchanging  the  indices  1  and  2,  represented  by  (12).  Clearly,  then,  there 
are  only  two  irreducible  representations,  and  the  character  table  is  as  fol¬ 
lows  : 


This  table  shows  that  the  irreducible  representations  are  one-dimensional, 
l.e.,  the  basis  functions  are  such  that  each  transforms  into  itself  or  into 
its  negative  upon  permutation  of  the  indices  1  and  2.  These  are  the  familiar 
symmetric  and  antisymmetric  representations.  Consider  now  the  effect  of 
these  operations  upon  the  nine  product  functions : 
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(yz)i(yz)2  (zx)i(yz)2  (xy)i(yz)2 

(yz)i(zx)2  (zx)i(zx)2  (xy)i(zx)2  (7-3) 

(yz)i(xy)2  (zx)i(xy)2  (xy)i(xy)2 

The  two  permutations  are  represented  by  9-x-9  matrices.  The  character,  or  the 
trace,  is  simply  equal  to  the  number  of  functions  that  are  not  affected  by 
the  permutations.  Consequently  for  the  class  (1^)  the  character  is  S,  since 
none  of  the  functions  change  when  left  alone.  For  the  permutation  (12)  only 
three  functions — (yz)i(yz)2,  (zx)i(zx)2,  and  (xy)i(xy)2 — transform  into  them¬ 
selves.  Consequently  the  character  for  class  (2)  is  By  inspection  (there 
is  no  need  here  to  appeal  to  the  general  theorem  enunciated  earlier.')  it  is 
seen  that  the  reducible  representation  contains  six  symmetric  and  three  anti¬ 
symmetric  irreducible  representations.  Since  these  representations  are  one¬ 
dimensional,  the  nine  functions  of  (7-5)  can  be  grouped  into  six  symmetric  and 
three  antisymmetric  functions.  The  symmetric  functions  are 

(yz)i(yz)2  (yz)i(zx)2  +  (yz)2(zx)i 

(zx)i(zx)2  (yz)i(xy)2  +  (yz)2(xy)i  (7-4) 

(xy)i(xy)2  (zx)i(xy)2  +  (zx)2(xy)i 

and  the  antisymmetric  functions  are 

(yz)  i(zx)2  -  (yz)2(zx)  1 
(zx)i(xy)2  -  (zx)2(xy)i 
(xy)i(yz)2  -  (xy)2(yz)i 
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(7-5) 


Next  consider  transformations  induced  by  the  operations  of  the  cubic  group. 


For  (7"^)  we  shall  obtain  2k  6-x-6  matrices  and  for  (7-5)  24  3-x-3  matrices. 
Our  task  is  to  decompose  these  into  a  sum  of  irreducible  representations. 

The  calculation  of  the  characters  would  be  rather  tedious  if  we  had  to  carry 
out  the  linear  substitutions  mentioned  earlier.  Fortunately,  however,  by 
using  the  result 


and 


X(R,S)  =  i  - 

[X(R) +  X(r2) . 

^  s. 

X(R,A)  =  i  ■ 

[X(R)  f  -  X(R^)  ■ 

(7-6a) 


(7-6b) 


the  necessary  decon^iosition  can  be  readily  carried  through  (see,  for  example, 
Ref.  12,  p.  134).  From  the  cubic  character  table  we  have 


To  the  characters  of  Tf,  we  need  to  add  or  subtract  those  of  X(R^) .  In  so 
doing,  we  note  that 


=  E 

ci  =  C3 
cl  =  E 
=  E 

cf  =  ci 


^3 


giving  the  set  of  numbers  indicated  in  the  row  X(R^).  Consequently,  char¬ 
acters  for  the  symmetric  and  antisymmetric  representations  are 


E 

8C.q 

3Ca 

6C2 

6C4 

s 

6 

0 

2 

2 

0 

A 

3 

0 

-1 

-1 

1 

so  that 

S  =  Ai  +  E  +  Ta  (7-7) 

and 

A  =  Ti  (7-8) 

As  a  check,  we  note  that 

(Ta)^  =  Ai+E+Ti+Ta  =  S+A  (7-9) 

According  to  Paul's  principle  the  overall  wave  function,  including  both 
orbit  Eind  spin,  must  be  antisymmetric  in  the  indices  1  and  2.  Therefore  if 
the  orbital  function  is  symmetric,  the  spin  functions  must  be  antisymmetric, 
and  vice  versa.  The  symmetric  and  antisymmetric  spin  functions  have  spins  1 
and  0,  respectively.  Consequently,  we  obtain  the  allowed  states 

^Ai  +  ^E  +  ^Ta  (7-lOa) 

and 

(7-lla) 

in  agreement  with  the  result  obtained  by  other  procedures. 

For  future  use,  we  shall  restate  our  results  in  terms  of  Young  tableaux 
(see  Refs.  12  and  IJ) .  Using  this  notation  we  have  shown  that 


1*4 


(7-lOb) 


=  ^Ai  +  +  ^Ta 

-  =  ^Ti  (7-llb) 

We  shall  need  these  results  to  determine  the  allowed  states  for  the  Jd^  con¬ 
figurations. 
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VIII.  THE  3<a^  CONFIGURATION 


The  aim  of  this  section  will  he  to  present  a  number  of  arguments  to  make 
the  Orgel-Sugano-Tanabe  diagram  of  Cr^'*’  understandable.  As  has  been  pointed 
out  already,  the  cubic  crystalline  electric  field  will  lead  to  the  splitting 
of  free  ion  energy  levels.  The  Orgel-Sugano-Tanabe  diagrams  show  how  the 
energy-level  splittings  are  affected  by  the  strength  of  the  crystalline  field, 
as  shown  in  Fig.  12. 

Cr  IV:  3d^rF)4S,^F~ll3  B  =  9I8 


Dq/B 


Fig.  12.  Splitting  of  states  of  the  ;3d  '  configuration 
by  an  octahedral  field  (Ref.  5>  Fig.  2). 


To  develop  the  qualitative  ideas,  on  the  basis  of  symmetry  arguments  we  shall 
first  consider  the  case  of  weak  crystalline  fiel'i  and  then  the  case  of  the 


strong  crystalline  electric  field. 


To  carry  out  this  program  we  shall  make  frequent  reference  to  the  group 
character  table  and  to  the  reduction  of  reducibility  representations.  For 
the  sake  of  convenience  the  group  character  table,  the  characters  of  the  re¬ 
ducibility  representation,  and  their  Irreducible  components  are  presented  in 
the  following  cubic  character  table: 


E 

8C3 

3C2 

6C2 

6C4 

A-l 

1 

1 

1 

1 

1 

As 

1 

1 

1 

-1 

-1 

E 

2 

-1 

2 

0 

0 

Ti 

0 

-1 

-1 

1 

Ta 

3 

0 

-1 

1 

-1 

P 

3 

0 

-1 

-1 

1 

Ti 

D 

5 

-1 

1 

1 

-1 

E+Ts 

F 

7 

1 

-1 

-1 

-1 

A2+T1+T2 

G 

9 

0 

1 

1 

1 

Aj+E+Ti+Ts 

H 

11 

-1 

-1 

-1 

1 

E+2T1+T2 

Ag  X  Ag 

1 

1 

1 

1 

1 

As  X  E 

2 

-1 

2 

0 

0 

E 

As  X  Ti 

3 

0 

-1 

1 

-1 

Ta 

As  X  Ts 

3 

0 

-1 

-1 

1 

Ti 

E  X  E 

4 

1 

4 

0 

0 

A1+A2+E 

E  X  Ti 

6 

0 

-2 

0 

0 

T1+T2 

E  X  Ts 

6 

0 

-2 

0 

0 

Ti+Ta 

Ti  X  Ti 

9 

0 

1 

1 

1 

A1+E+T1+T2 

Ti  X  Ts 

9 

0 

1 

-1 

-1 

Aa+E+Ti+Ts 

Ts  X  Ts 

9 

0 

1 

1 

1 

Ai+E+Ti-KCa 

h7 


The  terms  for  the  3<i^  electron  configuration  (see,  for  example,  Ref.  12, 

p.  423)  are  ^P,  (^D)®,  ^F,  ®G,  ^H,  ^D,  and  "^F.  According  to  the  Atomic  Energy 
6  l4 

Tables,  ’  these  levels  are  arranged  as  follows: 


^F  -  - 

- 36490 

-  - 

- 21078 

a^D  -  - 

- 20218 

-  - 

- 15064 

-  - 

- 14185 

■‘P  -  - 

- 14072 

^F . 0 


(The  numerical  indicates  the  term  value  for  the  lowest  J  value  of  the  multi- 
plet.) 

When  this  3“electron  system  is  placed  in  an  electric  field  having  cubic 
symmetry,  the  above  atomic  energy  levels  will  split,  as  Indicated  by  the 
character  table.  For  example,  since 

F  =  As  +  Ta  +  Ti 

the  ^F  groiind  state  will  split  into  three  levels,  as  shown  in  Fig.  13 . 

The  splitting  will  be  similar  for  other  levels.  The  complete  splitting 
scheme  of  the  terms  arising  from  the  3d^  configuration  is  given  in  Fig. 

14. 
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I 

I 


Fig.  15.  Splitting  of  ground  states. 


Fig.  l4.  Complete  splitting  scheme  of  the  terms  arising 
from  the  configuration. 
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Group  theory  does  not  give  the  sequence  of  terms  shown  in  Fig.  l4.  However, 
provided  the  constant  124  is  positive,  detailed  calculations  of  the  sequence 
can  be  obtained  from  pages  I3  and  I7  of  Ref.  I5. 

The  scheme  gives  an  explanation  of  the  left-hand  side  of  the  OST  diagram. 
It  should  perhaps  be  noted  that  Sugano  and  Tanabe^  have  indicated  the  split¬ 
tings  only  of  states  "^F,  ^G,  and  ^F. 

Consider  next  the  right-hand  side  of  Fig.  lU.  To  understand  how  these 
levels  came  about,  let  us  first  recall  an  earlier  remark  that  the  effect  of 
the  cubic  crystalline  field  is  to  remove  the  equivalency  of  some  of  the  d- 
electrons.  In  a  spherically  symmetric  electric  field  and  d-shell  can  accom¬ 
modate  as  many  as  ten  equivalent  electrons.  On  the  other  hand,  in  a  cubic 
field  six  electrons  will  be  affected  differently  from  the  remaining  four,  so 
that  often  the  terms  ta-subshell  and  e-subshell  are  introduced.  The  elec¬ 
trons  in  the  ta-subshell  will  be  equivalent  to  one  another,  but  not  to  those 
in  the  e-subshell.  It  is  easy  to  see,  then,  that  the  possible  electron  con¬ 
figurations  are 

(ta)^,  (ta)^e,  ta(e)^,  (e)^ 

The  right  superscript,  as  in  atomic  spectroscopy,  represents  the  number  of 
electrons  with  the  indicated  orbital.  For  example,  ta(e)^  means  that  there 
is  one  electron  occupying  the  ta  level,  and  two  electrons  occupying  the  e 
levels.  We  shall  discuss  the  allowed  states  associated  with  each  configura¬ 
tion. 

Consider  first  the  configuration  tf.  The  Hamiltonian  will  be  invariant 
under  the  pemiutation  of  the  indices  1,  2,  and  3>  so  that  the  wave  functions 
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must  constitute  the  basis  for  the  irreducible  representations  of  the  symmetric 


group  S3.  This  group  consists  of  six  permutations,  belonging  to  three  classes 
as  indicated  below: 


Classes 

(1^) 

(2,1) 

(5) 

Permutation 

E 

(12) 

(123) 

(13) 

(132) 

(23) 

There  are,  accordingly,  three  irreducible  representations  of  degrees  1,  2, 
and  1,  respectively.  These  representations  will  be  denoted  by  the  symbols 
[3],  [2,1],  and  [1^]  respectively.  The  group  character  table  (see,  for 
example.  Ref.  12)  is 


(1®) 

3(2,1) 

2(3) 

[3] 

1 

1 

1 

[2,1] 

2 

0 

-1 

[1^] 

1 

-1 

1 

ti 

27 

9 

3 

The  characters  of  the  reducible  representation  tf  given  in  this  table 
are  obtained  by  noting  how  the  three  electron  wave  functions  transform  under 
permutations  of  the  group  S3.  As  mentioned  earlier,  the  tg  wave  functions 
are  given  by 

I  =  (yz) 

II  =  ( zx) 

?  =  (xy) 
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The  3-electron  wave  functions  are  given  by  the  product  of  1-electron  wave 
function.  Since  the  wave  function  for  each  electron  can  be  any  one  of  the 

3 

three  functions  indicated  above,  these  are  3  =27  product  functions,  so  that 

the  character  for  the  identity  operation,  E,  is  27.  Consider  next  the  char¬ 
acter  for  the  permutation  of  class  (2,1) ,  and  in  particular  consider  the 
permutation  ( 12) .  The  character  for  this  operation  and  consequently  for  the 
class  (2,1)  is  given  by  the  number  of  functions  that  are  not  changed  by  the 
permutation.  This  will  occur  if  the  orbital  part  for  electrons  1  and  2  is 
the  same.  A  typical  function  is,  say, 

1(1)  X  1(2)  X  ti(3) 

The  orbital  of  electrons  1  and  2  can  be  one  of  the  three  functions,  and 
electron  3  can  also,  independently,  have  any  one  of  the  three  orbitals.  Con¬ 
sequently,  the  number  of  product  functions  remaining  unchanged  by  the  per¬ 
mutation  (12)  is  9>  which  means  that  the  character  for  this  class  is  9*  Fin¬ 
ally,  the  character  for  class  (3)  is  the  number  of  functions  that  are  not  af¬ 
fected  by  the  permutation  (123)  or  (132).  Clearly,  there  are  only  three  such 
functions,  given  by 

1(1)  X  ^(2)  X  1(3),  Tjd)  X  ti(2)  X  ti(3),  ^(D  x  ^(2)  x  ^(3) 

and  hence  the  character  is  3.  Using  the  theorem  enunciated  earlier,  we  find 
that 

(ta)^  =  10  [3]  +  8  [2,1]  +  [1^] 
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or,  in  terms  of  Young's  tableaux. 


(ta)^  = 

(10)  (8)  (1) 

The  dimensionality  of  the  representation  given  by  each  shape  is  given  by  the 
numbers  underneath.  The  representation  [3]  or  ,  for  which  the 

orbital  part  is  totally  symmetric,  is  forbidden  by  Pauli's  principle.  Or  to 
use  the  language  of  tableau  calculus,  in  order  for  the  total  wave  function 
to  satisfy  Pauli's  principle,  the  tableau  representing  the  spin  part  must  be 
the  transpose  of  the  orbital  part.  Furthermore,  since  the  electron  has  only 
two  possible  spin  orientations,  the  electron  spin  tableau  must  not  have  more 


as  the  only  allowed  orbital  tableaux.  Now 


belongs  to  of  the  cubic  irreducible  representation. 
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This  can  be  seen  as  follows.  The  tableau  represents  the  totally  antisym¬ 
metric  space  function  and  is  given  by 


ta 


Sp  P[(xy)i(yz)2(zx)3] 


(xy)  i(xy)2(xy)3 


N 


(yz)  i(yz)  2(72)3 


(zx)i(zx)2(zx)3 


where 


+1  for  an  even  permutation 
-1  for  an  odd  permutation 
Selecting  now  one  operation  out  of  each  class  of  the  cubic  group  and  operating 
on  the  we  get ; 


Operation 


Result  Character 

Eta  =  ta  1 


C3: 


C3ta  =  ta 


1 


C2 : 


-X 

■>  -y 

z 


C2ta  =  ta 


1 


/X  — =>  y  \ 

C2:  f  y - >  X  j  Cata  =  "'I'a  -1 


/X - >y\ 

C4:  (  y - C4ta  =  -ta  -1 

These  are  the  characters  of  the  irreducible  representation  A2. 


5^ 


The  irreduclhle  cubic  representations  contained  in  the  other  tableaux 


can  be  readily  obtained  as  follows.  The  outer  product  gives 


© 


Now 

belongs  to  T2 


belongs  to  ^T^ 


and  furthermore 

Ti  +  Ta  =  As  +  E  +  Ti  +  Ta 


so  that 

=  +  ^Ti  +  ^Ta 

Collecting  our  results,  we  have  shown  that  the  allowed  configurations  of  tf 
are 

ti  =  4^2  +  2e  +  2Ti  + 

This  is  in  agreement  with  Bethe's  result,  which  is  obtained  by  a  very  differ¬ 
ent  procedure. 
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that 


2 

The  states  for  the  configuration  (ts)  e  can  he  obtained  easily  by  noting 

(t2)^  =  "^Ti  +  +  ^Ei  + 

and  by  making  use  of  the  results  for  product  representations.  Thus,  the  ad¬ 
dition  of  an  e-electron  to  the  state  ®Ti  leads  to 

(^Ti)  (E)  =  ^Ti  +  ^Ts  +  “^Ti  +  ‘^T2 

The  others  are 

(^Ai)  (E)  =  ^E 

(^Ei)  (E)  =  ^Ai  +  ^As  +  ^E 

(^Ta)  (E)  =  ^Ti  +  ^Ta 

Carrying  through  a  similar  analysis  for  t2(e)  we  find  that 

t2(e)^  =  ^Ti  +  ^Ti  +  ^Ta  +  ^Ti  +  ^Ts 

And  finally  for  e^,  the  single  state 

e^  =  ^E 

In  summary  we  have  found  that  the  states  arising  for  the  weak-field  case 
are  as  follows : 
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I 

I 


From  the  strong-field  analysis  we  obtain; 


Total  12111455 


The  two  totals  agree,  as  they  should.  Each  of  the  states  ^Ta,  and 

^Ae  occur  only  once.  These  are  the  states  Indicated  by  the  straight  lines 
connecting  the  left-  and  right-hand  sides  of  the  OST  diagram.  Since  there 
are  two  "^Ti  states,  the  associated  energy  levels  are  given  by  the  solution  of 
a  2-X-2  determinant.  Furthermore,  according  to  well  known  results  in  quantum 
mechanics,  the  energy  states  of  different  symmetry  do  not  interact,  but  those 
of  the  same  symmetry  tend  to  repel.  Consequently,  on  the  left;-hand  side, 
the  two  '*Ti  states  should  start  out  parallel  to  each  other,  but  the  separa- 
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tion  should  increase  with  increasing  strength  of  the  crystalline  electric 
field. 

The  energy  levels  for  the  ^E,  ^Ti,  and  ^Ts  states  are  obtained  by  solv' 
ing  4-X-4  (for  ®E)  and  5-x-5  (for  ^Ti  and  ^Ta)  determinants. 

We  shall  now  analyze  the  energy  levels  in  greater  detail.  The  term 
values  for  the  d^  configuration  are  given  in  Ref.  11,  pp.  206  and  233  (see 
also  Ref.  l6) .  These  values  are: 

-15B 

0 

-IIB  +  3C 

2p,2H  -  6b  +  3C 

a^D  5B  +  5C  --^fl93B2TTBC~n;^ 

9B  +  5C 

b^D  5B  +  5C  +n/l95B^  +  8bc  +  4c^ 


Comparison  of  these  results  with  the  experimental  spectrum  given  earlier 
shows  very  clearly  the  inadequacy  of  even  the  atomic  theory.  If  the  term 
values  for  the  ^F,  '^P,  and  ^G  states  Eire  used,  the  following  numerical 
values  are  obtained : 

B  =  917.56  cm"^ 

C  =  3678,67  cm'^ 


and 


7 


C 

B 


4.01 


On  the  other  hand,  Sugano  and  Tanabe 
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have  used 
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B  =  918  cm' 


and 

7  =  5  =  ^.50 

for  their  analysis. 

For  the  octahedral  coordination,  the  crystalline  field  energy  of  the  ta 
and  e  orbitals  are  -iiDq  and  6Dq  respectively.  These  values  are  chosen  hy 
taking  the  separation  of  ta  and  e  to  be  lODq,  with  the  center  of  gravity  at 
zero  (or  in  matrix  language,  the  trace  of  the  cubic  field  energy  matrix  is 
taken  to  be  zero.)  Thus 


therefore 


3E(t2)  +  2E(e)  =  0 

E(e)  -  E(t2  =  lODq 

E(e)  =  6Dq,  E(t2)  =  -UDq 


Then  the  crystalline 


are  given  by 


E(t 


;i  e^"”)  =  (-i+Dq)n  +  (3-n)(6Dq) 


so  that  we  obtain: 


n 

3 

2 

1 

0 


Configuration 

(t2)^ 


(ta)^ 

ta(e)^ 


Crystalline  Energy 
-12Dq 

-  2Dq 

+  8Dq 

l8Dq 
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This  means  that  for  very  large  values  of  Dq  (i.e.,  for  a  large  crystalline 
field)  the  terms  arising  from  the  above  configurations  will  cluster  about 
the  values  indicated  in  the  last  coltimn. 
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PART  II 


X-RAY  PRODUCED  VO'^^  IN  SINGLE  CRYSTAIS  OF  ZINC 

AMMONIUM  SULFATE— ZnS04(NH4)  2304- 6H20 


i3y 


Robert  Borcherts 


Foreword 


The  purpose  of  this  report  is  to  present  the  results,  to  date,  of  a  set 
of  electron  spin  resonance  experiments  on: 

(a)  VO'*'^  and  in  single  crystals  of  zinc  ajnraonium  sulfate,  and 

(b)  X-ray  irradiated  VO  in  single  crystals  of  zinc  ammonium  sulfate, 
magnesium  ammonium  sulfate,  and  zinc  potassium  sulfate. 

Conclusions 

+2 

Upon  X-ray  irradiation  of  the  single  crystals  containing  VO  it  was 
found  that  the  X-ray,  or  the  subsequent  energetic  electron  that  it  produces, 
breaks  the  VO  bond,  removing  the  oxygen  from  the  VO'*’^  site  and  leaving  V'*’^ 
behind . 

Since  the  'V*'^  produced  by  Irradiation  and  the  crystalline  field  surround¬ 
ing  it  is  identical  ( as  far  as  EPR  measurements  can  determine)  to  that  found 

+2 

in  single  crystals  of  "grown"  V  in  zinc  ammonium  sulfate,  the  oxygen  must 
be  at  a  distance  so  removed  that  the  crystalline  field  is  unaffected  by  its 
presence. 

Figure  1  shows  the  EPR  spectrum  along  the  axis  for  the  VO  crystal, 

I  p 

the  irradiated  VO  crystal,  and  the  V  crystal,  depicting  the  results 
quite  vividly. 
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EPR  SPECTRA 
AL0N6  K.  AXES 


I.  BACKGROUND 


The  first  successful  application  of  microwave  techniques  to  the  study  of 
the  centimeter  and  millimeter  range  of  energy  level  separation  was  by 
Zavolsky^  in  19^5-  Since  then  the  study  of  electron  spin  resonance  spectra 
of  various  ions  in  crystalline  and  non-crystalline  solids  has  made  a  major 
contribution  to  knowledge  in  the  field  of  solid  state  physics. 

One  ion  of  particular  Interest  is  that  of  vanadium.  The  +2,  +3,  and  +4 
spin  states  of  (abundance  99*76^^  I  =  7/2)  have  electronic  spins  of  3/2, 

1,  and  1/2  respectively.  Consequently  the  electron  spin  resonance  signature 
of  each  oxidation  state  is  readily  recognized — in  non-cubic  crystalline  fields 
by  the  number  and  angular  dependence  of  the  fine  structure  groups  as  well  as 
by  the  separation  of  the  hyperfine  structure,  and  in  cubic  crystalline  fields 
by  the  separation  of  the  hyperfine  structure  (see  Figs.  2  and  3)* 

P 

In  1950 >  Bleaney^  Ingreun,  and  Scovll  reported  on  the  EPR  spectrum  of 

V^2 

in  ZnS04(NH4)2S04*6H20,  the  same  host  crystal  used  in  this  set  of  experi¬ 
ments.  In  this  crystal  the  six  water  molecules  form  a  distorted  octahedron 
surrounding  the  V^2  ion,  producing  a  rhombic  crystalline  field.  As  a  result 
the  2S+1  degenerate  ground  state  splits  into  two  spin  doublets,  Ms  =  ±3/2  and 
Ms  =  ±1/2  respectively,  separated  by  2D.  Applying  a  magnetic  field  splits 
the  remaining  degeneracy  and  application  of  microwave  energy  at  a  frequency 
V  gives  rise  to  three  fine- structure  groups  at  g  =  2,  corresponding  to  the 
selection  rule  AM  =  ±1.  Because  of  the  7/2  value  of  the  nuclear  spin  of 
each  of  the  fine-structure  groups  are  composed  of  eight  lines  (21+1)  with  a 
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Fig.  2.  and  V+2  energy  levels. 
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separation  of  approximately  90  gauss. 

In  1958,  Zverev  and  Prokhorov^  and  later  Lambe,  Ager,  and  Kikuchi^  reported 
on  in  sapphire  (a-AlsOs)  .  And  In  i960,  Lambe  and  Kikuchi^  found  a  small 
amount  of  in  AI2O3.  They  were  also  able  to  show  that  by  X-ray  irradiation 
or  gamma  irradiation  some  is  produced.  Also  in  i960,  Gerrltsen  and  Lewis^ 
and  Zverev  and  Prokhorov^  reported  on  in  rutile  (TiOa). 

In  sapphire  six  oxygens  surround  each  aluminum  site  in  a  distorted  octa¬ 
hedron,  giving  rise  to  a  trigonal  crystalline  electric  field.  If  sub¬ 

stituted  for  the  Al'''^  ion,  the  trigonal  component  of  the  crystalline  field 
splits  the  T2  energy  level  into  a  doublet  E  and  a  singlet  B.  This  splitting, 
acting  through  the  L-S  coupling,  causes  the  ground  state  to  split  into  a  lower 
spin  singlet  Ms  =  0  and  an  upper  spin  doublet  Mg  =  ±1,  separated  by  ~10/cm. 

An  applied  magnetic  field  causes  the  Mg  =  ±1  level  to  split  and  as  a  result  a 
forbidden  electron  spin  resonance  transition  corresponding  to  AM  =  ±2  can  be 
observed  at  a  magnetic  field  of  H  =  hv/2gp.  Since  the  ground  state  is  sepa¬ 
rated  by  ~10/cm  the  relaxation  time  of  this  transition  is  so  short  that  in 
order  to  observe  the  spin  resonance  transitions  the  experiments  must  be  per¬ 
formed  at  very  low  temperatures — approximately  h.2°K.  If  the  temperature  is 
too  low,  less  than  2°K,  the  Ms  =  ±1  state  becomes  so  depopulated  that  no  sig¬ 
nal  is  observed  at  all,' 

In  rutile,  six  oxygen  atoms  surround  the  Ti'*'^  site  in  a  distorted  octa¬ 
hedron,  producing  an  orthorhombic  crystalline  electric  field.  Since  has 
only  one  3<3  electron  and  thus  a  spin  of  l/2,  when  it  is  substituted  for  the 
Ti'*'^  ion  only  one  fine  structure  group  corresponding  to  a  AM  =  ±1  transition 
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Is  observed.  However,  if  tetragonal  symmetry  is  present  the  relaxation  time 
may  be  long  or  short  depending  on  whether  or  not  the  tetragonal  symmetry  is 
"squashed"  or  "elongated"  (see  Fig.  3).  For  rutile  the  tetragonal  field  ap¬ 
pears  to  be  elongated  since  the  experiment  has  to  be  performed  at  liquid  nitro¬ 
gen  temperatures  (77°K)  in  order  to  observe  the  EPR  spectrum.  The  observed 
hyperfine  spectra  is  very  anisotropic,  having  a  separation  of  157  gauss  along 
the  z-axis  and  35  gauss  and  48  gauss  along  the  x-  and  y-axes  respectively. 

The  g-values  (gx  =  1*915;  gy  =  1*913;  gz  =  1*956)  reflect  the  rhombic  sym¬ 
metry  of  the  crystalline  field. 

Although  the  "valence  state"  of  vanadium  in  the  vanadyl  radical  VO  is 
also  +4,  the  strong  cylindrical  crystalline  field  formed  by  the  vanadium  and 
the  oxygen  produces  an  energy  level  structure  quite  different  from  that  of 
in  a  cubic  field®  (see  Fig.  3)  • 

The  first  electron  spin  resonance  spectrum  of  the  vanadyl  radical  (VO"^^) 

was  reported  by  Garifianov  and  Kozyrev^  in  a  frozen  aqueous  solution  and  later 

by  Fake  and  Sands^*^  in  aqueous,  acetone,  and  ether  solutions;  by  O'Reilly^^ 

in  vanadyl  etioporphyrin  dissolved  in  benzene  and  high  viscosity  oil;  by  Faber 

and  Rogers^^  in  various  adsorbers  (charcoal,  Dowex-59;  IR-4B  and  IR-lOO) ;  and 

by  Roberts,  Koskl,  and  Caughey^5  in  some  vanadyl  porphyrins.  As  mentioned 
+P 

above,  the  VO  radical  produces  a  very  strong  cylindrical  of  axial  field  so 
that  it  matters  little  what  other  kind  of  field  surrounds  the  VO"^^.  This  is 
evidenced  by  the  similar  spectra  found  by  these  investigators. 9"13  The  re¬ 
sults  of  these  investigations  are  that  the  gj_  and  g,^  values  are  1. 98-1. 99 
and  1.88-1.93  respectively,  and  that  there  is  a  hyperfine  splitting  of  l60- 
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200  gauss  parallel  to  the  z-axis  and  of  60-85  gauss  perpendicular  to  the  z- 
+2 

or  VO  -axis.  The  variation  in  these  values  is  suggested  to  be  due  to  the 

extent  of  covalent  bonding  between  the  vanadyl  ion  and  the  surrounding  lig- 

andSo  Because  the  samples  used  in  all  these  investigations  were  either 

powders  or  solutions,  the  EPR  spectrum,  is  either  an  average  of  the  randomly 

oriented  VO  radicals  in  the  frozen  or  powdered  samples  (see  Fig.  4a),  or  a 

+? 

time-averaged  spectrum  as  a  result  of  the  motion  of  the  VO  in  solutions. 

Thus,  because  of  the  need  for  EPR  information  on  oriented  VO"^^  and  be¬ 
cause  of  the  success  of  Lambe  and  Kikuchi^  and  Wertz,  Auzins,  Griffiths,  and 
Orten^^  in  producing  the  various  oxldaxion  states  of  vanadium,  the  experiments 
reported  here  were  undertaken.  Underlying  these  immediate  reasons  is  perhaps 
a  deeper  one — that  of  a  systematic  study  of  the  solid  state  chemistry  of 
vanad ium. 
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II.  EXPERIMENTAL  METHODS 


A.  CRYSTAL  GROWTH 

Single  crystals  of  0.05,  0.1,  0.5,  and  1.0^  (vanadium  to  zinc)  concentra¬ 
tions  of  ZnS04(NH4) 2S04-6H20  were  grown  from  a  water  solution  into  which  re¬ 
agent  grade  ZnS04'7H20,  (NH4)2‘S04  and  V0S04-2H20  were  added  in  the  prescribed 
amounts.  Light  blue  crystals  up  to  3-^  in  length  were  usually  obtainable 
in  1-2  days.  Zinc  potassium  sulfate  and  magnesium  ammonium  sulfate  crystals 
containing  1.0^  VO  ^  were  also  grown  by  this  method. 

In  contrast  to  the  ease  of  growing  the  vanadyl  crystals,  the  growth  of 

the  1.0^  and  2,0'^  vanadous  ammonium  sulfate  crystals  required  a  greater  degree 

+P 

of  skill  due  to  the  rapidity  at  which  V  becomes  oxidized.  Cathodic  reduc¬ 
tion  of  the  vanadyl  solution,  in  which  the  cathode  and  the  anode  are  separated 
by  a  porcelain  cup,  followed  by  evaporation  at  5°C  in  a  carbon  dioxide  atmos¬ 
phere  was  the  method  used  to  grow  these  crystals. 

B,  CRYSTAL  STRUCTURE  AND  ORIENTATION 

Figure  5  shows  the  two  molecules  contained  in  the  unit  cell  of 
MgS04(NH4) 2S04'6H20.  The  position  of  the  atoms  is  taken  from  the  data  in 
Wykcoff.^5  The  octahedron  of  water  molecules  surrounding  the  divalent  metal 
ions  forms  the  crystalline  electric  field  which  determines  the  constants  in 
the  spin  Hamiltonian,  That  is,  the  data  obtained  from  paramagnetic  ions 
placed  in  these  positions  directly  reflects  the  environment  and  the  interac¬ 
tion  of  the  electrons  with  the  environment.  Measurements  of  EPR  experiments 
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Fig.  5«  Unit  cell  of  J^(NH4) 2(804)2*61^0. 


can  be  made  with  great  precision  so  that  very  minute  changes  in  the  environ¬ 


ment  can  be  observed. 

As  shown  in  Fig.  6  these  double  sulfates  or  Tutton  salts  grow  with  well 
recognized  faces-*-®  so  that  orientation  of  the  crystal  in  the  cavity  becomes 
a  task  accomplished  with  relative  ease.  Figure  6  also  shows  the  device  used 
for  positioning  the  22-inch  quartz  rod  to  the  crystal  to  within  a  degree  or 
two  of  the  desired  orientation.  After  placing  the  crystal  in  the  cavity  the 
final  adjustments  of  the  crystal  orientation  are  made  by  moving  the  quartz 
rod  until  the  desired  spin  resonance  spectrum  is  observed  on  the  recorder,  The 
experimental  data  that  can  be  recorded  are;  (a)  the  magnetic  field  at  reso¬ 
nance,  (b)  the  angle  of  the  magnetic  field  to  some  axis,  (c)  the  frequency  of 
the  klystron,  (d)  the  first  derivative  of  the  line  shape  as  presented  on  the 
recorder,  and  (e)  the  angle  that  the  c-axis  makes  with  the  horizontal.  The 
angle  that  the  c-axis  of  the  crystal  makes  with  the  horizontal  is  observed 
through  the  optical  hole  in  the  side  of  the  cavity  by  means  of  a  surveyors 
transit. 

C.  EPR  EXPERIMENTAL  SETUP 

Electron  spin  resonance  spectra  and  magnetic  field  measurements  were  made 
at  room  temperature  with  a  Varian  V-4500  EPR  Spectrometer  and  a  Varian  Model 
F-8  Fluxmeter  connected  to  a  Berkeley  7580  Transfer  Oscillator  and  a  Berkeley 
7370  Universal  E^ut  and  Timer.  Both  X-  and  K-band  klystrons  (~10  kmc  and 
'''2h  kmc  respectively)  were  used.  The  K-band  klystron  was  needed  for  the 
measurements  since  the  zero  field  splitting  is  98  kmc.^  Both  cavities  were 
a  silvered  ceramic  so  that  100  kc  modulation  of  the  DC  magnetic  field  could 
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be  enqployed.  X-ray  irradiations  were  made  for  approximately  15  minutes  in 
the  white  beam  of  a  GE  X-ray  machine  operated  at  1+5  kvp-40  ma  with  the  crys¬ 
tal  approximately  1  inch  fiom  the  tungsten  target.  Magnetic  field  measure¬ 
ments  could  be  made  accurate  to  0,2  gauss  and  0.1  gauss  for  K-  and  X-band 
frequencies  respectively.  Klystron  frequency  measurements  of  5  Mc/sec  and 
0.5  Mc/sec  accuracy  for  K-  and  X-bands  respectively  gives  a  lower  limit  of 
±0,0005  on  the  error  of  the  g-value. 
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III.  EXPERIMENTAL  RESULTS 


A.  CRYSTAIS 

From  the  spin  resonance  measurements  of  the  VO'*'^  crystals  the  following 
axial  spin  Hamiltonian*  constants  were  calculated: 


Present  Work 

Other  Results 

( Powder) 

Sll 

=  1.9328  ±  5 

1.947^^ 

1.948^^ 

1.88  -  1.93^^ 

Si 

=  1.9802  ±  5 

1.988 

1.987 

1.978  -  1.989 

|A| 

=  0.01824  ±  2/cm 

0.0158 

0.0159 

0.0158  -  0.0184 

|b| 

=  0.007162  ±  5/cm 

0.0054 

0.0052 

0.0061  -  0.0075 

IQ'  1 

=  0.00024  ±  5 /cm 

The  VO'*’^  radical  enters  the  crystal  with  the  VO'*’^  axis  in  one  of  four  different 

+P 

directions,  as  shown  in  Fig.  7.  A  pair  of  VO  axes  lying  in  the  same  plane 
belong  to  the  same  type  of  sites  (the  crystal  has  two  molecules  per  unit  cell), 
and  the  lower  concentration  position  (10^)  has  slightly  different  g-values 
(g,^  =  1.931^  -  5>  gj_  =  1.9812  ±  5)  a  line  width  of  4.7  gauss  compared  with 
5.7  gauss  for  the  40^  position.  No  variation  in  intensity  ratio  was  noted  for 
crystals  of  different  VO/Zn  concentration. 

* 

i'^axial  ~  ®j|  t  H^Sx+HySy  ]  +  AI2S2 

+  B(lxSx+IySy)  +  Q' [if  -  j  1(1+1)  ] 
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Fig.  7*  Experimental  positions 


FOUR  EQUIVALENT  V0*<  POSITIONS 
IN  A  REGULAR  OCTAHEDRON ,  THE  VO*^' 
AXES  POINT  TOWARDS  THE  CENTROIDS 
OF  THE  TRIANGULAR  FACES 


V0(NH4)2(S04)2-6Hb0  (1^  VD/Zn) 


The  angle  between  the  plane  formed  by  the  K1K3  axis  ( see  Fig.  7)  an<i  the 
c-axis  was  found  to  be  0.6®  ±  0.3°.  A  similar  measurement  for  the  crys¬ 
tals  could  not  be  made  since  it  was  found  that  the  yz-plane  of  the  one  site 
made  an  angle  of  ~6°  with  the  yz-axis  of  the  other  site. 

1.  Discussion 

The  axial  spin  Hamiltonian  constants  (gj,,  g^j,  |a|,  |b|,  and  |Q' ])  ob¬ 
tained  from  measurements  at  0  and  90"  fit  the  angular  part  between  0  and  90° 
to  a  remarkable  degree,  as  shown  in  Table  I.  In  Fig.  7  the  z-axis  for  the 
analysis  of  the  data  is  the  axis;  in  fact,  it  is  the  location  of  the  z- 

axis  and  the  result  of  g^^  >  (see  Fig.  3)  that  indicates  the  vanadyl  rad¬ 
ical  is  present  as  an  entity  in  the  crystal  and  that  its  axis  is  oriented  in 
the  directions  shown. 

One  can  attempt  to  explain  why  there  are  four  directions  for  the  VO'*'^ 
radical  and  only  two  molecules  per  unit  cell  on  the  following  basis.  The  VO"^^ 
ion  is  expected  to  enter  an  octahedron  of  water  molecules  in  one  of  four 
equivalent  directions  (see  Fig.  8)  in  order  to  achieve  minimum  energy.  How¬ 
ever,  in  the  Tutton  crystal  the  water  molecules  do  not  form  a  regular  octa¬ 
hedron  or  even  one  with  axial  symmetry,  but  rather  one  with  considerable 
rhombic  asymmetry.  This  is  seen  from  the  X-ray  data  and  from  the  gpp 
data  presented  in  Section  III-C.  This  rhombic  symmetry  makes  the  four  posi¬ 
tions  of  the  ion  non-equivalent  and  ordered  in  energy.  That  is,  one  of 

the  positions  (40^  -  see  Fig.  7)  is  energetically  more  favorable,  another 
(10^  -  see  Fig.  7)  is  energetically  less  favorable,  and  the  remaining  two 
positions  are  so  energetically  unfavorable  that  the  small  population  of  these 
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TABLE  I 

COMPARISON  OF  EXPER I  MENTAL  AND  CALCULATED  VALUES  OF 
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axes  In  regular  and  distorted  octahedra. 


positions  are  unobservatle  by  EPR  techniques.  The  assignment  of  a  10^  posi¬ 
tion  and  a  40^  position  to  one  of  the  two  types  of  sites  can  be  made  by  look- 
ing  at  the  positions  of  the  VO  ions  in  the  water  octahedra  and  noting  the 

energetically  favored  position.  This  assignment  is  then  not  arbitrary  to  the 

+2 

extent  that  the  X-ray  data  are  assumed  to  be  correct  and  that  the  VO  ion 
distorts  the  water  octahedra  a  negligible  amount. 

One  Interesting  result  lies  in  the  value  of  |Q' |.  Though  evaluated  at 
90®  where  its  effect  is  of  the  order  of  one  gauss,  the  term  involving  |Q' | 
Increases  to  4  gauss  at  6o°.  And  since  the  difference  between  calculated  and 
experimental  values,  as  shown  by  Table  I,  is  of  the  order  of  one  gauss,  this 
value  for  |Q' |  should  be  quite  close  to  being  correct.  Now  |Q' |  is  related 
to  the  quadrupole  moment  of  the  nucleus,  Q,  by  the  relation 


3eQ  S^V 
4I( I+l) 


where 


5z2 


sSz^/Qjys-tal 


+  <  l/r3  > 

'  electron 


Evaluating  this  expression  for  V^^,  (eQ  =  0.2  x  10"^^  cm^) 
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+  <  l/r^  >  1  =  5  X  10^^/cm3 

\8zVcrystal  '  electron  ^  ' 

The  value  of  <  l/r^  >,  if  evaluated  by  the  expression 


<  l/r^  >  = 


^  TPnP 
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where 


1)  =  (A-B)/5 

and 

y  =  Mn/Pnl 

is  1.5  X  10^5/cjii^,  Thus  it  appears  as  if  (S^v/dz^)  is  of  opposite  sign 

ci*ys"t8.x 

to  <  l/r^  >  and  equal  in  magnitude  to  three  places; 

Another  Interesting  experimental  result  is  the  splitting  of  each  of  the 
resonances  into  five  separate  resonances  when  the  magnetic  field  is  in  the 
direction  of  the  tetragonal  axis  of  the  surrounding  water  molecules  ( see  Fig. 
9) .  This  splitting,  approximately  4  gauss,  may  provide  some  information  on 
the  molecular  orbital  wave  function  of  the  vanadyl  ion  and  the  surrounding 
water  molecules. 

B.  POWDER 

•  p 

The  EPR  spectmra  obtained  from  the  powder  formed  by  crushing  the  VO 
crystals  gives  the  same  qualitative  spectrum  (see  Fig.  4a)  that  is  reported 
in  Refs.  9”13»  One  should  be  able  to  explain  this  spectrum  quantitatively 
from  the  spin  Hamiltonian  constants  of  the  single  crystal.  Since  the  result¬ 
ing  crystal  spectrum  shows  only  axial  symmetry,  the  spectrum  is  dependent 
only  on  0,  the  angle  between  the  magnetic  field  and  the  VO^^  axis.  Thus,  for 
a  given  magnetic  field,  only  those  crystals  or  crystallites  whose  VO^^  axes 
are  at  an  angle  0  to  the  magnetic  field  such  that  H  =  %(©)  will  contribute 
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to  the  absorption.*  Because  of  the  axial  symmetry,  the  number  lying  at  this 
angle  is  proportional  to  the  solid  angle  sinG  dG.  Thus  the  number  of  VO"''^ 
ions  contributing  to  paramagnetic  resonance  between  magnetic  field  H  and 
H+dH  is 

-7/2 

y  SinG 

Jr/2 

where  the  sum  is  over  the  eight  values  of  I,  the  nuclear  spin. 

Figure  4b  shows  a  plot  of  such  an  equation  evaluated  for  the  spin  Hamil¬ 
tonian  parameters  for  a  single  crystal  containing  VO'*’^.  Note  that  in  this 
figure  the  experimental  curve  is  the  first  derivative  of  the  absorption  curve 
and  that  it  shows  the  effect  of  the  finite  width  of  the  resonances.  The 
sharp  peaking  of  the  calculated  curve  corresponds  to  those  VO'*’^  oriented  at 
90°  and  corresponds  also  to  a  (cosG)”^  dependence.  From  the  powder  spectrum, 
* 

2g2p2Ho  \  K  /  \g2/  K^gSp 

X  A^B^g^Sgfm [41(1+1)  -  8m2  -  1]  -  i  (Q')^ 

X  [21(1+1)  -  2m^  -  l]m 

where 

=  g^A^cos^Q  +  g^B^sin^G 
g  =  g^cos^G  +  ^sin^G 
Ho  =  hv/ gP 
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the  separation  of  the  two  extreme  peaks  corresponds  to  7|A|/g||P,  and  if  one 
neglects  the  quadrupole  contribution,  the  separation  of  the  two  extreme  lines 
in  the  central  part  of  the  spectrum  is  7|B|/g^P  so  that  |a|,  |b|,  g^,  and  gj_ 
can  be  obtained  from  a  powder  sample — to  the  accuracy  limited  by  the  width 
of  the  peaks. 


C.  IRRADIATED  VO'^^  CRYSTALS 

Operation  at  K-band  frequenc ies  gave  the  following  values  for  the 
rhombic  spin  Hamiltonian  constants  (see  Appendix  B) : 


Sx 


V'*'^  ( Grown) 

g2  =  1.9718  ±  5 

gy  =  1.9750  ±  5 

|d|  =  0.15605  ±  5/cm 
|e|  =  0.02297  ±  5/ cm 
|a|  =  0.008270  ±  5/cm 

=  2°  ±  3.0'’ 

*a  20.5°  ±  0.5“ 


( Irradiated  VO'*'^) 
1.9719  ±  5 

1.91^5  ±  5 

0.15609  ±  5 

0.02503  ±  5 

0.008270  ±  5/cm 

5°  ±  1.0^ 

20.5°  ±  0.5° 


(Reference  2) 

1.951  ±  2 

0.158  ±  10 
0.0^9  ±  J+0 
0.0088 
2” 

22“ 


+P 

For  the  V  produced  by  X-ray  irradiation  in  the  three  different  Tutton  salts, 


the  following  rhombic  spin  Hamiltonian  constans  were  calculated: 


*a  and  t  locate  the  z-axis  with  respect  to  the  ac-plane  and  the  c-axis  (see 
Fig.  7). 
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Zn(NH4)2(S0j2-6H^0 

Mg(NH4)2(S04)2-6H£0 

ZnK2(S04)2'6H20 

g2  =  1.9718  ±  5 

1.9720  ±  5 

1.9722  ±  5 

=  gy  =  1.97^5  ±  5 

1.9733  ±  5 

1.9750  ±  5 

|d|  =  0.15609  ±  5/cm 

0.15795  ±  5  (+1.2^) 

0.15246  ±  5  (-2.5^) 

|e|  =  0.02303  ±  5/cm 

0.02456  ±  5  (-6.6^) 

0.02746  +  5  (+19.2^) 

|a|  =  0.008270  ±  5/cm 

0.008270  ±  5/cm 

0.008270  ±  5/Gm 

t  =  3”  ±  1“ 

1°  ±  1° 

11°  ±  2° 

a  =  20.5°  ±  0.5“ 

20.0°  ±  0.5° 

i4.7°  ±  0.5° 

The  angle  that  the  K1K3  plane  ineikes  with  the  c-axis  could  not  he  detemlned 
since  the  zy-plane  of  the  one  site  makes  an  angle  of  6°  with  the  zy -plane  of 
the  other  site. 

1.  Discussion 

As  seen  hy  the  close  agreement  of  the  V'*’^  (grown)  and  V'*’^  (irradiated 
VO"'’^)  spin  Hamiltonian  constants,  and  the  angle  a,  one  can  assume  that  when 
the  oxygen  moves  away  from  the  VO'^^  site  leaving  behind,  it  is  removed 
far  enough  so  that  its  influence  on  the  crystalline  field  is  no  longer  felt. 
This  is  also  implied  hy  the  data  on  the  irradiated  VO’*'^  in  the  Tutton  salts 
where  Mg  is  substituted  for  Zn  and  K  for  NH4.  That  is,  these  ions,  though 
several  angstroms  away  from  the  site  airi  its  associated  octahedron  of 
water  molecules,  distort  the  octahedron  by  an  amount  that  can  be  measured  by 
the  changes  in  the  values  of  |d|,  jE],  a,  and  \|f.* 

*In  experiments  on  Cu’*'^  in  vai-ious  Tutton  salts,  Bleaney,  Penrose  and  Plump- 
ton^9  noticed  that  the  changes  in  the  monovalent  ion  Influence  the  angles 
a.  and  \|r  to  a  great  extent.  They  also  noticed  large  changes  in  the  g-values 
(~1.0^)  which  were  not  observed  in  the  present  set  of  experiments. 
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If  one  performs  a  general  coordinate  rotation  on  the  spin  Hamiltonian 
for  rhombic  symmetry  (in  order  to  diagonalize  the  Zeeman  term — Appendix  A), 
it  is  seen  that  only  when  the  magnetic  field  is  along  the  x-,  y-,  or  z-axes 
do  the  cross  terms  become  zero,  simplifying  the  spectrum  and  enabling  one  to 
make  measurements  that  can  be  used  to  determine  the  spin  Hamiltonian  con¬ 
stants.  For  this  reason  no  attempt  was  made  to  measure  the  angular  part  of 
the  spectrum.  An  illustration  of  the  effect  of  this  mixing  can  be  seen  in 
Fig,  10.  Note  that  in  both  the  x-  and  y-directions  the  mixing  and  the  result¬ 
ing  spread  of  the  spectrum  from  the  second  site  is  so  complete  that  it  is  al¬ 
most  unobservable. 

The  |d|  value  can  be  calculated  from  the  data  along  the  z-axis  and  the 
|e|  value  from  that  along  the  x-  or  y-axis.  However,  the  plane  of  the  zy- 
axis  from  the  one  site  does  not  coincide  with  that  from  the  second  site.  For 
this  reason  the  center  part  of  the  spectrum  in  Fig.  1  showing  the  magnetic 
field  along  the  Ka^stxis  is  not  a  "pure"  spectrum.  The  Kg-axis  is  approxi¬ 
mately  5"  from  the  x-axls  of  each  of  the  two  sites. 

Note  that  the  value  of  |e1  reported  by  Bleaney,  Ingram  and  Scovll^  is 
almost  twice  the  value  found  in  these  experiments.  Since  they  used  the  in¬ 
formation  along  the  z-direction  to  evaluate  both  |d1  and  |e|,  and  |e|  enters 
as  a  second-order  effect  along  this  direction,  it  is  thought  that  their 
value  13  not  as  correct  as  the  value  reported  here. 


28 


29 


IV.  FURTHER  STUDY 


The  additional  structure  of  the  EPR  absorption  lines  shown  in  Fig.  9 
might  yield  some  information  about  its  origin  if  attacked  with  double  res¬ 
onance  (ENDOR)  techniques. 

Furthermore,  with  ENDOR  techniques  the  size  of  the  quadrupole  interac¬ 
tion,  which  may  be  present  to  some  extent  as  suggested  by  EPR  experiments, 
can  be  more  accurately  determined. 

Since  second-order  effects  indicate  that  A  and  D  have  the  same  sign,  a 
low  temperature  experiment  to  Indicate  the  relative  intensity  of  the  high  and 
low  field  transitions  should  determine  the  absolute  sign  of  D  and  thus  A. 

The  growing  of  both  V'*"*’  and  in  the  same  crystal,  though  the  ini¬ 

tial  attempt  has  been  unsuccessful,  would  serve  as  a  very  precise  indicator 
that  the  environment  of  the  (grown)  and  the  V*"*'  (produced  by  X-ray  ir¬ 
radiation)  are  identical. 

Similar  experiments  should  be  carried  out  to  determine  the  optical  ab¬ 
sorption  spectrum  for  determination  of  the  spin  orbit  coupling  constant  \  and 
to  measure  the  production  of  V"'"'’  ( irradiated)  as  a  function  of  the  radiation 
dose  to  the  crystal. 
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APPENDIX  A 


CALCULATION  OF  ANGULAR  VARIATION  OF  EPR  RESONANCE  IN  RHOMBIC  FIELD 

In  this  appendix  we  wish  to  obtain  formulae  so  that  evaluation  of  the 
constants  in  the  rhombic  spin  Hamiltonian  can  be  made  from  experimental  data. 
We  start  with  the  general  rhombic  spin  Hamiltonian 

^  rhombic  —  PCSz^z^z  Sx®x®x  Sy^^rSy]  ^[sf  “  ^  S(S+l)  ] 

+  E(S^-Sf)  +  AIzSz  +  BSxIx  +  CSyly  +  Q' [if  -  |  l( I+l) ] 

+  Q"[li  -  If] 

and  proceed  to  diagonalize  the  Hamiltonian  in  the  Zeeman  term  and  subsequently 
treat  the  off-diagonal  elements  as  perturbations. 

1.  FINE  STRUCTURE  TERM 

Since  we  would  like  to  have  gxHxSx+gy%Sy+gzHzSz  =  gS^H  we  proceed  to 
make  a  coordinate  transformation  with  the  new  set  of  coordinates  x'  y'  z'  re¬ 
lated  to  X  y  z  by  the  Eulerlan  angles  0,  and  \lf, 


y 
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then  Sx,  Sy,  and  Sz  are  related  by  the  following: 


Sjj  =  (cos  9  cos  \|f  -  sin  9  sin  \lf  cos  )Z^)Sx  -  (cos  9  sin  \|f 

+  sin  9  cos  cos  i6)S^  +  sin  9  sin  Sz 

Sy  =  ( sin  9  cos  \lf  +  cos  9  sin  \|;  cos  fi)  Si  -  (  sin  9  sin  \|f 

+  cos  9  cos  \Jr  cos  fi)Sy  -  cos  9  sin  j6  Sz 

S2  =  sin  ^  sin  \|/  Si  +  sin  ^  cos  \lf  S^-  +  cos 


Now  if 

~  H  sin  0  cos  6 
Hy  =  H  sin  0  sin  6 
Hz  =  H  cos  0 


and  if  we  require 


where 


then 


cos  9 
sin  9 
sin  ^ 
cos  f6 


-(gy/gx)sin  & 
(?x/8_i)cos  6 

(gj_/g)sin  0 
(gz/g)cos  0 


g|  =  gfcos^b  +  g^sin^S 
g^  =  g^cos^©  +  g^sin^© 


gz^zSz  gx^xSx  Sy%Sy  -  gHSz 
Thus  the  Zeeman  term  is  diagonalized. 

The  axial  and  rhombic  field  terms  become,  respectively, 
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The  diagonal  terms  of  the  fine  structure  part  involving  Sz,  Sf,  and  S+sl+ 
SiS-J-  evaluated  for  the  states  |s  Ms  >  with  the  selection  rule  M  ->•  M-1  yields 


D-3E  cos^S 

L  Vgf 

The  non-diagonal  terms,  i.e.,  the  terms  involving  (S+)^,  (sL)^ ,  (slSz+ 
SzSl)  and  (S+Sz+SzSi),  can  be  evaluated  by  second-order  perturbation  theory 
(there  are  no  first-order  terms),  give  the  following  correction  to  the  rhombic 
fine  structure  term; 


US(S+1) -24M(M-1)  -9 
2gPHo 


D  - 5  (gfcos26-gSsin26) 

-  &±  ^  _ 


sin^0  cos^0 


2  2 

+  E2sin20  sin226 

gfg^ 


2S(S+1)  -6m(M-1)-3 

8gPHo 


1“  g2  2 

D  “f  sin^  +  (gfcos^S 

L  g  gf 


(l  +  ^  008^0') 

4.  f^E  QOS  0  sin  26 

2"' 

\  g  y_ 

L  tfg  J 

s. 

2.  HYPERFINE  TERM 


In  attempting  to  carry  out  a  similar  procedure  on  the  rhombic  hyperfine 
terms  the  results  become  extremely  complicated.  Since  the  experimental  re¬ 
sults  indicate  little,  if  any,  rhombic  symmetry  in  the  hyperfine  terms  for 
VO"^  and  none  for  we  will  assume  that  3  =  C  and  Q"  =  0;  i.e,,  we  have 


AIzSz  +  B(lxSx+IySy)  +  Q"[lf  -  i  l(l+l)  ]  +  •  I 


Here  we  perform  a  rotation  on  the  nuclear  spin  as  we  did  previously  on 
the  electron  spin,  with  the  difference  that  where  the  rotation  matrix  for  the 
angle  ^  was 


1 


0 


0 


0  cos  ^  -sin 

0  sin  ^  cos  ^ 


it  now  becomes  (for  ly  and  I^) : 


where  K  is  defined  by  setting  the  dterminant  of  this  matrix  equal  to  unity, 
i.e.,  =  A^cos^^+B^sin®^.  Since  we  have  chosen  axial  symmetry  in  the  hyper- 

fine  tem  we  may  arbitrarily  choose  f  =  0.  The  hyperfine  term  then  becomes 


AIzSz  +  B( IxSx+IySy)  =  KIzSz  +  BSxIx  +  ^  Syly 

The  first  term  is  the  only  one  which  is  diagonal  in  M;,  m,  and  with  the 
section  rules  M  -►  M-1,  Am  =  0  its  contribution  to  the  separation  of  energy 
levels  is  fta.  Evaluating  the  other  terms  by  second-order  perturbation  theory 
(neglecting  the  terms  of  order  IGn  with  respect  to  gPH ),  we  obtain 


Bf  /A^+K^\[ 1(1+1) -m^]  2 

^  \  K?  J  gpH 


_ m 

gPH 


3.  QUADRUPOLE  TERM 

Since  the  quadrupole  term  Q' [if  -  ^  l(l+l) ]  transforms  similarly  to  the 
fine  structure  part,  one  can  quickly  obtain  the  quadrupole  contribution  to  the 
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difference  in  energy  levels  as 


6±gf(  sin^Bcosge 
g2  2KM(M-1) 


[4l(l+l)  -  8m2  -  1] 


sln^0 

8KM(M-1) 


[21(1+1)  -  2m^  -  1] 


These  are  also  obtained  by  second-order  perturbation  theory  and  then  by 
applying  the  selection  rules  M  -►  M-1,  Am  =  0.  Although  there  is  a  diagonal 
contribution  of  the  quadrupole  part  to  the  energy  level  of  the  form 


IS 


cos^^  m2  +  — ^  sin^^  [l(l+l) 


m2]  -  -  1(1+1) 


it  is  the  same  for  each  value  of  M,  so  that  the  selection  rule  Am  =  0  results 
in  cancellation  of  this  term. 
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APPENDIX  B 


EXACT  ENERGY  LEVEL  CALCULATIONS  FOR  V+2  ALONG 
X-,  y-,  AND  z-AXES  FOR  RHOMBIC  FIELD 

Given  the  rhombic  spin  Hamiltonian  with  axial  hyperfine  term 


iky. 


=  Pl'i'S  +  D[Sf  -  I  S(S+1)  ]  +  E(S|-Sf)  +  I-A-S  , 


Rhombic 

we  obtain  for  the  fine  structure  term  for  V'*’^  (S  =  3/2)  and  for  9V- parallel 
to  z-axis 

=  gzpHSz  +  D[sf  -  i  S(S+1)  ]  +  I  [S?  +  S?]  . 


Degenerate  perturbation  calculations  lead  to  the  following  secular  deter¬ 
minant  : 


M  = 

-3/2 

-1/2 

1/2 

3/2 

3 

"  2 

-  J  gzPH+D 

-X 

V?  E 

_  1 

-  1  g,.PH-D 

^/3  E 

2 

-X 

1 

2 

V3  E 

1  SzPH-D 

-X 

3 

2 

■S  E 

1  gzPH+D 
-X 

This  determinant  may  be  solved  for  the  following  energy  levels: 
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Sz  =  5  gzPH  +  ^3E2+(gzPH+D)‘ 


i.  x"*" 

2  J 


S,  =  - 


b  ^ 


=  -  gzPH  +  V3E2+(g2PH-D)' 


=  -  gzPH  -  \/3E2+(gzPH+D)‘ 


2  j 


-  i  gz^H  -  N/3E2+(g„PH-D)‘ 


For  the  transitions  M  ->•  M-1  the  change  in  energy  is; 


23/2,1/2  -  ^  -^1  -  gz^Ha/a  +  '^32^+( gzPHa/a+D) ^  -  s/3E^+( gzPHa/a-D)* 


=  X1-X2  =  -gzPHi/a  +  'j'3E2+(gzPHi/2-D)2  +  N/3E2+(gzPHi/2+D)‘ 


E-1/2,.3/2  =  \2-Xi  =  -  '^3E2+(gzPH.i/2+D)2  +  N/3E2+(gz&H.i/2-D)‘ 


Making  use  of  a  Taylor's  expansion  for  each  of  the  square  roots  and  solv¬ 
ing  for  the  magnetic  field  (in  the  major  term)  yields  the  following  equation 
(to  fourth  order) : 


( A)  H3/2  =■-  —  -  2  ( 

'  gzP  \gzP/  (gzP)^Hf 


HS/2 


U  -  (— ^)"J 

'  —  Ci  TT  !  • 


‘gzPH3/2 


SI  g*D  1 

^  (SzP)^  (Hg/a)^ 


hv 


(B)  H1/2  = 

/  ^  Q 


3E2 


gzP  ( SzP) ^Hi/2 


1  - 


SzPHi/a 


Er*D 


(gzP)^Hl/2 


1_  g^PH  ^  Vg^PHi/a/  _ 
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t 


(c)  H.i/2 


3E^D 

(gzP) 


( 


-)■ 


^  /  1  Y 

^  (gzP)^  \H.i/2/ 


For  the  magnetic  field  parallel  to  the  x-dlrectlon,  the  Zeeman  term  Is 
not  diagonal.  To  diagonalize,  a  coordinate  rotation  about  the  y-axls  by  90° 
Is  performed,  resulting  In  the  fine  structure  spin  Hamiltonian 


gxPHS^  +  D 


1 

4 


[(s;)^  +  (si)^  +  sis:  +  sisi] 


5 


+ 


[(si)^  +  (Sl)^ 


sis; 


which  leads  to  the  following  secular  determinant: 


M  = 

-3/2 

-1/2 

1/2 

3/2 

_  3 

2 

-  1  fePH 

r  D+E 

s  — 

1 

2 

-  \  SxPH 

1 

2 

V3^ 

2 

1  SxPH 

2 
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This  is  the  same  secular  determinant  for  parallel  to  the  z-axis  with  -► 
-(D-3E)/2  D,  and  (d+E)/2 E.  Thus  with  these  changes,  Equations  (A),  (b), 
and  (C)  also  hold  for  the  x-axis.  Similarly,  for  the  magnetic  field  parallel 
to  the  y-axis  the  following  substitutions  are  made : 


gx  ^  gz.  -(^^)  ^  D,  and  ^  ^  E  . 


D-E 
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PART  III 


Mn++  EPR  RISULTS  IN  AnByi  COMPOUNDS 


^>y 

G.  H.  AzarBay-ejani 


I 


I.  INTRODUCTION 

In  earlier  reports  on  ESR  measurements  in  AjjByj  compounds  doped  with 
Mn  ^  some  Interesting  results  on  the  ground  state  splitting  factor  Ja  have 
been  obtained.  Further  experiments  which  have  been  carried  out  pertaining 
to  the  cubic  ground  state  splitting  will  be  described  here.  In  Section  II 
crystal  preparation  and  doping  with  impurities  of  CdTe  and  ZnTe  will  be  con¬ 
sidered.  In  Section  III  ESR  measurements  will  be  reported. 
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II.  CRYSTAL  PREPARATION^ 


CdTe  axid  ZnTe  crystals,  which  we  have  used  in  our  experiments,  are  made 
from  the  elements  in  fused  silica  quartz  tubes  under  vacuum  and  heated  in  R-F 
furnaces.  The  following  steps  are  taken  in  preparing  these  crystals. 

A.  PREPARATION  OF  QUARTZ  TUBES 

All  sample  tubes  are  made  of  10  mm  (ID)  tubing  with  4-10  in.  of  added 
tubing.  In  one  end,  these  tubes  are  closed  in  regular  round  shape,  and  in 
the  other  are  connected  to  a  quartz  tube  stub  (see  Fig.  1).  The  empty  part 
is  4-10  in.  long. 


SLAB 


SAMPLE 


Fig.  1.  Tube  preparation. 


The  tubes  are  checked  for  leaks  by  using  the  vacuum  station  and  Tesla  coil. 
Good  tubes  are  cleaned  in  the  following  sequence: 

(1)  Benzene  rinse 

(2)  Acetone  rinse 

(5)  Tapwater  rinse  with  brushing  at  each  stage 

(4)  Warm  concentrated  nitric  acid  bath  for  2-4  hr 

(5)  Deionized  water  rinse  (20  times) 

(6)  Drip  drying  by  holding  the  tubes  vertically  with  the  open  end  down 
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(7)  Heating  the  tube  at  quartz  fusing  temperature  so  that  any  foreign 
impurity  remaining  on  the  walls  of  the  quartz  tube  will  he  fused 
into  the  quartz 

(8)  Cooling  by  air  flow 

(9)  Closing  the  tube  with  a  cork  wrapped  with  glassine  powder  paper 

(10)  Labeling  the  tube  prepared  in  this  fashion  with  the  letter  "C, "  in¬ 
dicating  that  it  is  clean  for  use  in  crystal  preparation. 

B.  PREPARATION  OF  SAMPLE 

The  following  steps  are  taken  for  weighing  the  saii5)le. 

(1)  Recording. — On  a  special  data  sheet  the  elements,  properties,  and 
weight  of  each  are  recorded.  The  base  elements  such  as  Cd  and  Te  in  CdTe 
should  have  a  purity  of  99^999i>>  whereas  in  the  elements  of  intended  impurity 
such  as  Mn  or  Cu  it  can  be  99.99^- 

( 2)  Treatment  of  the  Elements. — Before  weighing  Cd  and  Zn,  the  metal 
bars  are  cut  into  pieces  less  than  10  mm  long  and  are  etched  with  concentrated 
nitric  acid.  Then  they  are  rinsed  thoroughly  with  deionized  water,  covered 
with  high-purity  metanol,  agitated,  drained,  and  allowed  to  dry  on  paper 
towels.  Tellurium  is  usually  used  without  etching  but  always  in  pieces  so 
that  its  oxide  content  will  be  small. 

(5)  Weighing  of  the  Elements . — The  amount  necessary  for  each  element  is 
calculated  by  desk  calculator  up  to  7  figures  and  then,  with  use  of  the  ele¬ 
ments  etched  as  described  in  paragraph  (2),  the  weighing  is  carried  out.  The 
elements  belonging  to  Group  (Cd,  Zn)  must  be  weighed  before  those  belong¬ 
ing  to  Group  and  Group  VI;  therefore  for  CdTe,  Mn  Cd  is  weighed  first 
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and  imnediately  added  to  the  tube,  with  tellurium  and  Mn  then  weighed  and 
added  to  the  tube. 

(4)  Evacuation  and  Seal-Off  of  Sample . — The  stube  in  Fig.  1  connects 
the  tube  to  a  vacuum  station  which  has  a  fore  pump  for  initial  evacuation  and 
a  diffusion  pump  for  further  evacuation.  After  the  desired  vacuum  is  achieved 
the  tube  is  heated  in  the  "c"  region  to  a  length  of  1 — l-l/2  in.  with  about 

2  in.  above  the  top  of  the  sample  so  that  the  walls  of  the  tube  collapse  com¬ 
pletely  for  a  satisfactory  seal-off. 

(5)  Reaction  of  Elements  in  R-F  Furnaces. — To  ensure  a  uniform  distri¬ 
bution  of  the  impurity  ions  in  the  compound  the  tube  is  placed  vertically  in 
a  graphite  cylinder.  This  graphite  cylinder  is  placed  along  the  axis  of  a 
vertical  solenoid  of  l/8-in.  (ID)  copper  tube  connected  to  a  20-KW  Lepel  R-F 
generator.  The  temperature  is  raised  above  the  crystalline  melting  point  and 
held  there  for  several  minutes;  then  the  power  is  turned  off,  which  in  5"10 
minutes  lowers  the  temperature  to  room  temperature.  The  quartz  tube  of  the 
sample  is  then  Inspected  for  cracks.  If  single  crystals  are  desired,  the 
polycrystalline  obtained  from  the  R-F  furnace  is  annealed  for  several  weeks 
under  a  temperature  slightly  lower  than  the  melting  point.  The  same  technique 
is  used  for  both  ZnTe  and  CdTe  doped  with  Mn. 
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III.  ESR  MEASUREMENTS 


15 

The  first  experiment  on  AjjByi  compounds,  performed  by  Van  Wieringen'^  on 
ZnS :Mn,  raised  a  great  demand  for  similar  ESR  experiments.  Since  then  many 
other  materials  with  structures  close  to  cubic  have  been  investigated.  In 
Table  I  some  of  the  recent  ESR  results  on  Mn'*'"*'  which  have  been  obtained  in 
our  laboratories  or  by  other  investigators  are  given. 


TABLE  I 

Mn'^+  ESR  EXPERIMENTAL  RESULTS 


Material 

Crystalline 

Structure* 

Coordination 

So 

(A“) 

Co 

(A') 

g 

3a 

10"4cm”i 

D 

lO'^cm"^ 

Ref. 

MgO 

R 

6 

4.24 

55 

0 

a 

CaO 

R 

6 

4.81 

17.7 

0 

b 

aZnS 

Z 

h 

27.7 

0 

a 

pZnO 

W 

4 

5.24 

5.18 

6 

-216.9 

c 

ZnSe(  P) 

Z 

4 

5.67 

2.010 

- 

- 

d 

ZnTe 

Z 

4 

6.12 

- 

a 

CdS 

W 

4 

4.15 

6.64 

4.2 

8.2 

c 

CdSe 

W 

4 

4.50 

7.01 

47 

? 

a 

CdTe 

Z 

4 

83.1 

- 

a 

*W  stands 

for  wurtzite  ( BZnS) ,  R  for 

rock  salt  (Nad) 

,  and 

Z  for  zincblende  ( ZnS) 

structure. 

Refs. 

a.  Results  of  our  laboratories. 

b.  Shuskes,  Phys.  Rev.  I27,  I529  (I962). 

c.  Doraln,  Phys.  Rev.  2^,  IO58  (I958). 

d.  Matsumara,  J.  Phys.  Soc.  Japan  14,  I08  (1959). 
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£. 

According  to  Watanabe® 


5a  Qc  ( Dq)  ^ 

(1) 

Dq  ccr'^/^ 

(2) 

where  r  Is  the  position  of  an  electron  with  respect  to  a  paramagnetic  ion  and 
R  is  the  position  of  the  nearest  neighbors  of  paramagnetic  ions.  This  theory 
was  originally  developed  for  ions  in  an  octahedral  field  produced  in  crystals 
with  the  structure  of  NaCl.  In  Table  I,  MgO  and  CaO  have  this  structure,  and 
if  we  assume  that  r^  of  Eq..  (2)  remains  the  same  for  both  materials  we  have: 


The  experimental  results  (Table  I)  gives 


^  _55_ 
VfcaO^  17.1 


5.1 


Also 


where 


and 


E  =  experimental 
T  =  theoretical 


(5b) 


(4a) 


(4b) 
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Comparison  of  and  gives 


where 


Also 


where 


0.312 


(5a) 


(5b) 


^;aiTe\ 

.^CdTe/g 


1.07 


(6a) 


AznTeN 

VCdTe/p 


1.75 


(6b) 


Equations  (4a)  and  (4b)  reveal  a  sharp  discrepancy  (as  much  8.s  35^)  between 
experimental  result  and  theoretical  prediction  and  indicate  a  great  need  for 
testing  other  materials  having  a  cubic  field.  This  is  why  many  studies  have 
been  carried  out  so  far.  As  reflected  by  Eqs.  (5a),  (5b),  (6a),  and  (6b)  the 
discrepancy  of  35^  in  MgO  and  ZnS  is  raised  by  1000^  in  the  case  of  ZnS  and 
ZnTe,  and  by  75^  la  the  case  of  ZnTe  and  CdTe.  This  demonstrates  that  al¬ 
though  crystalline  field  theory  gives  a  satisfactory  result  in  light  AuByi 
compounds  and  especially  in  octahedral  coordination  such  as  CaO  and  MgO,  it 
fails  to  account  for  heavier  compoiinds  belonging  to  this  group,  especially 
those  in  tetrahedral  coordination.  In  the  following  sections  the  experiments 
carried  out  on  ZnTe  and  CdTe  doped  with  Mn  will  be  considered. 


A.  EXPERIMENTAL  METBDD 

The  measurements  at  3OO,  78,  and  4.2°K  were  carried  out  in  an  X-band 


9 


(9. 2-9. 4  kmc/sec)  magnetic  resonance  spectrometer.  The  ZnTe  single  crystals 
containing  0.1^  MnTe  were  prepared  as  discussed  in  Section  II-B.  Being  Brittle 
and  of  low  resistivity,  the  ZnTe  single  crystals  are  very  difficult  to  cleave 
and  orient  properly.  Therefore  we  had  to  select  many  samples  and  design  a 
new  cavity  which  permitted  the  simultaneous  and  orthogonal  rotation  of  Both 
DC  magnetic  field  and  crystal  in  the  cavity.  This  cavity  is  shown  in  Fig.  2. 
The  teflon  pulley  and  nylon  string  are  used  particularly  to  facilitate  the 
rotation  of  the  crystal  at  low  temperatures.  With  this  device  we  were  aBle 
to  orient  the  crystal  as  well  as  find  the  temperature  dependence  of  the  ESR 
parameters  in  ZnTe  and  other  similar  materials.  Considering  that  cuBlc  ZnS 
and  ZnTe  have  exactly  the  same  crystalline  structure  and  that  the  meiximum 
splitting  Between  the  fine  structure  components  of  Mn"*”*"  ESR  spectrum  in  the 
cuBlc  ZnS  occurs  when  the  magnetic  field  and  [OOl]  axis  are  coincident,'^  one 
simply  looks  for  the  maximum  separation  in  the  spectrum  of  Mn’^’^'  as  a  function 
of  angles  9'  and  These  are  the  angles  through  which  the  DC  magnetic 

field  and  crystal  in  the  cavity  are  rotated  from  arbitrary  positions.  After 
maximum  splitting  is  achieved  and  hence  the  location  of  one  of  the  crystalline 
axes  is  determined  further  rotations  regarding  this  axis  are  made  to  ensure 
the  correct  assignment  of  the  axis.  Before  making  an  ESR  measurement,  it  is 
of  great  importance  to  ensure  that  the  piece  of  crystal  is  not  polycrystal- 
llne  or  twin.  Metallurgical  microscopes  with  magnification  ranging  from  3 
to  ICX)  have  Been  used  for  closer  inspection  of  the  faces,  which  are  not  shiny. 
After  this  stage  X-ray  photographs  from  Both  powdered  and  solid  pieces  of 
crystal  are  made  to  identify  the  structure  as  well  as  the  singularity  of  the 
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pulley  ateembly, 
showing  sample 
in  place 


Fig.  2.  Rotating  8aii5)le  holder  and  resonance  cavity  used  for  measurements 


crystal.  These  steps  are  very  important  in  the  measurement  of  cubic  field 
parameters  because  the  structure  of  the  crystal  depends  a  lot  on  the  method 
by  which  they  are  prepared.  Fortunately,  in  the  case  of  the  tellurides  of 
Zn  and  Cd  the  structure  is  zincblende  alone.  For  the  selenides  of  these  el¬ 
ements,  however,  structures  of  sphalerite  (cubic  ZnS)  and  of  wurtzite  (hex¬ 
agonal  ZnS)  are  observed,  and  sometimes  it  is  possible  for  both  structures  to 
be  present.  A  proton  magnetometer  and  Hillard-Packard  frequency  counter  have 
been  used  to  obtain  the  line  positions.  A  hydrazyl  marker  is  used  for  g  meas¬ 
urements  due  to  the  fact  that  in  each  measurement  there  might  be  slight  shifts 
in  klystron  frequency,  which  is  an  especially  important  factor  for  ten^jerature- 
dependence  measurements. 

B.  EXPERIMENTAL  RESULTS 


1.  ZnTe;Mn 

In  Fig.  3  the  spectrum  at  500°K  is  shown.  At  the  top  is  the  spectrum 
corresponding  to  0  =  0,  and  at  the  bottom  the  spectrum  corresponding  to  0  = 
30°.  (0  is  the  angle  between  [OOl]  and  the  DC  magnetic  field  H. )  This  phe¬ 

nomenon  of  dependence  of  intensity  on  0  is  more  pronounced  in  the  materials 
with  larger  a  (cubic  field  splitting  parameter)  because  the  farther  the  fine- 
stinicture  components  a  and  S  are  from  the  main  line  y  ( see  Fig.  17  in  Section 
IV-D)  and  the  higher  the  temperature,  the  more  destructively  they  tend  to 
combine.  At  0  =  30°  the  separation  of  the  fine -structure  components  reduces 
nearly  to  zero;  therefore  the  five  transitions  corresponding  to  =  -3/2, 
-l/2,  1/2,  3/2,  and  5/2  occur  at  about  the  same  magnetic  field.  Figures  3a 
and  3h  clearly  manifest  the  correctness  of  orientation  assignment  giver  to 
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P/VP 


YTTTl 

(b)  9  =  SO* 

MAGNETIC  FIELD,  H 

Fig.  5,  ESR  spectrum  at  500°K:  (a)  6  =  0";  (b)  9  =  50®. 

The  spectra  for  T  =  78°  and  ^^.2®K  at  9  =  0°  are  given  in  Figs,  L  and  5 
respectively. 

2.  CdTe:Mn 

Measurements  were  made  at  5OO,  78,  and  4.2®K;  the  corresponding  ESR  ab¬ 
sorption  spectrum  for  78®K  is  shown  in  Fig.  6,  and  the  spectra  for  1+.2®K  (9 
50°  and  9  =  0°)  are  shown  in  Figs.  7  8  respectively. 
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The  superhyperflne  stinicture  corresponding  to  Cd  ions  was  observed  both 
in  CdS  and  CdTe  by  other  investigators.^ 
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IV.  THEORY 


A.  GENERAL  HAMILTONIAN  OF  Mn++ 

To  explain  the  spectra  of  Mn'*"*'  shown  in  Figs,  3-8  we  should  obtain  a 
relation  between  (l)  the  DC  magnetic  field  which  is  measured  experimentally, 
and  (2)  the  electron  spin  magnetic  quantum  number  M,  the  nuclear  magnetic 
quantum  number  m,  the  polar  angle  ©,  and  the  azimuthal  angle  Such  a  rela¬ 
tion  is  self-evident  only  if  one  derives  it  directly  from  the  following  total 
Hamiltonian  of  Mn"*"^; 


X  t  (*»+*)  -  I X » 


where 

K. 


7l  71 

Z  (Pf/2mi-Z^ri)  +  Z 

i=l  i>J=l 

Tl 

^  +  bij£^-s>5  +  Cij/-sJ) 


n 

i>J 


71 

Z  (eV2mc)(_£V2s^) -H  -  S{  L+2S)  •  H 


(7) 


(8) 


(9) 


(10) 


(11) 


^  6nPn2P 


% 

Z  |^rJ^[rf(i_^-£^)  •I+5ri- slri- 1]  +  —  6(rl)£l*I 
i=l  ^  ”  3  -  - 


(12) 
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f 


n 

6  =  — — -  y  rT^[l*2rf  -  3(ri.l)^]  with  I*  =  sJl(l+l)  (13) 

21(21-1)  Y  ”  ” 


-  XI  V(r 


(14) 


j=l  i=l 


8  -  -gNPltS-I 


(15) 


"ifti  [Eq.  (8)  ]  is  the  s\m  of  the  kinetic  energies  of  electrons;  the  coul¬ 
omb  interaction  of  the  electrons  with  the  nucleus  of  the  ion  and  each  other, 
the  TV  is  the  number  of  electrons  around  the  Mn"*"*'  nucleus;  yi  =  Z-2  =  23. 

y/ 2  [Eq.  (9) ]  is  the  spin  orbit  interaction.  The  proportionality  con¬ 
stants  a^j,  bj^j,  and  Cj^j  are  in  general  functions  of  quantum  numbers  (n^,i^, 
m^) .  The  expression  for  instance,  refers  to  the  energy  associated 
with  the  force  which  the  orbital  angular  momentum  of  the  ith  electron  exerts 
on  the  spin  of  the  jth  electron.  In  general 


Hi  I 


>  a 


ij 


It  can  be  shown  that' 


.8 


2m^c®  <i'^ 


(16) 


where  r  is  the  position  of  the  electron  with  respect  to  the  nucleus,  c. .'s 

1 J 

are  primarily  due  to  Heisenberg's  exchange  effect,  which  appear  in  the  inter¬ 
action  terms  between  the  two  electrons  i  and  j : 
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=  Kij  -  (l/2-2si.sJ)Jij 


(17) 


where 


and 


Kij  =  <  ilVijIi  >  =  <  JI'KIJ  > 


Jij  - 


(ifl) 


(19) 


Though  the  dependence  of  on  the  angle  between  the  spins  is  similar 

to  that  of  the  mutual  potential  energy  of  two  magnetic  dipoles 

r-5(r2j^i.^J  _  3jii- r^ij*r) , 
they  are  two  entirely  different  things, 

(lO)]  is  the  magnetic  dipole-dipole  interaction  and  it  can  he 
simply  derived  hy  finding  the  magnetic  field  produced  by  ith  dipole  at  the 
position  of  Jth  dipole; 

=  V^xA(rij)  =  -V^x[|iix7^(l/rij)  ] 


=  V-^x 


r  •  ^ 
ij 


=  V^x(M.^xr^j)  -  —  x(u^xr^j) 


=  r 


■ij 

-3 

Ij 


r? . 
ij 


Jr^*)  .^Ir^*) 


ij  -J 


-  -rf  -iiM  j) 


end 


K 


iJ 

3 


=  = 


r"  j  ( r|  j^l  •  ji  J -5^1 .  rl  J  •  rl  J ) 
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Let 


=  and 


2PsJ 


Then 


t 


ij 


3  =  .(r?  -  Jrij 

iJ 


and 


^  0-- 

^  ^a'  =  ^  "  3£^'^*s^£^j-_sj)  .  (20) 


l>j=l  i>j=l 

"^4  [Eq.  (11)  ]  is  again  the  interaction  of  the  magnetic  moments  of  all 
ion  electrons  with  the  external  DC  magnetic  field,  and  it  is  evident  that  if 
the  DC  magnetic  field  is  not  strong  enough  to  decouple  the  angular  momenta 
of  the  electrons  then  the  result  of  interaction  on  the  closed  shells  is  zero. 
For  the  last  shell  this  result  is  as  shown  in  Eq.  ( 11) . 

^^5  [Eq.  (12)  ]  is  exactly  like"^^,  the  only  difference  being  that  the 
total  number  of  Z  protons  and  A-Z  neutrons  is  assumed  to  be  a  particle  with 
angular  momentum  I  at  the  origin;  therefore  the  contribution  of  the  orbital 
angular  momentum  of  the  electron  is  Just 


(21) 

where  the  contribution  of  the  ith  electron  spin  is 

=  -2gj^Pjj3  -  5r^-siri-I  -  8jr6(£^)  ks^- l/j ]|  (22) 

The  first  pant  of  Eq.  (22)  is  due  to  dipole-dipole  interaction  and  the  second 
to  the  electrons  present  at  origin.  Therefore  Eq.  (23)  results  (see  Fig.  9): 
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Fig.  9.  Schematic  representation  of  electron- 
nucleus  magnetic  interaction. 


[Eq.  (13)  ]  is  the  electric  quadrupole  moment.  Consider  that  the  nth 
nucleon  potential  on  the  ith  electron  is  ( see  Fig.  10) : 


and 


where 


Eind 


-  I 


£=0 


% 


A  00 


''NE 


I  I  =  I  I  ^Pni^i(^nl) 


i=l  n=l 


^  ^®n  r_o 


i=l  n=l  t=0 


—  fPni^oClnl)  +  Pni^i(lni)  +  Pni^aC ini) ^ 

i,n 

=  Vo  +  Vi  +  V2 

ini  “  ®ni  “  (^i^n''’yi^n‘''^i^)/^i^n 

P^i  =  (RnAi)^ 


(24) 


(25) 


(26) 
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I 


Fig.  10.  Schematic  representation  of  electron 
nucleon  coulomb  interaction. 


(27) 


(28) 


The  matrix  elements  of  Vi  correspond  to  the  nuclear  electric  dipole 
moment  which  vanishes,  and 

Va  =  Z  ^I^eie^p^iPad^i)  =  ^  (29) 

i,n  i,n 

Substituting  Eqs.  (25)  and  (26)  we  have; 
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^2  =  i  X  ^i'^^^eien[3rJ®R-^(xfxg+yfY^+zfzg+2xiyiXnYn+2yiZiYnZn+2ziXiZnXn)-l] 
i,n 

=  I  ^  ri^eien[xf5^+yfY2+zfZg+2xiyiXnYn+2yiZiYnZn+2ziXiZ^X^ 

i,n 

-  (xf+yf+zf)(xg+Y^^)/3] 

=  ^  ^  rJ^ej^ej^[9(x^)^+y^Y2+z^^)  +  l8(x;j^yj^XjjYn+  '••)  ■  5riRn] 
i,n 

=  1  I  rl'®i«nUrf-3xf)(R^-3)^)  +  (rf-3yf ) (R^-3Y^)  +  (rf-3zf) (Rg-3Z§) 

i,n 

+  l8(xiyiXnYn+  •••) 1 


5  I  '"'1  I 


Ur 


+  l8(XiyiXj,Yn+ 


Consider  now^ 


f-3xf)(  1^-33^)  +  (rf-3yf)(I^-3Yg)  +  (rf-3zf) (R|-3Zg) 
...)]  (30) 


I  enU^-5^) 


eS  - 
1(21-1) 


(I*^-3I^) 


(31) 


®n^^n 

n 


eQ  1 
1(21-1)  ^ 


( ^x^y'*’^x^y) 


(32) 


Thus  we  obtain: 


I*®( 3rf-3rf ) -3rf ( I^+I^+I^) +9( I^xf  • • • ) +l8( Xiyilxly+- • * 

(33) 

For  all  electrons  the  ei's  are  the  same,  ei  =  -e;  therefore  Eq.  (33)  reduces 
to: 


Va 


i  y  £i  . 

6  Zj  r? 


-®s. 


f  1(21-1) 
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-M' 


'^6  =  V2  =  V  =  y  r-5[l(l+l)r|  .  3(ri.l)^]  m 

^  21(21-1)  Y  - 

^^7  [Eq.  (l4)  ]  describes  the  effect  of  the  sum  of  the  coulomb  potentials 
produced  by  the  neighboring  ions  at  the  ith  electron.  In  the  case  of  approx¬ 
imation  of  the  closest  neighbors  in  ZnTe  and  CdTe,  this  is  the  sum  of  the 
potentials  produced  by  the  four  closest  Te  ions.  is  discussed  in  Section 

IV-C. 

^^8  [Eq.  (15)  ]  is  the  direct  interaction  of  the  DC  magnetic  field  H  and 
nuclear  spin  I  with  the  gyromagnetic  ration  g^^.  guPijI'H  is  considered  as  a 
constant  and  neglected. 


.  A>t>  t  J. 


B.  SPIN  HAMILTONIM 

In  case  of  II-VI  compounds  and  for  3d^  ions,  J  >  /  Is  +0^3'.  there¬ 
fore  the  ground-state  electronic  configuration  of  Mn'^'*'  is  still  3d5.  The  terms 
associated  with  this  configuration  are: 


S,P,D,F,G,H,  and  I.  (35) 

One  way  of  studying  the  effect  of  different  Hamiltonians  on  these  terms 
is  to  transform  them  into  functions  of  angular  momentum  operators,  L,  S,  and 
I,  which  facilitates  the  theoretical  calculations  considerably  and  does  not 
affect  the  ESR  experimental  calculations.  This  transformation  is  carried  out 
in  Appendix  B.  The  result  [Eq.  (B-14)  of  Appendix  B]  is; 
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Kl-S  -  M[(L"S)^  +  I  i'S  ”  L*^S*^/3] 


+  p(  L+2S)  •  H  -  1 


+  N 


-  3(L*^-I+L-^-S)/2]  +  L-I  -  kS*i| 


+  p[(L*l)2  -  L-l/2 


(36) 


vhere 


L*  =  n/l(L+1)  j  S*  =  n/S(  S+1)  ;  I*  =:  s/ 1(1+1) 

and  K,  M,  N,  p,  and  k  are  independent  of  L,  S,  I,  Mj^,  Mg  or  Mj,  The  energy 
of  each  of  the  terms  of  Eq,.  (35)  is  calculated  in  Appendix  A.  Equations 
(A-28)  of  this  appendix  give  for  Mn"*"*"  (3<i^) 


E(  3d^ ,  ,  B/C=.  3 .  C)  -lOA 

C 


0.5 

for 

0 

a 

for 

-  T.o 

for 

4p 

-  7. 50 

for 

\^-17.30 

for 

(37) 


where  A,B,C  are  Racah's  coefficients. 

Thus  the  ®S  state  lies  the  lowest;  considering  Eqs.  (36)  and  (37)  one 
finds  that  in  the  case  of  Mn'^'*"  the  eigenvalues  for  operators  including  L  will 
he  zerc ,,  Thus : 


<  I^L 


IMl> 


2&S- H  +  S»I 


(38) 
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where 


A  =  -Nk 

Therefore,  for  Mn"^  in  the  intermediate  crystalline  field  appropriate  to  II- 
VI  compounds  such  as  MgO,  ZnS,  ZnTe,  and  CdTe,  the  spin  Hamiltonian — which  is 
derived  from  the  total  Hamiltonian  after  the  latter  has  been  transformed  into 
functions  of  orbital  and  spin  angular  momenta  and  then  integrated  over  orbital 
angular  momenta  L — is : 

(5.1)  +  =  2PS-H  +  AS-I  +  (39) 

where 

A  =  -Nk 

and  %PiiH*I  is  assumed  to  remain  constant.  On  the  basis  of  the  assvimptions 
made  to  derive  ^ ,  In  Eq.s.  (7)"(15)>  one  findsthat: 

(1)  The  gyromagnetlc  ratio  of  Mn’*"*'  outer  shell  electronic  cloud  should 
be  the  same  as  that  of  free  electrons. 

(2)  Any  difference  between  electronic  spin  states  Mi( -5/2, • • • 5/2)  should 
result  from  ^7(S)  and  AS-I  in  Eq.  (39) . 

C.  CRYSTALLINE  FIELD 

The  crystalline  field  in  II-VI  compounds  of  zincblende  structure  can  be 

obtained  from  Fig.  11  and  Eq.  (4o). 

00  -<1 

''('■)  •  I  I  (*") 

n=0  m=-n 

Symmetry  considerations  will  reduce  the  number  of  terms  in  Eq.  (4o). 
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Fig.  11.  Schematic  representation  of  crystalline  field. 

(A  is  the  ion,  B  is  its  closest  neighbor,  and  C  is  its 
electron. ) 

For  d  electrons  one  has  to  consider  only  expansions  up  to  n  =  t.  No 
terms  with  odd  n  in  a  given  configuration  will  contribute  in  ground-state 
splitting,  regardless  of  whether  the  center  of  symmetry  is  present  or  not. 

The  term  A^  is  a  constant  shifting  all  of  the  levels  of  a  given  configuration 
by  the  same  amount.  Therefore: 

=  Aao^^Yao  A44Y44  +  A4_4Y4_4 


The  choice  of  2  axis  parallel  to  the  fourfold  rotary  reflection  axis 
S4  facilitates  the  determination  of  the  coefficients  in  the  above  relations. 
Furthermore  consider  that 

V(r,9,«  -  X 

nm 

where 

Onm  =  f(r,n,m)  only 

and 

i  =  cos  9 
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Thus 


results  in 


or 


S4V(r,e,^)  =  V(r0+rt,  ,!$  +  |)  =  ^r,Q,^) 


[cos(0+ir)  ]'^ 


im(  ^+n/2) 
e 


n-m 

(cos  0) 


m  =  0,  ±U  . 


Now  consider  that  reality  condition  requires  that 


V 


* 


V 


or  each  term  having  the  same  condition,  i.e,, 

(A44Y44)  =  A44Y44  =  A4-4Y4-4 

since 


y:4  =  (-1)^Y4.4  =  Y4.4  =IDA44  =  A4.4  (4la) 

The  presence  of  reflection  planes  o  requires  that 

“  ^-m  (^1^) 

Conditions  (Ula)  and  (4lh)  result  in: 

~  '\i-m  “  ^^nm  > 

namely,  Aj^j^'s  are  real.  In  our  experiments  no  axieJ.  dependence  of  the  field 
on  0  was  observed;  therefore  the  field  set  up  by  four  tellurium  ions  around 
Mn'*’’^  electronic  cloud  is 
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V  -  A40  ^^^^40  (^4/^40)  (^44"*'Y4-4)  ] 

=  A4or‘^(  35  cos‘^0-30  cos^©+3)  +  r'*A44  n/to  sin'*©  cos'^^ 

Now  consider  Fig.  12.  We  have 

V^Cr,©  =  0,^  =  0)  =  VB(r,©  =  n/2,^  =  O) 

or 

8r^A4o  =  3A4oI-^  +  sItO  A44 

01 

( A44/ A40)  ~  'J^/lh 


Fig.  12.  A  unit  cell  of  cubic  ZnTe  single  crystal  showing  a 
Mn"*"*"  ion  surrounded  by  four  Te  ions  and  the  reflection  planes 

Therefore 


1/2 

7(v,9,^)  =  V  =  A4or‘*[Y4o+(5/l'+)  (Y44+y4.4)  ] 


and  Eqs.  (B-28)  of  Appendix  B  give: 


In  the  ®S  state  of  Mn"’"*',  for  H  parallel  to  z  axis,  one  finds  from  Eqs.  (39) 
and  (44)  that 

<  M,n|^^|M,m  >  -  <  M-l,m|^  |M-l,m  >  =  2PH(M,m)  +  mA  =  hvo  =  2PHq 

(45) 

and 


H(M,m)  =  Ho  -  mA*  A*  =  A/2P  (46) 

where 

m  =  nuclear  magnetic  quantum  number. 

Equations  (45)  and  (46)  indicate  that  there  should  be  six  absorptions  only 
( see  Fig.  13) . 


-  5/2  |A|  -3/2  |a|  -1/2  |a|  1/2  |a|  3/2  |a|  5/2  |A| 


i 

1 

Ho 

H 


Fig.  13.  Theoretical  prediction  of  Mn’*"*'  EBR  spectrum. 


These  absorptions  are  observed  only  in  powders  and  we  realize  that  to  de 
scribe  the  spectra  shown  in  Figs.  3-8  a  further  refinement  of  theory  is  neces 
sary.  Bleaney  and  Stevens,^®  after  a  phenomological  argument,  conclude  that 
a  more  reedistic  representation  of  the  crystal  potential  for  Mn"*^  in  II-VI 


compounds  may  be  written  as: 


1^7[(r,r')^,L,S]  =  B[I^  +  +  l4  -  L*^(3L*^.l)/5] 

+  P[S4  +  S4  +  Sf  -  S*^(3S*^-l)/5] 


(47) 


Therefore 


-  <o|j(fjo>+^  l<0|K'ln  >|^/[E(0)-E(n)  ] 


=  g3H-S  +  AS- I  +  P[S$  +  +  S|  -  S*  (3S*  -l)/5]  (48) 


Owing  to  the  fact  that  P  is  a  measure  o: 


-  'Jt 


cubic 


when  H  =  0  it  is  called  the 


zero  field  splitting  factor.  Bethe^^  shows  that  the  state  ®S  ,  will  be  split 

5/2 

into  and  ‘^Fg  under  the  effect  of  a  tetrahedral  field  created  by  the  four 
closest  tellurium  ions  of  Zn,  Cd,  and  Te.  To  correlate  P  and  AE  =  “^Fg-^F,, 
use  is  made  of  T  tensors  given  as;^^^^5 


-i/a 

’’fqil  " 


(49) 


S+  =  Sx  ±  iSy 


-4+4 


'Pjo/ilG)  S± 


(50) 


Sx  =  (S4  +  S.)/2 


Sy  =  (S+  +  S.)/2i 


(51) 


Consider  that 


-  (3^.  +  S.)  /l6  =  (S4  +  Sl)/l6  +  f{s3  +  --Sf) 


=  (S+  +  S.)Vl6  =  (S^  +  Sl)/l6  +  h(S?  +  . . -Sf) 
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Then 


P(S^  +  sj) 


=  ^{St  +  Si)  +  f  4-  h-  = 


2P 


^  (T44  +  T4.4)  +  f"  +  h" 

VTO  (52) 


Evaluating  S+  from  Eq.  (U9)  and  substituting  in  Eq.  (52)  one  gets: 


(Mn++) 

cubic 


gpS-H  +  M.S  +  [T40  +  n/5M  (T44  +  T4-4)  (53) 


where 


T40  =  I  (35S|  -  30S*^S|  +  25Si  -  6s*^  +  3S*^)  S*  =■  ^S(S+l)  (54) 

The  matrix  elements  of  the  2P[T4o+  ^ 5/l4  (T44+T4_4) ]/5  can  be  found  very 
easily,  and  the  result  is  as  shown  in  Table  II. 


TABLE  II 

MATRIX  ELEMENTS  OF  |£  [T40+  \/5/l4  (T444T4_4) ] 


M'  = 

1 

+  — 

2 

~  2 

4  5 

2 

M  =  ±  - 
2 

6p 

0 

0 

+  5 

2 

0 

-9P 

3  s^P 

45 
~  2 

0 

3  n/5P 

3P 

In  the  secular  equation 

113^1  j-6ije II  =  (6P-€)[(-9P-e)(3P-e)  -  =  0 
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giving 


ei(±  |)  =  6P,  €2  =  6p,  €3  =  -12P  (55) 

€i  and  e2  correspond  to  fo\ir  eigen  factions  with  energy  6P  above  the  ground 
state  of  the  free  ion.  Thus  the  situation  shown  in  Fig.  14  results: 


1 - r-r 

a  6p 


^Sp;/g(Mn++  free  Ion)  / _ f 

\  ~  T  ” 

\  2a  12P 

_ L^r 


Fig.  l4.  Effects  of  cubic  field  on  levels. 


The  splitting  of  state  in  a  tetrahedral  field  is  given  by 


'‘Fg  -  %  =  18p  =  5a 


where 


a  =  6p 


(56) 


Substituting  for  Pin  Eqs.  (52)  and  (53),  one  finds  expressions  similar  to 
those  given  by  Bleaney  and  Stevens^^  and  by  Lambe  et  al.^: 


'75 


=  gPS-H  +  a[s^+S^+S4-S*®(5S*^-l)/5]/6 


cubic 


cubic 


+  AI-S  - 


(57) 


or 


^  =  gpS.H  +  AI.S  +  a[T4o+N/5/l4  (T44+T4.4)  ]/l5  (58) 

^  ^  cubic  ”■ 


To  find  the  energy  corresponding  to  electronic  and  nuclear  spin  levels 

S 

f . s,  " 


we  set  up  the  secular  equation  for  fine-structure  Hamiltonian^ 


XL  ■  X  -  M-S  =  gPS-H  +  a(35S|  -  30S*^S| 


+  25S|  -  6S*^  +  3S*^)/120  +  a(S+  +  S^)/U8 


(59) 


For  the  DC  magnetic  field  H  parallel  to  one  of  the  principal  axes  of  the  cube, 
namely  the  z  axis,  we  will  find  the  matrix  elements  of 


(HIIZIISJ 

f*  •  S « 


gPHS2  +  a(35S|+25S|-30S*^S|-6s*^+3S*^)/l20  +  a(S++sl)/U8 

(60) 


The  first  and  second  term  are  diagonal  while  the  third  term  is  off-diagonal; 
therefore 

<  M|gPHS2  +  a(35S|-30S*^Sf +2531-63*^+33*“*^) /120  +  a(  3++S.)  A81m'  > 

=  <  MigPHSz  +  a(353^-303*^S|+253|-63*A33*^)/l20|M'  > 

t  <Mla(sj4Sl)/1.8|M'  >  -  (61) 

with 

A  =  [2Me  +  a(lUM^ -951^+189) /48] 

B  =  V5  a/2  6  =  gPH/2 
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The  secular  equations  (determinant)  can  be  constructed  from 


MM' 


0 


(62) 


Thus 


\  m' 

M  \ 

5/2 

3/2 

1/2 

-1/2 

-3/2 

-5/2 

5/2 

5€+a/2-E 

0 

0 

0 

>/5  a/2 

0 

3/2 

0 

3e-3a/2-E 

0 

0 

0 

n/5  a/2 

1/2 

0 

0 

e+a-E 

0 

0 

0 

-1/2 

0 

0 

0 

-e+a-E 

0 

0 

-3/2 

V5  a/2 

0 

0 

0 

-3e-3a/2-E 

0 

-5/2 

0 

V?  a/2 

0 

0 

0 

-5€+a/2-E 

(63) 


The  determinant  is  an  equation  of  6th  power  in  E.  To  solve  consider 


A  = 


a  0  0  0  7  0 
0  p  0  0  0  7 
0  0  6  0  0  0 


0  0 
7  0 

0  7 


\  0  0 

0^0 
0  0  Tl 


=  6\(a^-72)  (611-72)  =  0 


(6M 


giving 


El  =  a+e 


Eg  =  a-€ 


or,  provided  e  »  a, 
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Eg  =  e  -  ^  +  ^{he+B.)^+5a.^/k  ^  e  .  |  +  Ug  +  a  +  ^/Se 


M  - 


1L5  - 


^  ■  2  ■ 


Ee  — 


2€+a 

2 

2e+a 


s/(Ue+a)^+5a^/U  ^  e-^-4e-a  -  ■^—/Se 


+  *y  ( 4€-a)  ^+5a^/4  =  -e-^+4e-a  +  ^^/Se 


(65) 


N/(4€-a)^+5a^/4  =  -e-|^-4€  +  a-  ■^^/Se 


It  is  obvious  that  the  states  M  =  ±l/2  do  not  mix  and  that  E^  and  E2  there¬ 
fore  belong  to  M  =  l/2  and  M  =  -l/2.  For  other  m' s  one  has,  provided  a/e  «  1. 


E^*^‘±  i  =  a  ±  € 
2 


rX.S. 


(M  =  ±  i)  =  -  ±  3e  ±  ^/8e  =  -  ±  Je  ±  5a732€ 


(66) 


.f.s. 


(M=±^)  =|±5e±  ^/Se  =  |  ±  5€  ±  5aV526 


Thus,  the  energy  diagram  of  spin  levels  for  H|1z  and  in  the  absence  of 
nuclear  spin  effect  is  as  shown  in  Fig.  15  and  Table  III. 


TABLE  III 

VARIATION  OF  SEPARATION  OF  Mn’^  ESR  FINE- STRUCTURE  COMPONETS 
AT  0  =  0  AS  A  FUNCTION  OF  p  =  gPH/2a  =  e/a 


p 

E-s/a-Ejr 

E-3/2-Ep 

1 

CM 

H 

• 

Ei^2-Ejp 

l\3 

1 

Es/a-Ep 

0 

-  2.00 

-  2.00 

1.00 

1.00 

1.00 

1.00 

1 

-  4.71 

-  4.62 

0.00 

2.00 

1.71 

5.62 

2 

-  9.59 

-  7.57 

-  1 

3 

4.59 

10.57 

3 

-  14.56 

-  10.55 

-  2 

4 

7.05 

15.55 

4 

-  19.54 

-  13.54 

-  3 

5 

10.54 

20.54 

5 

-  24.53 

-  16.53 

-  4 

6 

13.53 

25.54 

6 

-  29.52 

-  19.53 

-  5 

7 

16.53 

30.52 

7 

-  54.52 

-  22.52 

-  6 

8 

19.52 

35.52 

10 

-  49.52 

-  31.51 

-  9 

11 

28.51 

50.5 

100 

-499-5 

-301.50 

-99 

101 

298.5 

500.5 
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For  the  case  where  H  has  direction  cosines  a,  p,  and  7  with  coordinate 


axes, 

gpH.S  =  gp[HzSz  +  I  {I4S.+H.S+)  ]  =  2e(7S2+VS_+\_S+)  (6?) 

where 

2e  =  gPH  and  \+  =  (a  ±  iP)/2  (68) 

and  the  energy  levels  are^^ 

E±i/2  =  pa  ±  €  +  ria^/e^  ±  q^a^/e 

E+3/2  =  -3pa/2  ±  3e  +  rga^/g^  ±  q^a^/e  (69) 

E±5/2  =  5pa/2  ±  5e  +  raa^/e^  ±  qga^/e 

where 

p  =  [1  -  (a2p24.p2^+^Q2)  ]  =  1  _  5 

=  -56(7-256) /6 

(70) 

=  5 [1  +  6(22-756) ]/32 

qg  =  5[1  +  6( 50-1136) /31/52 

and 

ri  =  -56( 196-16356+31256^) /lh4 

ra  =  5 [1  +  36(79-6156+11256^) ]/l28  (71) 

rg  =  -5[1  -  56( 113-7056+10756®) /91/128 


D.  HYPERFINE  SPLITTING 

The  hyperfine  splitting  Hamiltonian  for  a  cubic  field  is: 


AS-I  =  ASzIz  +  A(S+I_+S_I+)/2 


(72) 
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The  exact  energy  to  he  assigned  to  each  state  M,m  requires  solution  of 
a  36  X  56  secular  equation.  An  approximate  method  is  to  use  time-independent 
perturbation  theory  up  to  a  sufficiently  high  order  with  respect  to  the 
diagonalized  Hamiltonian  for  h1|z||S4. 

Let 


t  - 

1/0  •  ePHSj 


Z  - 


AS^Iz  +  A(S+I.+S.I+)/2 


crystalline 


Then,  taking  free  ion  energy  as  zero,  one  has : 


M,m  M 


(0)  (1)  (2) 

E  +  E  +  E 
",m  M,r 


,m  M,m 


where 


e(°^  =  <  M,m|  ]/^iM,m  > 

M,m  ^  » ''O  '  ' 

4^1  =  <M,m|JfjM,m> 


4!i  =  I 


m'  ,m' 


M.“iXix'’"'>i^/[4?i  •  4>' 


(75) 

m 

(75) 


(76) 


(77) 

(78) 

(79) 


Then  the  energy  associated  with  each  level  M,m  is  given  by 


where 


^f.S. 

M,m 


w  p(  1) 

-  \ 


m 


+  E 


;(2) 

‘M,m 


(80) 


ho 


(81) 


The  values  of  E^’^‘  are  given  by  Eq.  (66)  and  and  eJ^^^  are 


:(^)  ^(2) 


4^1  - 


I 

=  Z  I  <M,m|A(S^I.+S.lJ/2|M',m'>lt[4°i  -  4'lm’] 


M'  ,m' 


=  (<  M,m|s+I.|M' ,m'><M^S_I^|M' ,iii'>) 

(x)  <  M'  ,m'  |S+I.  +  S.I+|M,m  ^ 

=  |-  <  M,m|S+I.|M-l,m+l  >  <  M-l,m+l lS_I+ lM,in  >/[4°ii  “  4-1, “+1^ 


+  ^  <  M,m|S.I+lM+l,m-l  >  <  M+l,Tn-l  |S+I.  lM,m  >/[4°ili  ’  4+1, m-l^ 


Consider 


S+I.  lM-l,m+l>  =  [S*^  -  (M-l)M]^/^[l*^  -  (in+l)in]^/^|M,in  > 

S.I+lM,in  >  =  [S*^  -  M(M-1)  -  m(m+l)  ]V2lM-l,m+l  > 

S.I+lM+l,m-l  >  =  [S*^  -  (M+1)M]^/®[I*^  -  (in-l)m]^/^|M,m  > 

S+I_|M,ni  >  =  [S*^  -  M(M+1)  -  m(m-l)  ]^/^|M+l,m-l  > 


Consider  also 


Jo) 

M,m 


gPHM, 


Jo)  rg 

Mtl,mTl 


gPH(M±l), 


e(o) 

M,m 


M+l,m±l 


±g6H  (8U) 


Substituting  Eqs.  (82)  and  (85)  in  Eq.  (8i+),  one  finds: 


gPH  gPH 

■^—  [M(l*^-m2)  -  in(S*^-]yP)  +  M(l*^-m2)  -  m(S*^-M2)  ]  (85) 

4gPH 

-^'|m[I(I+1)  -  -  m[s(S+l)  -  M^]  ' 

2gpH  L  J 


Considering  Eqs.  (67),  (75)>  (80),  (81),  and  (85)  one  finds  for  (H||z||S4)  : 


m 


^3/2, m 


^-3/2, m 


E_ 


5/2 


e  -  I  +  sl(ke+a)  ^+5a.^lk  +  5Ain/2  +  A^[5(35A-ni2)/2-5m/2]/2gPH 
-€  -  I  +  '«/(4€-a)2+5a2/U  +  3Am/2  +  A2[3(35A-m2)/2-13ra/2]/2gPH 
e  +  a  +  Ara/2  +  A^[(35/4-m^) /2-17m/2 ]/2gPH 


(86) 


-€  +  a  -  Ara/2  +  A^[-(35/4-m2)/2-17m/2]/2gPH 
^  ~  2  ~  •  3Am/2  +  A^[-3(35/4-m^/2-13m/2  ]/2gpH 

-€  -  I  -  V(4e-a)2+5a2/4  -  5Ara/2  +  A^[-5(55/^-ni^)/2-5m/2]/2gPH 


For  a  given  m  there  are  six  energy  levels;  therefore  five  absorption  transi¬ 
tions  will  occur,  giving  rise  to  five  resonance  lines  corresponding  to  a 
given  ra  and  to  different  M  values.  The  ESR  spectrum  of  Mn'*"*'  clearly  shows 
this  feature  (see  Fig.  l6) . 

Considering  Fig.  l6  one  notes  that  associated  with  each  ra  the  spectrum 
at  9  -=  C  has  the  form  shown  in  Fig.  I7.  To  identify  the  M  values  to  which 
the  lines  a,  P,  7,  6,  and  \  are  associated  we  must  find  the  intensity  of 
each  line.  This  is  proportional  to  magnetic  dipole  absorption  probability; 
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Fig,  l6.  ESR  spectrum  in  ZnS  (cubic)  at  3CX)‘’K,  0  =  0°. 


a  p 


6  \ 


Fig.  17.  Mn’*"*'  ESR  due  to  a  given  m. 


and  the  ratio  of  intensities  of  two  lines  is  as  follows 


Im 


<  M-1  rK-r  M  > 


(87) 


M'  |<M'-l|rK.r|M'>|^ 

But,  from  replacement  theory,  vectors  of  position  and  linear  mcanentum  can  he 
substituted  by  angular  momentum  operators  L  or  S,  namely 


(88) 


where 


^  1  T,  ^ 

r  ~  P.r  oc 


•li 


AA  A 
LL.L 


A 

or  SS 


*2 


A  A 

S  =  Sz  +  |_S+  +  l+S.  u 


(l±l)/2  i  =  -PL 


and 


Ij^  <  M-lll+S.  |M  >  <  mIi.S+Im-I  >  S(S+1)  -  M(M-l) 

<  m' -i|i+s. |m'>  <  m' U_s+|m'-i>  s(s+i)  -  m'(m'-i) 


Thus 


=  =  5:8;9:8:5  (89) 

Note  also  that  the  separation  of  the  measured  magnetic  field  correspond¬ 
ing  to  each  line  with  intensity  ,  p  for  0  =  0  can  be  obtained  as 

follows : 


^M,m  -  ^-l,m  =  eP«M,m  ^  ^M,m  =  "  ^Vo  for  all  M,m  (90) 


Thus 


%,m"%-l,m  %',ra”%'-l4n  _ 


(91) 


or 


*^MM'  ,m  ^,m  ^M,m  ,m 


f  =  f7g3 


(92) 


where  f'  (A,a)'s  are  given  in  Eq.  (86)  and  one  considers  the  fact  that  in 
M,m 

compounds  ja^^  <  |A^^  ,  f^^  ^  must  be  found  up  to  a^,  i,e., 


Mn 


Mn 


1+1+ 


^/(4€±a)  2+5a2/U  ^  1+6  ±  a  +  5a2/l6gPH  gPH  =  26 


Substituting  Eq.  (93)  into  Eq.  (86)  one  finds  a  ,'s,  e.g., 

MM 


where 


'i/2,5/2  ■  ^1/2  ^l/2  ■  ^5/2 


V  =  -  V 


Equations  (86)  give 


\/2  -  ^-1/2  ^  2€+inA+A2(33/^-m2)/2gpH  .  /gp  =  h 

-Q  1 


gP 


gP 


Thus 


^1/2  =  mA'  +  A'  [(35/i+-in2)  ]/2H  A'  =  A/gp,  a’  =  a/gP 


and 

=  hiA'  +  2a'  +  A'^(35/^-in2+8in) 

Therefore 

=  -2a'  -  2mA' ^/h 

Similarly 

S/2  =  5a’ /2  -  5a'2/l6H  -  mA'^/H 

S/2  =  0 

=  -5a' /2  -  5a'2/i6H  +  mA'^H 

Sa/a  =  2a'  +  2mA' ^/H 

Considering  Eqs.  (89)  and  (96),  and  also  the  fact  that  the  two  fine 


(93) 


(94) 


(95) 


(96) 

structure 


components  tend  to  merge  (see  Figs.  I7  and  I8)  at  higher  field  one  can  deter- 
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Fig.  l8.  Schematic  plot  of  fine  structure  of  ESR  spectrum  of 
Mn++  in  cubic  ZnS  for  m  =  -5/2  and  5/2. 

mine  that  A  and  a  have  opposite  signs.  According  to  Watanabe5  a  oc  (Dq)^  and 
thus  A  <  0.  Since  ( |5a'/2 |)/( |5A'^/h|)  is  8,  28,  and  29  for  ZnS,  CdTe,  and 
ZnTe  respectively,  the  measured  DC  magnetic  field  determines  the  M  to  be  as¬ 
signed  to  the  lines  a,  P,  y,  6,  e.g. ,  a  corresponds  to  the  M  giving  max¬ 
imum  separation,  etc.  The  result  (see  Fig.  18)  is 

a  M=  -1/2;  p  5/2;  y  l/2;  6  -5/2;  \  5/2  (97) 

Consider  Eqs .  ( 9I)  and  ( 95) ,  as  well  as  the  fact  that  A*  <  0  and 


M,m 


-  H 


M,m 


=  f 


M,  m  M,  m 


Then 


^1/2, m'  ^1/2, m 


r 


A'(m-m')  +  A'^(m'2-m2)/2Ho  ] 


>  0 


<  0 


m'  >  m 
m'  <  m 


(98) 


Therefore,  the  lowest  group  of  five  fine-structure  components  a,  P,  7,  6, 
and  \  (see  Fig.  18)  belongs  to  m  =  -5/2,  where  the  group  at  the  high-side 


1*6 


field  belongs  to  m  =  +5/2. 


E.  DETERMINATION  OF  3a  =  PARAMETERS 

The  electronic  level  of  the  Mh"'”''  ion  is  split  into  two  levels,  ^Fg 

p 

and  F^,  due  to  the  presence  of  a  cubic  field.  This  energy  difference  is 
easily  obtained  from  Eq.  (96)  : 

I  [ |a(M  =  -1/2, m  =  -5/2)  I  +  |a(M  =  3/2,m  =  5/2) |  ] 

=  5  [|-5a'/2  -  5a'2/l6H  -  5A'V2H|  +  |5a'/2  -  5a'^/l6H  -  5A' V2H|  ] 

=  3[|5a'/2  +  5a'Vl6H  +  5A'^/2H|  +  |5a’/2  -  5a'^/l6H  -  5A'^/2H|  ]/5  =  3a' 
and 

3a  =  gP  3a'(erg) ,  3a'(gauss),  or  (cm-^)  (99) 


F.  DETERMINATION  OF  HYPERFINE  SPLI'TTING  CONSTANT  A' 
Equation  (98)  gives 


H, 


L/2^5/2  “  ®i/2-5/2  ^  ^ 


(100) 


A  =  gpA'(erg)  or  a'(  gauss)  or  ggA'/lic  in  (cm"^) 


APPENDIX  A 


DETERMINATION  OF  Mn++  GROUND  STATE 

The  electronic  configuration  of  Mn"*"^  is 

Is^  2s^  2p®  3s^  3p®  (A-1) 

According  to  Griffith-^®  the  relative  energy  levels  of  the  ion  can  be  obtained 
by  considering  only  the  last  five  electrons  in  the  partially  filled  shell  of 
3d.  The  energy  terms  corresponding  to  these  electrons  can  be  found  simply  by 
the  "counting  method."  This  method  consists  of  setting  up  a  table  with  pos¬ 
sible  positive  Mg  and  Mj^  values  in  rows  and  columns^  respectively,  with  deter¬ 
mining  the  number  of  times  that  each  Mg  occurs  for  a  given  M^^  by  constructing 
all  possible  wavefunctions  corresponding  to  M^.  In  Eq.  (A-2)  a  set  of  wave 
functions  of  M^  =  2  is  given. 


|2+  2-  1+  -1+  -2'b> 

with 

Ml  =  2 

and 

Mg 

=  5/2 

|2+  2"  1+  -1+  -2"> 

with 

OJ 

II 

and 

% 

=  1/2 

|2+  2“  1+  -1-  -2‘b> 

with 

Ml  =  2 

and 

Mg 

=  1/2 

|2+  2-  1-  -1+  -2+> 

with 

Ml  =  2 

and 

Mg 

=  1/2 

One  can,  with  the  same  procedure,  find  that  there  are  nine  other  wave 
functions  with  Mg  =  l/2  and  two  more  with  Mg  =  3/2.  The  result  for  other  M^ 
values  can  be  obtained  much  easier;  they  are  given  in  Table  A-I.  Among  the 
sixteen  terms  in  Table  A-I  the  sextet  ®S  ajid  quartets  ^G,  *F,  ^D  and  are 

the  most  responsible  terms  which  account  for  the  observed  optical  and  EPR 


TABLE  A- I 


MgM^  TABLE  FOR  3^^  CONFIGURATION 


\ms 

Ml\ 

Total  Nvimber 

Difference 

Energy  Terms 

1/2 

3/2 

5/2 

1/2 

3/2 

5/2 

6 

1 

0 

0 

1 

0 

0 

^I 

5 

2 

0 

0 

1 

0 

0 

4 

5 

1 

0 

3 

1 

0 

^G,^G  and  ^G 

3 

8 

2 

0 

3 

1 

0 

^F,^F  and  "^F 

2 

12 

3 

0 

k 

1 

0 

^D,^D,^D  and  “^D 

1 

14 

h 

0 

2 

1 

0 

^P  and  ^P 

0 

l6 

5 

1 

2 

1 

1 

and  ®S 

1  P 

experimental  results.  According  to  Racah-^'^  these  terms  can  be  calculated 
from  the  terms  of  because  Jd^  is  an  configuration  and  because  the 

sextet  and  quartet  terras  have  the  spin  S  =  /+l/2  and  S  =  i-l/2.  The  5d^ 

MgMj^  table  (Table  A-II)  reveals  that  for  Jd^  there  are  five  energy  terms, 

^G,  ^F,  ^D,  and  ,  (A-5) 


whose  energies  can  be  obtained  by  using  the  diagonal-sum  mle.  In  this  way 
one  finds  the  sum  of  the  energies  of  several  terms  from  the  relation: 


(A-4) 


where  is  a  polynomial  of  m^^th  power  provided 

P 


a=l 
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TABLE  A- II 


MqMl  table  for  3d^  CONFIGURATION 


\  Ms 
Ml\ 

Total  Number 

1  0 

Difference 

1  0 

Energy  Terms 

4 

0 

1 

0 

1 

3 

1 

2 

1 

1 

^F 

2 

1 

3 

0 

1 

^D 

1 

2 

h 

1 

1 

3p 

0 

2 

5 

0 

1 

with  n' 

the  n-umber  of 

electrons 

in  the  partially 

filled 

shell  (n'  <  2£+l) . 

Each  P(j  satisfies  the  relation 


P 


a 


=  0 


(A-5) 


i=l 

The  coefficient  of  -E  in  Eq.  (A-5)  is  the  desired  diagonal-sum  which  we 
must  find.  It  can  be  obtained  from  the  relation 


m. 


(A-6) 


i=l 

where  J(k,\)  is  the  coulomb  integral  defined  as 


J(i,j)  =  Yj  j) 


(A-7) 


k=0 


F^  is  the  so-called  Slater-Condon  parameter  which  for  a  pair  of  electrons 
i  and  j  with  the  same  n  and  i  takes  the  form 
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00  00 


(1,J)  =  (A-8) 

0  0  O  U 


and  a^(i,j)  in  this  case  (same  n  and  i)  can  be  given  as: 


=  c^(mi,m^)c^(mJ,mj) 


(A.9) 


with 


c^(m,m') 


4it 


2}r  n 


L(2k+1). 


im 


(A-10) 


From  the  diagonal-sum  rule  which  is  described  in  Eqs.  (A-4),  (A-5),  and 
(A-6)  we  can  write  the  energies  of  the  terms  of  3d^  as  follows: 


E(^G)  -  J(2,2) 


( A- 11a) 


E(3F)+E(1g)  =  J(2,l)+J(l,2) 


(A-llb) 


E(  ID) +E(  ^F) +E(  ^G)  =  J(2,0)+J(l,l)+J(0,2) 


(A-llc) 


E(  ^P) +E(  ^D) +E( ^F) +E( ^G)  =  J(2, -l) +j( 1,0) +j(0,l) +j( -1,2)  (A-lld) 


and 


E(  ^S)  +E(  3p)  +E(  ^D)  +E(  ^F)  +E(  ^G)  =  J(  2 ,  -2)  +j(  1 ,  -1)  +  J(  0, 0)  +j(  -1 , l)  +j(  -2,2) 


(A-lle) 


Equations  (A-6) -(A- 10)  reveal  that  j(i,j)'s  of  Eqs.  (A-ll)  can  be  obtained 

from  the  much  simpler  relation 

i  I 

J(i,j)  =  Y.  =  Z  ^k^2k  =  aoFo  +  aaFa  +  a4F4  (A-12) 


k=0 


krO 
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In  Table  A-III  the  Integrals  j(i,j)  of  Eqs. 

of  the  coefficients  A,  B,  and  C  defined  as 

(A-ll)  are  given  in  terms 

A  =  Fq  -  U9F2 

(A- 13a) 

B  =  F2  -  5F4 

(A-13b) 

c  =  35F4 

(A-13c) 

TABLE  A-III 

INTEGRAIB  IN  TERMS  OF  A,  B,  MD  C 


m 

m' 

J(m,m' ) 

2 

±2 

A  +  1+B  +  2C 

2 

±1 

A  -  2B  +  C 

2 

0 

A  -  4B  +  C 

1 

±1 

A  +  B  +  2C 

1 

0 

A  +  2B  +  C 

0 

0 

A  +  4B  +  5C 

Substituting  for  in  Eqs.  (A-ll)  one  finds  the  relative  energy  levels 

as  follows : 


E(  iG)  =  A  +  4b  +  2C 

( A-l4a) 

E(3f)  =  A  -  8B 

(A-l4b) 

E(  ^D)  =  A  -  3B  +  2C 

( A-l4c) 
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(A-lUd) 


■1^  ^ 

E(  =  A  +  7B 

E(  ^S)  =  A  +  ikB  +  7C  ( A-lUe) 

Equations  (A-14)  reveal  that  the  relative  energy  levels  for  the  terms  with 

1  ^ 

S  =  n'/2  (n'  being  the  number  of  electrons)  can  be  express  as"^ 

E[d^' ]  =  n'(n'-l)(A-8B)/2  +  3[6n'  -L(L+1)]b/2  (A-I5) 

2  3 

and  therefore  the  energy  of  the  triplets  of  3d  >  quartets  of  3d  ,  quintets  of 
3d'^  and  sextets  of  3d^  can  be  found  easily,  e.g., 

E(dS,®S)  =  lOA  -  35B  (a- 16) 

The  energies  of  the  quartets  of  3d^  cannot  be  obtained  as  easily,  it  can, 
however,  be  shown  that  they  may  be  derived  from  Eqs.  (A-lU).  Consider  that 
each  quartet  of  d^  has  four  electrons  with  a  spin  and  one  with  p  spin.  We 
construct  the  wave  functions  ^i(m,m')  of  the  quartet  of  d^  so  that,  com¬ 
pared  with  the  wave  functions  \|fi(mi- • -m®)  of  the  quintet  of  d® 

\(fi(mi- •  •m6,d®,^D)  =  |2'^  2‘  l"^  O’*"  -l"^  -2'^>  ,  (A-I7) 

it  will  have  the  form 

(!5i(ra^.  .  .m®,d®,^F)  =  j!^i(2,l)  =  12'*'  2"  l"^  -2‘^>  .  (A-I8) 

In  Eq.  (A- 18),  m(2)  stands  for  the  m^  of  one-electron  wave  function  with  p 
spin  and  m'(l)  stands  for  the  -m^  of  the  extra  electron  in  the  quintet  of 
d®  configuration.  The  diagonal  elements  of  ^(m,m')  can  be  obtained  by  relating 
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In  Table  A- III  the  integrals  J(i,j)  of  Eqs.  (A-ll)  are  given  in  terms 


of  the  coefficients  A,  B,  and  C  defined  as 

A  =  Fq  -  I+9F2 
B  =  F2  -  5F4 

c  =  35F4 


TABLE  A-III 

J(m,m')  INTEGRALS  IN  TERMS  OF  A,  B,  AND  C 


m 

m' 

J(m,m') 

2 

±2 

A 

+  4b  +  2C 

2 

±1 

A 

-  2B  +  C 

2 

0 

A 

-  4B  +  C 

1 

±1 

A 

+  B  +  2C 

1 

0 

A 

+  2B  +  C 

0 

0 

A 

+  4b  +  3C 

Substituting  for  j(l,j)  in  Eqs.  (A-ll)  one  finds  the  relative  energy 
as  follows : 


E(  iG)  =  A  +  4B  +  2C 
E(3f)  =  A  -  8B 

E(^D)  =  A  -  5B  +  2C 


(A-lJa) 

(A-13b) 

(A-13c) 


level s 

(A-l4a) 

(A-l4b) 

(A-14c) 
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( A-l4d) 


E(^P)  =  A  +  7B 

E( ^S)  =  A  +  14B  +  7C  ( A-l4e) 

Equations  (A-14)  reveal  that  the  relative  energy  levels  for  the  terms  with 

1  6 

S  =  n'/2  (n'  being  the  number  of  electrons)  can  be  express  as"^ 

E[{jn' ^n'+l^(L)  ]  ^  n' (n' -l)  (  A-8b) /2  +  3  [6n'  -L(L+1)]B/2  (A-15) 

2  3 

and  therefore  the  energy  of  the  triplets  of  ,  quartets  of  ,  quintets  of 
3d‘^  and  sextets  of  3<i^  can  be  found  easily,  e.g., 

E(d=,®S)  =  lOA  -  35B  (A-16) 

The  energies  of  the  quartets  of  3<3^  cannot  be  obtained  as  easily,  It  can, 

however,  be  shown  that  they  may  be  derived  from  Eqs.  (A-14).  Consider  that 
each  quartet  of  d^  has  four  electrons  with  a  spin  and  one  with  p  spin.  We 
construct  the  wave  functions  )^i(m,m')  of  the  quartet  of  d^  so  that,  com¬ 
pared  with  the  wave  functions  \lfi(mi' • -m®)  of  the  quintet  of  d® 

ti(mi- •  .m6,d®,^D)  =  12'^  2'  l"^  O'^  -l’^  -2^>  ,  (A-17) 

It  will  have  the  form 

^i(mi.  ..m®,d®,'‘F)  =  ji^i(2,l)  =  |2'^  2'  l"^  -2+>  .  (A-18) 

In  Eq.  (A-18),  ra(2)  stands  for  the  m^  of  one-electron  wave  function  with  p 
spin  and  m’(l)  stands  for  the  -m^  of  the  extra  electron  In  the  quintet  of 
d®  configuration.  The  diagonal  elements  of  ^(m,m')  can  be  obtained  by  relating 

54 


them  to  the  corresponding  elements  of  the  terms  of  and  d°  configurations. 

These  terms  can  he  found  from  the  relation^^ 

6 


E(d®,®D)  =  ^  [J(m'^,m^)  -  K(m^,m^)  ] 


k<\=2 

5 


^  [j(m'^,m^)  -  K(m'^,m^)  ]  +  ^  [j(m'^,m®)  -  KCm'^^m®)  ] 


k<X=2 


ic=0 


(A- 19) 


where 


=  E[(!^(m,m')]  +  J(m,m')  -  K(m,m')  +  ^  f(m'^,m^) 

K=1 


C. 

K(m,m')  =  6(s,s')  ^  =  0  £  ^ 

k=0 


(A-20a) 


and 


f(m'^,ra^)  =  J(m'^,m^)  -  K(m'^,m^)  . 


(A-20h) 


Thus  Eq.  (A-I9)  can  he  rewritten  as 

5 

E(d®,^D)  =  E[^(m,m')]  +  j(m,m')  +  ^  f(m'^,-m') 


K=1 


We  have  also 


(A-21) 


55 


E(d^,®S) 


y 


k<X=2 


5  5 

=  ^  f(ni'^,m2)  +  ^  f(m'^,m^) 


K=1 

5 


ic<\=3 

1+ 


K=1 


K<x=2 


(A-22) 


J 


•m')  +  E(d'^,^D) 


K=1 


A  comparison  of  Eqs.  (A-21)  and  (A-22)  gives: 


E[j!$(m,ra')]  =  E(d®,®D)  +  E(d^,^D)  -  E(d®,®S)  -  j(m,m')  (A-23) 


Equation  (A-23)  gives  the  energy  levels  associated  vith  the  sextet  and  the 
quartets  of  Mn++  and  Ee’*^"''  within  a  constant  which  can  he  found  hy  consider¬ 
ing  Eqs.  (A-lh) ,  (A-16),  and  (A-23)  for  the  case  where  Ml  =  m+m'  =  0: 

E[jzS(m,m'),  Ml  =  0]  =  E(d®,®S)  =  C  -  E(d^/s)  (A-21+) 

where 


C  =  E(d®,®S)  +  E(d^,^S)  =  lOA  -  35B  +  A  +  IUB  +  JC 

=  llA  -  21B  +  7C  .  (A-25) 


Equations  (A-24)  and  (A-25)  reveal  that 


E[dS,^T(L)  ]  =  llA  -  21B  +  7C  -  E[d®,^®'^^T( L)  ] 


(A-26) 
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.*JF' 

2S+1 

Substituting  from  Eqs.  (A-14)  for  E[d^,  T(L) ]  we  find  the  relative  en¬ 
ergies  of  the  quartets  and  of  Jd^  configuration 


E(^G)  =  lOA  -  25B  4 

■  5C 

(A-27a) 

E(^F)  =  lOA  -  13B  4 

-  7C 

(A-27b) 

E(^D)  =  lOA  -  18B  4 

■  5C 

(A-27c) 

E(^P)  =  lOA  -  28B  4 

■  7C 

(A-27d) 

The  relative  energies  of  these  levels  for  the 

case  b/c 

=  .5 

and  in 

units  of  C 

can  be  computed  with  the  result  as 

f  0.5 

for 

4f 

(A-28a) 

E[dS,^®’*’^T(L)  ;  B/C  =  .5]  -  lOA 

C 

-  4.0 

for 

( A-28b) 

-  7.0 

for 

4p 

(A-28c) 

-  7.5 

for 

^G 

( A-28d) 

\  -17.5 

for 

(A-28e) 

A  careful  study  of  Fig.  4.4,  Ref.  l6,  indicates  that  only  at  b/c  <  .08 
do  the  quartets  lie  lower  compared  with  the  doublets;  but  even  at  this  region 
Hund's  Rule  does  not  hold  completely,  and  the  is  lower  than  ^F. 
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APPENDIX  B 


REDUCTION  OF  THE  TO  A  FUNCTION  OF  L,  S  AND  I 

To  derive  the  spin  Hamiltonian  Eq.  (55))  use  is  made  of  the  replacement 
theory 


{a'K||L.l||aK) 

(a'»n'  |l|aKm)  =  - =-^= -  (Km'  |L|Km) 

L(L+1) 

=  B(a' ,q:)(Kii' |L|Km)  . 


(B-1) 


where  I  can  represent  such  vectors  as  r  or  p.  Thus  we  can  replace  the  posi¬ 
tion  coordinates  with  angular  momentum  components  as  follows: 


X  Lx,  xy  =  ( LxIy+LyLx) /2,  x”^  I^  , 


(B-2) 


[Eq.  (9)  ]  is  already  independent  of  position  vector  and  assuming 
that  Russel-Saunders  coupling  holds  and  that  it  will  take  the 
form: 


=  K  L-S 


(B-3) 


where  K  is  a  constant,  L  is  total  angular  momentum  and  S  is  the  total  spin. 

[Eq.  (lO)  ]  is  one  of  the  Hamiltonians  containing  position  vectors 
r^"^  in  conjunction  with  _s^  and  _sJ.  To  replace  position  coordinates  with  the 
components  of  angular  momentum  we  expand  ^  as  follows: 


3(3  =L  ■  5(|s^+Tisj+^sJ)(|sJ+T)sJ+5sJ)  ] 


'3  ^  iJ 

1>J 

=  Rx  +  Ry  +  R2 


(B-4) 
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where  and  C  are  the  Cartesian  coordinates  of  r^J  and 

'Yl. 

fi*  ■  Z 

i>j 


i>J 


L(L+l)sisj  -  3[l|sxS^  +  Ly  VL^Iy)  ( Sysi 


i-i)  ] 


■*"Sz®y' 


Recalling  that 


2x4 


Pij^x  > 


hy  substituting  for  s^s^  in  Eq.  (B-5)  one  finds: 


TL 

Y  Ci(iJ)Pij  '  UL+1)S§_  -  3[l|S^  +  I  (I^V^zV^ 


1>J 


p[3  +1  (lyL2+L2ly)(SySz+SzSy)  -  L*^s|] 


^  -  r 


^  p[(3l|-L*^)S|  +  |(iIoc+2L^Iy)(-iSx+2SyS^)  ] 

-  r^  p[(3l|-L*^)Sx  +  2  (LxSx+6LzLySySz)  ]  +  ••• 

The  result  for  Ry  and  Rz  can  be  obtained  from  Rx  [Eq.  (B-8)  ]  and 
expressed  as 


R 


^  (L,S)  =  -M[(L.S)^  +  I  (L-S)  -  L(L+l)S(S+l)/3] 

Since  [Eq.  (ll) ]  is  already  given  in  terms  of  L  and  S,  ^ 
will  be  considered  now.  Each  of  the  three  parts  rf j( -I,  3r 
0jjg,(ri)si-^3  will  be  treated  separately  as  follows; 


(B-5) 

(B-6) 

V”} 

(B-8) 

,  can  be 
3 

(B-9) 

[Eq.  (12)] 

3 

. s^ri.I  and 
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(B-10) 


r 


3 

i 


Dn'(l.i 


kS-I) 


=  Nafg  ( IocSx^Ix"^^IxIocSx'*'^SxLzIz+L2;IzIocSx''"^SxIyIy+IyIyL2Sz) 

3  3 

+  2  y>y  ^  ->  z)  +  5  (x  -V  z,  z  >  x,y  -<■  y)  ] 

=  I  n"  [  ( IxSx+IySy+Iz^X  L^Ix+Iyly+L^I^) 

+  (  I^lx+Iyly+L2lz)  (  Lx^x''‘^®y'^^z®z)  ] 

=  i-i  +  L-S]  (B-11) 


^  r  -  I  r 


I  = 


n’”  S-I 


(B-12) 


and  finally  we  have 


^5  ^  N|a[rt-I  -  I  (L.S  L.I  +  L.I  L.S)  ]  +  L-I  -  k  S-iJ  (B-13) 

where  N  and  a  are  constants.  In  the  same  way,  [Eq.  (15)  ]  can  he  changed 
into  a  function  of  L,  S  and  I.  The  final  result  is 


6l 


j.,'  (L,s,i)  =  -  :}i^ 

=  K  L-S  -  M[(L-S)^  +  L-S/2  -  L*^S*V5 1 


+  p(L+2S)-H  - 

+  N{a[L*^S  1  -  3(L-S  L-I+L-I  L.S)/2]  +  L-I  - 


(B-14) 
K  S-I] 


+  p[(L-l)^  +  L-l/2  - 

7 


where  K,  M,  N,  p,  a,  p  and  y  are  independent  of  operators  L,  S,  and  I,  and 
^  [Eq.  (14)]  is  to  be  obtained  for  each  particular  crystalline  field.  ZnTe 

7 

and  CdTe  single  crystals  have  cubic  structure  with  each  Zn  or  Cd  surrounded 
by  four  tetrahedrally  positioned  Te  ions.  The  crystalline  field  appropriate 
to  this  structure  is  [see  Eq.  (43)]: 


'll  (r,9,^)  =  +  (5/l4)  V2[y44(0,^)  +  Y4.4(0,)!5)  ]) 


A4o(«4o/8(35z^-30r2z2+5r4)  +  sj’y/lk  a44(3r^ 


-4r4cos  20+r'^cos  40)  cos  4^] 


(B-15) 


=  A4o[o:4o/8  (33z'^-30r2z^+3r‘^)  +  V 5/l4  lO5r‘*sin'*0  cos  4^]. 


Recalling  that 


cos  4^  =  8  -  8  cos^^  +  1 

r4sin49  =  (x^+y^)^  =  -  2r^z^  + 
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one  finds  that 


r^sin^Q  cos  -  6x^y^ 


(B-16) 


Equations  (B-I5)  and  (B-I6)  reveal  that 


^  =  ^o/^Ci:4o(35z'^-30r2z2+3r4)  +  168O  n/5/1^  a44(x'*+y4-6x2y2) 


where 


Thus 


and 


«n,m  =  (-1) 


n 


(2n+l)(n-|m|)  I 
2(n+m) : 


1/2 


^44  =  (9/2(8)  :2:t)"/^ 


a40  =  (9/2rt)^/® 


(B-17) 


(B-18) 


(B-19) 


Substituting  for  Oji^ni  (B“1T)  one  finds: 


=  ^A4q/8  •/2rt)[(55z‘^-30r^z^+5r'*)  +  5(x^+y‘^-6x®y^)  ] 
=  (15A40/  N/2^(x‘^+y‘^+z'^  -  5r‘^/5)  . 


(B-20) 


The  next  step  is  to  transform  r'*  into  a  function  of  angular  momentum  com¬ 
ponents.  Consider  that 


r'^  =  x*^  +  y*  +  z*^  +  2(x^y2+y^z2+z^x2) 


(B-21) 


and 


x®y^  _ — [^g  ( LxIy+IyI^+I<x%%^'^%^^%''‘LxIyLxIy+IyLxIyI«) 
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(B-22) 


where 


LxIyLyLx  +  LyLxLxIy  =  2L|  +  2l|Ly  -  Lf  -  l|  (B-23) 

and 

LxLyLjc%  +  LyLxLyLx  =  Lf  +  2l|l^  -  -  l|  (B-24) 

Equations  (B-22),  (B-23),  and  (B-2i+)  reveal  that 

x2y2  IZ:D>i  (2I§l|+Ul2l|+3L|-2Lg-2L§)  =  |  (6l|l|+3Li-2l|-2l|)  (B-25) 

The  expressions  corresponding  to  and  can  be  obtained  from  Eq.  (B-25) 

by  proper  substitution  of  x,  y,  and  z.  The  sum  of  these  expressions  which  is 
needed  for  evaluating  Eq.  (B-21)  is; 

6(x2y24.y2z2+z2x2)  6(  L^Ly+IyLf+Lf  I^) 

+  3(Ll+l|+lf)  -  2(l|+l|+Lf) 

(B-26) 

-  2(i5+L^I^) 

=IC>  6(  l|l2+l2Li+L|l2)  -  (I|+I|+Li) 

Now  note  that 

1|I|  +  ifLf  +  L|I|  -  -  (L^+i4+L|)  -  (Lgl|+l|Lf+L|l|) 

or 

6(I|I|  +  I|L|  +  Lfl|)  -  3[L^^  -  (I^+I^+L4)]  (B-27) 

Equations  (B-20)  through  (B-27)  reveal  that 


61+ 


(B-28) 


^  (L)  =  B[]^  +  +  l|  -  I  (L*^  -  i  L*^)  ] 

'"7  ^  ^ 

=  B[I^  +  -  i  L*^(3L*^-1)  ] 

=  B[l^  +  i  L(L+1)(3L*^+3L-1)  ] 

Equation  (B-28)  is  similar  to  the  Bleaney's^®  results  for  this  case  and  B  here 

is  a  constant  obtained  as  a  result  of  transforming  9,^)  into  j/ (L)  . 

7  7 
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PART  IV 


SOLID  STATE  INSTRUMENTATION 


■by 

Glenn  G.  Sherwood* 


♦Captain,  USAF,  in  Civilian  Institutions  Program,  AFIT,  at  The  University 
of  Michigan. 


This  research  project  is  concerned  with  the  design  and  development  of 
a  solid  state  charged  particle  dE/dx  detector  to  replace  the  traditional 
gaseous  ion  chamher.  The  advantages  of  the  solid  state  counter  over  the  gas 
counter  are: 

(a)  Resolving  times  of  less  than  0.1  microsecond  due  to  short  collec¬ 
tion  times  in  the  semiconductor. 

(h)  Less  scattering  of  particles  out  of  the  "beam  "by  multiple  coulomb 
Interactions,  and  hence  greater  counting  efficiency. 

(c)  Ionization  efficiency  which  is  10  times  greater  than  that  of  the 
gas  counter,  and  hence  some  probable  improvement  in  resolution. 

(d)  Fast  rise  time,  which  permits  high  counting  rates  and  the  poss¬ 
ibility  of  using  pulse  shaping  techniques  to  discriminate  against  background 
gamma  radiation. 

Some  work  on  semiconductor  dE/dx  detectors  was  done  by  H.  E.  Wegner  of 
Los  Alamos  Laboratory  in  i960.  He  developed  a  diffused  junction  dE/dx  de¬ 
tector  with  a  thickness  of  50  microns  and  a  diameter  of  3/16-in.  As  Wegner 
reported  out  of  six  original  detectors  this  was  the  only  one  which  survived 
the  difficult  grinding  and  etching  process  to  give  a  workable  thin  detector. 
However,  this  one  was  enough  to  demonstrate  the  feasibility  and  advantages  of 
semiconductor  dE/dx  detectors. 

Our  research  project,  wherein  we  are  attempting  to  measure  the  spatially 
dependent  fast  neutron  fission  spectra  throTigh  several  types  of  semi-infinite 
media,  requires  a  proton  recoil  telescope  spectrometer  which  has  a  resolving 


3 


time  of  less  than  1  microsecond  and  is  relatively  insensitive  to  background 


neutrons  and  gammas.  Previous  telescopes  with  gaseous  counters  had  extreme 
limitations  imposed  by  the  large  resolving  times,  and  this  restricted  their 
performance  so  greatly  that  experiments  with  fission  sources  were  virtually 
impossible,  as  explained  by  Johnson. ^ 

Our  prescription  for  a  possible  method  of  measuring  the  spectra  is  shown 
in  Fig.  1.  Here  is  depicted  a  proton  recoil  telescope  with  a  semiconductor 
detector  preceded  by  a  semiconductor  dE/dx  detector.  The  proton  recoils  in 
energy  range  1  to  12  Mev  must  pass  the  dE/dx  detector  and  stop  in  the  thick 
E  detector.  Only  those  pulses  which  pass  through  both  detectors  and  are  hence 
in  time  coincidence  are  allowed  to  enter  the  pulse  height  analyzer.  This  re¬ 
duces  the  large  neutron  background  produced  by  the  Si^®  (n,p)  and  Si^®  (n,a) 
reactions  as  well  as  betas  from  the  walls  of  the  instrument. 

The  portion  of  this  telescope  spectrometer  project  covered  herein  in¬ 
volves  the  development  of  the  dE/dx  detector.  Three  factors  dictate  the 
size  of  the  thin  detector:  noise,  counting  statistics  and  thickness.  The 
energy  of  the  protons  require  a  detector  of  about  25  microns  so  that  1.5  Mev 
protons  will  pass  into  the  second  detector.  So  that  the  efficiency  can  be 
kept  in  the  order  of  magnitude  of  10"®  or  greater  (corresponding  to  an  energy 
resolution  of  15-20^)  for  reasonable  counting  statistics,  the  diameter  of 
the  finished  detector  should  be  as  close  to  2  cm  as  possible.  Noise  will 
limit  the  achievement  of  the  ultimate  resolution  and  efficiency  .  Theoret¬ 
ically,  the  noise  of  this  detector  should  be  about  85  Kev  if  predicted  by  the 
noise  analysis  of  Goulding.^  This  would  be  acceptable  because  the  12-Mev 
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Telescope  Housing 


Oscilloscope 

Printout 

Record 


Dual 

Beom 

Scope 


Fig.  1.  Block  diagram  of  telescope  spectrometer  and  necessary 
circuitry  for  coincidence  counting. 


proton  deposits  about  I50  Kev  in  25  microns  of  silicon. 

In  attempting  to  develop  high  quality  semiconductor  detectors  the  surface 
barrier  technique  appears  superior  to  the  diffused  junction  method  as  used  by 
Wegner  because: 

(a)  The  surface  barrier  method  does  not  require  excessive  heat  which 
degrades  the  performance  of  the  completed  diode. 

(b)  The  p-n  junction  diode  has  a  "window"  of  several  microns  which  de¬ 
grades  the  energy  of  slower  protons  without  collecting  charge. 
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.-lat 

With  these  advantages  in  mind  work  was  started  in  February  1962  to  de¬ 
velop  these  solid  state  devices. 

A  2000  ohm-cm,  500-microsecond,  vacuum  zone  floated  N-type  silicon,  2h 
ram  in  diameter,  is  waxed  into  a  circular  pyrex  boat  which  is  in  turn  waxed 
to  a  flat  aluminum  plate.  The  silicon  is  then  cut  into  discs  .055-in.  thick 
with  a  high  speed  diamond  impregnated  blade  on  a  Di  Met  No.  120BQ  cutoff 
machine.  The  silicon  disc  is  then  waxed  to  a  pyrex  optical  disc  5  in.  in 
diameter  ( see  Fig.  2) .  Three  small  blocks  of  silicon  are  waxed  120  degrees 
apart,  as  shown  in  the  figure,  and  it  is  the  purpose  of  these  blocks  to 
maintain  the  parallelism  between  the  face  of  the  silicon  crystal  and  the 
pyrex  blocks. 


Fig.  2.  The  silicon  wafer  is  shown  mounted  on  a  pyrex 
polishing  disc.  The  optical  glass  flat  and  the  iron 
flat  are  used  in  the  lapping  process. 

Two  types  of  optical  lapping  flats  are  used  in  this  process,  a  12-in. 
iron  flat  and  a  6-in.  glass  flat.  The  iron  flats  are  used  with  silicon  car¬ 
bide  abrasive  and  the  pyrex  flats  with  aluminum  oxide  abrasive.  The  final 
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polishing  is  done  with  Linde  A  alumina  on  Buehler  AB  Metcloth  with  the  Buehler 


polishing  wheel.  Both  iron  and  pyrex  flats  are  checked  periodically  to  ensure 
that  they  do  not  loose  their  flatness. 

The  silicon  ciystal  is  first  lapped  with  No.  400  silicon  carhlde  in  water 
on  the  metal  flat  to  remove  the  silicon  damaged  by  the  saw.  The  crystal  is 
then  lapped  with  No.  600  silicon  carbide  and  then  with  No.  1200  aluminum  oxide 
in  water.  The  polishing  blocks  are  measured  regularly  during  this  process 
with  a  micrometer  to  maintain  them  at  an  even  height  with  respect  to  the  back 
of  the  polishing  block.  In  this  way  the  removal  of  silicon  by  lapping  can  be 
controlled  so  that  the  crystal  is  ground  with  parallel  surfaces.  When  the 
pits  have  all  been  removed  by  the  No.  1200  abrasive,  as  observed  with  a  metal¬ 
lurgical  microscope  of  power  110,  the  crystal  is  transferred  to  the  glass 
optical  flat  and  lapped  with  No.  3200  aluminum  oxide.  The  crystal  is  then 
polished  to  a  mirror  finish  with  the  Linde  A  until  all  pits  and  scratches  are  re¬ 
moved.  At  this  point  the  surface  should  have  minimum  damage  (about  1  micron) 
but  should  still  exhibit  a  relatively  high  surface  recombination  velocity. 

The  crystal  is  then  turned  over  on  the  3-in.  pyrex  flat,  polished  side 
down,  and  the  reverse  side  of  the  crystal  is  lapped  and  polished  in  the  manner 
explained  above.  The  crystal  is  polished  down  to  45-micron  thickness,  with 
particular  care  being  taken  to  insure  that  the  crystal  is  planar.  This  can 
be  controlled  by  measuring  the  edges  with  a  microscope. 

The  crystal  is  removed  from  polishing  by  gentle  heating.  It  is  washed 
in  successive  baths  of  trichloroethylene,  alcohol,  and  distilled  water,  dried 
in  a  clean  glove  box,  and  finally  cemented  to  a  lavlte  type  A  ring  with 
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Araldite  No.  951  amine-type  epoxy  (see  Fig.  2). 

The  back  of  the  crystal  is  protected  by  sax  while  the  entire  crystal 
and  holder  is  immersed  in  an  etch  bath  (nitric  acid,  hydrofluoric  acid,  and 
acetic  acid  in  volume  proportions  2:1:1)  for  two  minutes.  After  washing  in 
distilled  water,  the  wax  is  removed  from  the  back.  After  24  hours,  75  micro- 
grams/cm^  of  gold  is  evaporated  on  the  front  and  the  same  quantity  of  aluminum 
is  evaporated  on  the  back.  The  semiconductor  dE/dx  detector  is  allowed  to 
age  for  several  days  and  is  then  conpleted. 

Experimentation  with  several  etches  showed  that  the  most  controlled  etch 
producing  a  chemical  polish  on  the  l/2-micron  polished  surface  of  the  silicon 
was  of  the  above  noted  composition,  maintained  at  0°C  for  about  2  minutes. 

This  procedure  removes  evenly  about  10  microns,  leaving  a  polished  surface 
which  appears  to  produce  the  surface  barrier  required. 

The  important  aspect  of  this  technique  is  that  we  have  apparently  been 
able  to  control  the  etch  rate  without  any  preferential  action  on  the  part  of 
the  etch  solution.  Using  this  procedure  on  the  face  of  the  silicon  crystal 
apparently  provides  the  density  of  fast  states  necessary  for  the  inversion 
layer.  With  evaporation  of  a  gold  layer,  the  rectifying  Junction  appears  to 
be  complete. 

We  are  still  experimenting  with  back  contacts  to  achieve  a  stable  or 
ohmic  contact  which  is  free  from  injection.  We  have  only  partially  achieved 
this  because  the  back  contacts  still  provide  some  rectifying  action.  The 
tests  to  show  the  amount  of  minority  carrier  injection  are  still  inconclusive. 
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The  latest  dE/dx  detectors  have  progressed  to  the  point  where  we  are 

now  producing  detectors  J-cm^  in  finished  area,  and  30  microns  thick,  with 

210 

stable  characteristics  over  periods  of  weeks.  Resolution  for  Po  is  5.8^ 

( see  Fig.  3)  the  peak  noise  extends  to  channels  corresponding  to  about 

280  Kev.  Improvement  in  the  back  contact  should  help  to  decrease  the  noise 

to  the  acceptable  level  of  I50  Kev.  Figure  4  shows  the  pulse  shape  versus 
PI  n 

bias  for  Po  alphas. 

The  dE/dx  detector  in  conjunction  with  a  thick  detector  is  now  operating 
in  a  proton  recoil  test  spectrometer  with  coincidence  equipment.  This  tele¬ 
scope  spectrometer,  shown  in  Fig.  5,  is  being  used  to  help  find  those  sources 
of  noise  which  can  be  removed.  (Since  tests  have  just  started  there  are  no 
conclusive  results  as  yet.)  The  mechanism  shown  in  Fig.  5  is  the  telescope, 
which  operates  in  a  vacuum  housing.  The  front  wheel  holds  the  polyethylene 
radiators,  which  are  of  different  thickness.  Behind  the  front  wheel  is  the 
dE/dx  detector,  and  behind  this  detector  is  the  E  detector.  The  rear  of  the 
device  holds  the  motor  for  positioning  the  different  radiators  corresponding 
to  various  parts  of  the  neutron  energy  spectrum  being  sampled. 

The  instrument  has  been  Initially  tested  with  a  weak  Pu-Be  source.  Fig¬ 
ure  6  shows  the  response  of  the  dE/dx  detector  and  E  detector  to  recoil  pro¬ 
tons.  The  picture  is  poor  because  the  count  rate  is  low  and  a  30-45  second 
exposure  was  necessary.  The  coincidence  counting  rate  is  about  1  per  minute. 
In  January  the  instrument  will  be  tested  with  the  Michigan  Van  de  Graaff. 

Initial  results  indicate  at  least  partial  success  with  the  new  semicon¬ 
ductor  dE/dx  detectors  and  some  success  with  the  telescope.  As  far  as  we  can 
determine  no  other  group  has  achieved  success  as  yet  in  producing  these  solid 
state  counters. 
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CHANNEL  NUMBER 


Fig.  4.  Pulse  shape  of  large  area  dE/dx  detector  to 
alphas  vs.  bias  across  the  detector  and  22 -meg  load  resistor. 
Voltage  is  0,  25,  75,  and  150  volts;  the  leaJcage  current  is 
over  1  microampere. 


Pig.  5«  The  telescope  spectrometer.  The  detector  nearest 
to  the  proton  radiator  wheel  is  the  dE/dx  detector.  The 
second  detector  is  a  semi-conductor  E  detector. 
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Fig.  6.  The  top  line  shows  the  response  of  the 
dE/dx  detector  and  the  bottom  the  E  detector  to 
proton  recoils  using  a  Pn-Be  source. 
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